^{1}

^{*}

^{2}

^{*}

According to the production-specific SPQR axiom, which addresses the observance of business requirements of speed, punctuality, quality, return [

Just in time (JIT) production systems are mainly taktrate (TR)-governed manufacturing systems with SD[TR] = 0 where SD denoting the standard deviation and TR representing the pace of ordering. The necessary and sufficient conditions for a TR-governed manufacturing system to observe OTD is given by the Equation system 1 [

{ ∀ i : inf { E R i } ≥ T R Z : P L T Z ≤ E D T (1)

The first condition for OTD requires that the capacity is sufficient—capacity given by the exit rate (ER) which is determined by the smallest ER_{i} of each process step. The ER has to be large enough to process the arriving orders expressed by the taktrate (TR). The second condition to be observed requires that the process lead-time (PLT) is shorter than the customer’s expected delivery time (EDT) confirmed at the moment of the ordering. Having a TR-governed production system with a deterministic call-off characteristic of production orders within a defined timespan, we can assume an even spreading of incoming orders within this timespan. With spreading we intend here the time-dimension of distribution.

For non-deterministic order-based JIT production systems, or more generally, for batch& queue (B&Q) production systems based on a push manufacturing principle—which is a manufacturing system modus operandi usually applied in the case of erratic order entry—Equation (1) has to be modified. Indeed, the order arrival rate is stochastic and we should not talk any more about a TR, but about an order rate (OR) with the characteristic of SD[OR] > 0. To observe OTD in this case, the more general Equation system 2 is valid

for : S D [ O R ] > 0 { ∀ i : inf { E R i } > E [ O R ] Z : B W T Z + P L T Z ≤ E D T (2)

The variable OR will sometimes exceed ER, and therefore according to the Theorem of WIP build-up a work in progress (WIP) stock at the process entry point Z which we will name backlog (BL). Usually the production order input rate (IR) to the shopfloor is reduced to match the ER, i.e. IR = ER. The production system to be modelled is schematically shown in

The customer visible time (CVT) will therefore comprise the PLT plus the WIP-generated backlog waiting time (BWT) which has to be shorter than EDT; for simplicity reason, we are not considering the shipping time, which is not an

issue of manufacturing. Equation (2) is directly derived from the following important theorem, the Theorem of General Production Requirements or OTD Theorem, which sets the necessary conditions for OTD [

Theorem of General Production Requirements (or OTD Theorem): The necessary and sufficient conditions to supply a customer with OTD, i.e. with the right quantity at the right time, is that first the capacity requirement and second the lead time requirement have to be satisfied simultaneously, independent of the applied transfer principle, i.e. SPF or B&Q. The capacity requirement is given by the Corollary of Weak WIP Stationarity and the lead time requirement necessitates that MLT plus BWT is shorter than EDT.

As mentioned in the OTD theorem, it would be more appropriate to talk about MLT (Manufacturing Lead Time) instead of PLT, which indicates the timespan one single piece needs to transit the process. The MLT comprises the lead time to produce the entire batch B(n) with the transfer principle of n pieces, i.e. when the first piece enters the process and the last piece of the batch leaves the process. In Western production systems, with its dominant batch-culture and B&Q production systems making the lion’s part, the difference between MLT and PLT can be considerably big. The application of the MLT metric will be mandatory in multi-product manufacturing systems where order-batches of different products are produced and have to be delivered again on time as whole batches to the customer. In this introductory paper of conceptualization regarding post-optimality for OTD we will refrain from considering batches and will therefore view each order as a single product to be delivered.

These are the general, theorem-based requirements to observe OTD. However, such as Equation (1) for the TR, also Equation (2) for the OR applies the discrete Poissonian view of average occurrences per timespan represented by λ. If λ is an average, the same average can be the result of different distributions of the phenomenon to be modeled. The two extreme distributions are on the one hand the ideal situation of the even spreading over time of occurrences with constant inter-arrival time (TR-like) and on the other hand the bulk occurrence, i.e. the verification of concentrated arrivals. In this context we rather prefer the word spreading instead of distribution, indicating the focus on the random sequence of arrival orders (or the sequence of changing patterns which would lead to a non-ergodic process) in the time-dimension t of the distribution τ of the stochastic variable x(t). However, we use the appropriate word distribution to characterize the sequence-independent view of arrival orders represented by the aggregated shape of the inter-arrival time τ = x − 1 ( t ) of the random variable. If the process is ergodic, it is not necessary to make this difference because the process is time-invariant. It seems to be clear that in the case of an ideal situation of even arrival sequence of orders, Equation (1) can be observed theoretically depending in reality only on the characteristics of underlying manufacturing systems to transform efficiently “order arrivals” into “product deliveries”. However, if the arrival spreading of orders is aggregated (bulk arrivals), the instantly needed capacity is not enough to process all the orders at the same instant forcing them to be queued and compromising eventually the observance of OTD. This is true even though the average condition of Equation (2) may be observed. Indeed, whereas usually the constant ER for a mono-product manufacturing system always presents an even deterministic spreading over time, this might not be the case for the OR. In the case of non-deterministic arrival rate the question now is, in which post-optimal range of a non-evenly time-distributed OR, Equation system 2 will be observed.

The treated topic is complex and until presently it has not yet been investigated although queuing theory application is wide-spread. In fact, the findings have far-reaching implications in production economics and manufacturing theory as well as optimization under uncertainty. Therefore, this paper does only have introductory character. The complex character of the topic requires a swing back to basic theory. This entails rather a didactic connotation of the presentation and might be seen as an unnecessary step in scientific papers. However, in this context for sure it increases understanding of the whole topic allowing to reaching a larger readership.

In the following, we will use for the solution finding a new production mathematics’ reasoning based on a Cartesian scientific approach [

Let us have a closer look at two important production topics, on the one hand the production capacity and on the other hand the manufacturing elasticity. We prefer to call it manufacturing elasticity and not production elasticity in order not to confuse it with the economic output elasticity. In this paper, we are explicitly not interested in the topic of production flexibility, i.e. the capability to transform different type of products resulting in a high mix variability. Not because this topic is not of importance, on contrary, however, because it adds a further complexity to the problem statement, which is a variable ER instead of a deterministic ER with a mix-specific distribution. Indeed, the concurrent optimization of two target functions, such as elasticity and flexibility, is a problem of Pareto efficiency [

Production capacity defines the amount of products, or generally, the volume of produced items, which a certain production infrastructure can transform from raw materials or components into intermediate or finished goods. The expressions throughput or exit rate (ER) define the concept of capacity figuratively. The concept of production capacity is directly linked to the concept of bottleneck, modelled in Equation (3)

{ C T b = sup { C T i } E R = E R b = 1 C T b (3)

which is a mathematical representation directly derived from the following Theorem of Throughput [

Theorem of Throughput (or Bottleneck Theorem): Given is a sequence of production steps with each process step having a deterministic, but different cycle time CT. The maximum throughput, i.e. the maximum exit rate ER, of a process is given by the slowest process step, i.e. the process step with the longest cycle time CT; this process step is called bottleneck.

The ER is a performance concept, defined as number of pieces per time unit and represents the productivity of the process constituting the nominal specific capacity. The absolute production volume V is given by the timespan t the ER is applied, i.e. V = E R ⋅ t . The production capacity is therefore also a concept of the applied work shift regime.

Manufacturing elasticity, however, deals with the cost of a variable production regime subject to load changes for the quantity to be produced. Usually, production capacity is defined for an optimal volume range, production minimum given by the long-term minimal economic sustainable production level V_{BEP}, i.e. by the break-even point (BEP). In this paper we are not primarily interested in the BEP, but rather in the concept of insensitivity to cost change for various production level regimes. The BEP is a function of fixed cost K_{fix} of investment depreciation as well as fixed operating cost (Equation (4))

V BEP : = K f i x p − k v ( V ) (4)

where p is the price and k_{v}(V) represents the variable unitary specific cost. Therefore m ( V ) = p − k v ( V ) is the margin contribution per unit sold. Please note that the variable cost show in reality a variable unitary characteristic with cubical exponent in function of increasing production volume V due to inefficiencies at low volume and increased expenses for additional shifts. Moreover, the production volume where the variable cost k_{v} are minimum is given by

V min { k v } : = { V | min { k v ( V ) } } → { V min { k v } | d d V k v ( V ) = 0 }

This leads to a non-linear model with an optimal finite operational production level V_{opt} above the BEP and not in the infinite like a linear model would suggest. The optimal production level is the associated production volume where the marginal price decrease equals the marginal variable cost increase

V o p t : = { V | d d V p ( V ) = d d V k v ( V ) }

or

V o p t : = { V | max { ( p ( V ) − k v ( V ) ) ⋅ V } } → { V o p t | d d V ( m ( V ) ⋅ V ) = 0 }

In addition, it is evident from the equation that the optimal production regime V_{opt} is not at the volume level where k_{v}(V) is minimal, but where the margin m(V) is maximal, i.e. marginal cost equal marginal price (

The “art” of maximizing manufacturing elasticity is therefore not reduced to the problem of minimizing investment in order to lower the BEP, but to have an extended range of volume with reduced k_{v}(V). How industry 4.0-type production system, allowing batch size one (one-offs), will solve this contradiction, is left out in this paper. However, this is very important. Because to define correctly elasticity of production we have to distinguish two different elasticities: we can call it improperly the average elasticity based on a full costing model and the classical is the marginal elasticity based only on variable cost. Please note that

the rarely precisely defined word elasticity in production theory of manufacturing, if ever defined, is intended of having low sensitivity of cost change to volume change ε = min { ∂ k / ∂ V } . However, nota bene, this corresponds in economics to inelasticity. In production theory elasticity has a different meaning. Production systems with high elasticity allow the usage of the equipment over a large production volume range without significant cost changes. Elasticity is linked to the kind of invested equipment and work shift model. Nevertheless, automation generally comprises high investment increasing BEP and therefore lowering, i.e. worsening the average elasticity, presenting a high average elasticity only in the high volume range (

Theorem of Optimal Production Range (or Limited Elasticity Theorem)

Given is a production system characterized by fixed cost K_{fix} and non-constant specific variable cost k_{v}(V) in function of production volume. The marginal elasticity is the insensitivity of variable cost in presence of changes in production volume. Due to the variability of the specific variable cost, the optimal manufacturing elasticity is confined within a limited range around minimal cost volume V min { k v } and with economic meaningful extension towards the volume around the optimal production level V_{opt}.

Therefore, this theorem reveals that the exploitation of optimal manufacturing elasticity is not only a function of cost, but also of revenue. It is directly linked to the Postulate of Incompatibility (or Flexibility-Elasticity Contradiction Postulate) enounced in [

We are used to apply statistical parameters, especially such as the mean value, everywhere we have to describe synthetically the characteristic of a physical or social phenomenon. The mean is the average or better the expected value of a stochastic data set of a random variable X, represented with E[X]. Moreover, many think that the mean is a good representation of this random variable X and don’t further care about the dispersion of the data around the mean described as the variance Var[X] or the standard deviation SD[X]. However, the “mean” does not exist physically, but is basically an artificial value—a “fiction” if you want—derived from a dataset of variable X with a mathematical description (the first statistical moment). In fact, a random variable is usually represented by its real distribution. Generally, statistical parameters are not sufficient to define a distribution unequivocally [

Let us suppose a random process O R = X ( τ , t ) with an average arrival rate of call-offs (representing potentially immediate executable manufacturing orders) E[OR] of λ = 12 orders/day during a certain timeframe t and that the supplier has agreed with the customers to respect a commonly accepted EDT. These 12 orders per day might already be an average scale down from 60 weekly orders. We are used to transform weeks in days, days in hours, or even hours in minutes for high performance lean-takted lines or paced manufacturing cells, to bring it to shopfloor cycle time (CT) dimensions shown in Equation (5) which looks logic and apparently not questionable.

λ = 12 / day = 1.5 / h (5)

with 8 working hours/day

The mean is an interesting parameter because it is a so-called non-biased estimator of a sample statistics for the population with a reduced SE of the mean

given by σ E [ X ] = σ X n . Nevertheless, Equation (5) makes implicitly an important simplification of an even distribution of the order arrivals.

If we are bringing the average consideration of 12 orders/day to a proportional value of 1.5 order/h, we imply an even uniform spreading of arrivals over the time span of one day, which could be represented again by a TR. Then this TR can be approximated to a “deterministic” ergodic process. With an OR of 1.5/h, most probably, it is potentially possible to fulfill the OTD theorem of Equation 2 for all orders. It will, however, most probably not be possible to fulfill it with the mean logic of considering E[X] = 12/day, because it depends on the underlying order arrival distribution X. Indeed, for a potential bulk arrival of all 12 orders within one hour, only for some of the first orders of the bulk OTD may be respected, but not for all orders. This leads directly to the apparent paradoxical contradiction

λ = 12 / day ≠ 1.5 / h

Therefore, we can generalize and let evolve the attention catching title of this section “Paradoxon of the mean” into the more appropriate statement “paradoxon of proportions”. While an infinitesimal approach applied to the distribution would not rise any concern, the so often applied discrete view, using for simplification an average approach based on the mean, however, may become a concern—in certain cases.

Indeed, we have to analyze the arrival process by taking into account its real distribution. Therefore, we will look at the OTD performance of a defined production system under different arrival random variables. For this purpose, we define the following queuing system taking alternative arrival processes into consideration. We consider a M = Markovian (i.e. exponential), U = uniform (even), N = normal distribution, and D = deterministic according to Dirac’s δ-function, all distributions with the same mean inter-arrival time E[τ(OR)] = 1/λ where τ(OR) denotes the inter-arrival time of the order rate distribution. In our example 1/λ = 0.66 h. Furthermore, let us take a deterministic D characterized server process, with only a single server, which in our case is a machine or workstation (no parallelization of processing capacity) and unlimited buffer capacity b, queuing system formalized as

{ ( X i / D / 1 / ∞ ) X i : = M , U , N , D

If the machine s is not a single stage server but a process, we get the case of a manufacturing process. Let us compare E_{X}[τ(OR)] = 1/λ of the four arrival rate distributions X to the deterministic CT of the workstation or the CT_{b} of the bottleneck of the process, which determines the ER of the process, supposing to be CT_{b} =1/λ. The distributions M, U, N are shown in

The probability of symmetric distributions P_{X} (the third statistical moment is µ_{3} = 0) for shorter or longer arrival intervals than the cycle time of the bottleneck is intrinsically given by the definition of the mean (and median) and equals 0.5. The random variable X stands for the four specific distributions and P_{X} is the cumulative distribution function (i.e. it is the integral of the probability density function p_{X} of the arrival intervals).

Nota bene, the OR cannot really be influenced and hardly a company decides to refuse an incoming order. Therefore, to maintain an efficient and lean manufacturing process, the input rate (IR) of the order release to the shopfloor will be adapted down to the existing ER capacity, with the consequence that orders exceeding instant capacity will become backlog (BL) and therefore experience a BWT. This has not to be a concern—it is just a natural consequence of the capacity limit of the production line. However, OTD might be compromised. Nevertheless, this BL could then be rescheduled according to customer-specific

urgencies with different scheduling principles, for instance LIFO instead of FIFO or a preferential “VIP” scheduling might be applied to observe OTD.

Now, transforming the Poissonian rate-view of the first condition of Equation 2 into the Markovian time-view becomes

inf { E R i } : = sup { C T i } = C T b < 1 E [ O R ] = 1 λ

We can therefore assume that for inter arrival times τ = 1/x longer than CT_{b} (i.e. 1 / x > 1 / λ = C T b ) the first condition is observed, being possible to process immediately arriving orders if no WIP exists. On the other side, for the contrary with short arrivals (i.e. 1 / x < 1 / λ = C T b ), the capacity (i.e. ER) will not be enough to process the orders immediately. Indeed, according to the Theorem of WIP [

For the exponential distribution the probability P_{M}(τ) of arrival times up to CT_{b} and above CT_{b} is given by the Equation system 6

{ P M ( 0 ≤ τ ≤ C T b ) = 1 − e − λ ⋅ C T b P M ( C T b < τ ≤ ∞ ) = e − λ ⋅ C T b (6)

The Markov distribution is very interesting, because it is characterized by only one parameter λ. This reality-matching probability density function is often used to model queuing systems, covering partially also the bulk arrival characteristics. It has the advantage of an easy to calculate integral of the cumulative distribution function, not necessitating statistical tables, which makes it the preferred stochastic function to model queuing systems. Without showing the calculus, assuming C T b = E [ τ ] = 1 / λ , where τ =1/x then

{ P M ( 0 ≤ τ ≤ C T b ) ≈ 0.63 P M ( C T b < τ ≤ ∞ ) ≈ 0.37 (7)

meaning the probability of short inter-arrivals is higher than the long inter-arrivals above the mean, asymmetry which is typical of skewed probability density functions. Again, without showing the calculus, the median results to be

E ˜ [ τ ] = 1 λ ln 2 < E [ τ ] = 1 λ

what demonstrates that the median is below the mean being

{ P M ( 0 ≤ τ ≤ E ˜ [ τ ] ) = 0.5 P M ( E ˜ [ τ ] < τ ≤ ∞ ) = 0.5

Remember the median corresponds to the 50% percentile. This is an indication that the mean arrivals should not be applied for the first condition of Equation (2); the median of arrivals would be more appropriate, because the number of arrivals is relevant. If the order arrivals are not single orders but each order has an own batch-size of its own, the situation changes further. This situation has to be taken into consideration in the case of a multi-product manufacturing system. In this case here, we can neglect this further complication which would require an adaptation of the OR and ER.

For the uniform distribution, the probability P_{U}(τ) of arrival times up to CT_{b} and above (assuming τ_{max} = 2CT_{b}) is shown in Equation (8)

{ P U ( 0 ≤ τ ≤ C T b ) = 1 − 1 τ max ⋅ ( τ max − C T b ) = 0.5 P U ( C T b < τ ≤ τ max ) = 1 τ max ⋅ ( τ max − C T b ) = 0.5 (8)

The uniform distribution is simple, but doesn’t reflect well the reality of arrival rates. Nevertheless, through its simplicity, it is very easy to analyze the post-optimal conditions to observe OTD.

The normal cumulative distribution function (we do not apply the better suited lognormal distribution for simplicity reason, although being conscious of non-existing negative arrival times) is shown in Equation (9). If we assume CT_{b} = 1/λ we obtain

{ P N ( − ∞ ≤ τ ≤ C T b ) = 1 − Φ ( z ) = 0.5 P N ( C T b < τ ≤ + ∞ ) = Φ ( z ) = 0.5 (9)

where Φ(z), the standardized integral of the probability density function p_{N}(τ), is a non-elementary function difficult to calculate

Φ ( z ) : = { 1 2 π ∫ 0 z e − 1 2 τ 2 d τ | z = τ − μ σ }

The symmetric normal distribution is made-up of two independent parameters, the mean and the variance. We will conceive the variance as small as necessary in order not to get negative arrivals. This parametric analysis of variance leads at the limit to Dirac’s δ-function

lim σ → 0 N ( μ , σ 2 ) = δ ( τ )

which we model as deterministic arrivals in the following with S D X [ τ ] = 0 . For deterministic arrivals following Dirac’s δ-function, the probability density function is summarized in Equation (10)

{ p D ( 0 ≤ τ < C T b ) = 0 p D ( C T b ) = ∞ p D ( C T b < τ ≤ + ∞ ) = 0 (10)

with P D ( C T b = τ < ∞ ) = 1 . Being X in this case a deterministic arrival process, we can substitute the concept of OR with that of TR = λ for which the Equation (1) is valid building-up no backlog with no BWT. Assuming no saturation effects of the process, the production process will therefore be able to process all arriving orders on-time. This is the ideal case for a deterministic TR-based JIT-production system. As this case is deterministic, no post-optimality analysis is necessary, because the TR has to be forcedly below the nominal ER to have a sustainable solution. However, we will analyze the post-optimality conditions for a non-ergodic type process within a restricted time-period for which the deterministic OR is above the ER.

To analyze the post-optimality conditions to respect OTD we have to distinguish different cases. The simplest case is an arrival order process of a single product and stage-changing deterministic OR (“deterministically non-ergodic” process). A more general case, always based on a mono-product manufacturing line, assumes the OR as a stochastic variable of ergodic process characteristic. For these two cases, the Second Corollary to the Theorem of Throughput (Corollary of Bottleneck Time-invariance) applies which states that the ER for a mono-product manufacturing line is time-invariant [

In our case, the post-optimality conditions are determined by the deviation of the OR from the nominal capacity of the production system, which, according to the Theorem of Throughput, is given for a mono-product manufacturing line by the deterministic cycle time CT_{b} at the bottleneck. According to the OTD theorem and Equation system 2, the solution space is also limited by EDT. Whereas OR < ER is not critical, OR > ER might become a concern. The post-optimality question is: what is the sup{OR} and the maximum timespan sup{Δt} of a temporarily higher order-entry regime while still observing the second condition of Equation system 2. During the Δt the OR behaves like a TR. We are expressly writing sup{OR} and not using the diction max{OR}, because we are not intending to maximize a target function, i.e. solving an optimization problem, but we are interested in the largest admissible value of the random-deterministic variable OR.

In this case of non-ergodicity we will assume for OR(t) and especially for the sup{OR} deterministic values over the timespan Δt. To develop the post-optimality conditions for observing OTD we have to take Equation system 2 and solve it in view to determine the solution space for the OR. Please note that apart from the SPQR axiom, if the Theorem of General Production Requirements would not have been enounced, the problem of post-optimality of order entry cannot be defined, hence the fundamental importance of this OTD theorem for production systems. Notice that the compliance of OTD refers directly to the punctuality axiom.

For developing the post-optimality conditions, we can rewrite the second condition of Equation (2) as

B W T ≤ E D T − P L T = Δ T

where ΔT allows to let aside the real values of EDT and PLT. In any case it has to be BWT < ΔT. The backlog waiting time (BWT) can be modelled by setting OR = x and ER = λ

B W T ≤ B L E R = B L ⋅ C T b = B L ⋅ 1 λ = ( x − λ ) Δ t ⋅ 1 λ

and where the backlog (BL) has been computed as

B L = ∫ t 1 t 2 ( O R − E R ) d t = ( x − λ ) ⋅ t | t 1 t 2 = ( x − λ ) ⋅ Δ t

and where Δt is the timespan of duration the higher order regime (x > λ) is active. The resulting post-optimality condition to determine the admissible range of data for x and Δt is shown in Equation (11)

( x ( t ) λ − 1 ) ⋅ Δ t ≤ f ( Δ T ) (11)

Equation (11) is a parametric relation of implicit type, such as f ( x , y , k ) = 0 , with the variables x(t) and Δt in parametric function of ΔT, i.e. the difference between the customer’s expected delivery time and the process lead time. Therefore, given EDT and being ΔT = EDT-PLT, the faster the throughput speed PLT is, the production systems can process a higher order rate during a certain timespan while respecting OTD. This intuitive, but now finally derived and explicitly formulated conclusion is very important, because it allows for non-TR governed production systems to play with the backlog. That means high order entry can be managed by buffering temporarily instant-capacity exceeding orders as BL and reschedule them accordingly to observe OTD. These findings are compliant to the Lemma of Flexible Scheduling Principle [

From

observe OTD. If ( x λ − 1 ) > 0 then BWT will increase, else decrease (

right side). The virtual elasticity can be exploited until the BWT has reached ΔT, which represents the maximum intrinsic virtual capacity reserve of fast production lines. When this virtual production capacity has been exploited, backlog has to diminish with an OR below ER to be exploited again. The maximum reduction

rate of BL is when lim x → 0 ( x λ − 1 ) = − 1 representing the case of no-order entry.

These time-bound discrete values of x refer to deterministic order entry, an OR situation which corresponds seldom to reality. Such a changing deterministic process behavior with E D [ τ 1 , Δ t ] ≠ E D [ τ 2 , Δ t ] and S D D [ τ , Δ t ] = 0 within defined timespans Δt can be assimilated to a non-ergodic process.

If the order rate x is not a deterministic value, as assumed before, but more realistically a random variable X, the probability to experience a higher order rate x than λ depends on the probability density function and is (written in the Markovian time-dimension of inter-arrival time τ)

P X ( 1 x < τ < 1 λ ) = ∫ 1 / x 1 / λ p X ( λ , τ ) d τ

and depends not only on the increased OR, but also on the type of probability density function p_{X}(τ). Please note that here τ stands for the inter-arrival time. Do not confuse it with the timespan Δt during which a higher OR is present. Let us take as an example an OR higher up to 30% respect to λ = 1.5/h resulting in extremis x = 1.95 orders/h. This corresponds for the Markovian time view to a shorter inter-arrival time of 1/x = 0.51 h compared to the CT_{b} = 1/λ = 0.66 h resulting into a 23% shorter inter-arrival time.

The computation reveals for a Markovian exponential process a probability of 9.4%, for an even distribution a probability of 11.4%, and for a normal distribution a probability of 25.2%. The big difference between 9.4% and 25.2% shows that the assumed underlying distribution of arrival rate OR(τ) is of utmost importance. Therefore, for non-deterministic arrival rates, i.e. for a random variable X, Equation (11) has to be rewritten resulting in Equation (12)

∫ t 1 t 2 ( x ( τ , t ) λ − 1 ) ⋅ d t ≤ f ( Δ T ) (12)

Equation (12) models the generic condition for post-optimality of a stochastic process to observe OTD. If the stochastic variable x ( τ , t ) becomes time-invariant x ( τ ) , i.e. E[X] = const and SD[X] = const, then the process is ergodic. The computational solution of this integral equation for different OR distribution x ( τ ) is left to the interested reader.

In Equation (11) we have assumed a deterministic OR temporarily higher than the nominal capacity calculating the timespan during which a higher OR regime might be “digested” by the manufacturing system. However, this is not the typical case of a random variable. Indeed, how a random variable is defined, also lower OR than the mean E[X], i.e. the assumed capacity given here by λ, have a certain probability to occur reducing the BL-WIP. Equation (12) has a general validity covering also this real situation of a random variable. Equation (12) gives the general post-optimality condition for stochastic OR in presence of a fixed capacity characterized by the ER in parametric function of faster PLT than the requested EDT. This difference ΔT between PLT and EDT can be considered to be the built-in virtual elastic capacity of a production system.

We can therefore enounce the following Corollary of Post-Optimality (or Corollary of Virtual Elasticity). This new corollary belongs to the mathematical category of problem existence and in this case, not giving the solution, it is of non-constructive type for which Equation (12) is a non-constructive proof for Equation (2) (OTD theorem):

First Corollary to the Theorem of General Production Requirements (Corollary of Post-Optimality or Virtual Elasticity)

Given is a production system with a deterministic exit rate ER and a stable process lead-time PLT. The market is characterized by a stochastic order arrival rate OR each single order with expected delivery time EDT. The post-optimality range to comply to on-time-deliveries OTD depends on the difference ΔT between EDT and PLT as well as from the distribution of OR. This virtually increased capacity is called virtual elasticity of a production system.

This corollary seems to be in contradiction to the OTD theorem itself because according to the corollary it can temporarily also result OR > ER while still respecting OTD. This has not to be a contradiction because the OTD theorem is generally valid and over a long time-span whereas this corollary is valid for a limited time-span (post-optimality condition) according to Equation (12) (indeed, by increasing virtually production capacity) and constitute generally only a weak solution for OTD. However, we can therefore take a practical recommendation, synthesized in the following

Lemma to the Corollary of Post-Optimality (Lemma of Increasing Speed)

Based on ΔT = EDT − PLT, accelerating PLT of the production process increases the manufacturing-inherent reserve ΔT of production capacity, i.e. the virtual elasticity, in order to observe OTD for stochastic OR.

Therefore, it is always recommendable to speed-up processes and to reduce PLT. This is in perfect concordance with the Lemma to the Theorem of Generalized Lead Time (Lemma of SPF Desirability) [

The new coined term virtual elasticity is opposed to costly conceived physical elasticity of production equipment, which is depending on the BEP. Elasticity, i.e. the non-sensitivity of production cost to load-changes, will become very important for so-called Industry 4.0 type production systems allowing batch size one [_{b} i.e. the cycle time at the bottleneck.

Furthermore, symmetric probability density functions, such as uniform (U) or normal (N) distributions, show the same cumulative distribution probability on the left and on the right side of the mean or median, levelling out in the average higher and lower OR than the ER. This is not the case for a skewed probability density such as the exponential Markovian (M) density function, i.e.

∫ 0 1 / λ p M ( τ ) d τ ≠ ∫ 1 / λ ∞ p M ( τ ) d τ (13)

where E [ τ ] = 1 / λ , and referring to the case before with λ = 1.5

P M ( 0 < τ < 0.66 ) > P M ( 0.66 < τ < ∞ ) i.e. 0.63 > 0.37 (QED)

which confirms the statement of Equation (7). Due to the asymmetric probability of the probability function with regard to the mean (Equation (13)), this would imply an increasing backlog. Excluding seasonality effects and economic growth, this suggests that not the mean, but the median has to be taken as reference for capacity considerations.

Now, the asymmetric extension of inter-arrival time of probability P_{M} (0.19 < x < 0.45) with higher density and P_{M} (0.45 < x < 0.89) with lower density have the same probabilities to occur which is 0.25 corresponding to the definition of quartile.

These quartile/median figures are quite different from the figures computed with Equation (13) referring to the mean. These differences do not exist in symmetric distributions.

Conclusion: If symmetric distributions show random arrivals, they represent stochastically a strong solution for OTD depending only on Equation (2a). If they show an enhanced but non-random behavior above the mean over a certain timespan, they can only have a weak solution for OTD. Asymmetric distributions always represent only weak solutions. A strong solution means that OTD can probabilistically always be fulfilled, weak solution means that it can be fulfilled under certain circumstances given by the Corollary of Post-Optimality just enounced before. We can therefore enounce the

Second Corollary to the Theorem of General Production Requirements (Corollary of Strong and Weak OTD Solutions)

Given is an order rate with generally random arrival character within a fixed-capacity regime. Necessary and sufficient condition for a strong OTD solution is the ergodic characteristic of a stochastic order process paired with asymmetric distribution of order rates of the random OR variable and if they comply with the requirements of the OTD theorem. However, non-ergodic processes with arbitrary arrivals, or ergodic processes with asymmetric random distributions of order rate, constitute only a weak solution for OTD observance that is depending on the specific circumstances. The circumstances are the distribution of arrivals, the sequence of arrivals, and timespan of non-randomness, as well as dependency on the virtual elasticity of the production system.

Please note that this just enounced corollary constitutes a constructive proof to Equation (12), whereas the Corollary of Post-Optimality belongs to the mathematical category of existence problem with non-constructive proof. The Corollary of Strong and Weak OTD Solutions states clearly that a Poissonian-distributed order rate process in this case might only constitute a weak solution to respect OTD.

An increased order entry is beneficial for every enterprise. With the utilization of the virtual elasticity and despite limited capacity, in certain cases, it is possible to observe the mandatory OTD in order to maintain a high related customer satisfaction referred to punctuality of deliveries. However, in extreme cases it might be necessary to increase capacity by additional shifts. In this case, additional aspects have to be taken into consideration to benefit from increased order intake:

Two aspects have to be satisfied

· Maximizing marginal manufacturing elasticity

· Maximizing margin contribution m(V)

These two conditions translate into (see

d d V k v ( V ) = 0

and

d d V m ( V ) = 0

The concept of optimal cost-level and optimal production-level are schematically shown in

However, the virtual elasticity allows to process increased order entry while still satisfying the OTD customer requirements. Interesting is that the variable production cost are not those of higher production volume, because the temporarily increased volume are produced at the same cost of the actual (not speed-up produced) volume V_{act} within the ΔT timespan, if staying in the same shift regime. Given that the maximum production level depends also on the price curve, the maximum volume V_{max} can go theoretically beyond the optimal volume V_{opt} to maximize margin contribution if the ΔT timespan is large enough (right graph of

From an economic point of view, the additional order entry is limited to the condition of marginal sales equaling marginal virtual elasticity variable cost, condition given in Equation (14)

V max : = { V | d d V p ( Δ V OTD ) = d d V k v ( V e l a s t ) } > V o p t (14)

Although Equation (14) might be of limited applicability in high performance production systems given by reduced ΔT it is for sure of theoretic importance. Nevertheless, we can come back to the Theorem of Optimal Production Range and therefore enounce the following lemma:

Lemma to the Theorem of Optimal Production Range (Lemma of Hidden Potential Exploitation)

Based on the realistic BEP model determining an optimal production volume to maximize margin contribution with an existent equipment, exploiting the property of virtual elasticity it is possible to extend maximization of margin contribution theoretically even beyond the common optimal production volume. This is possible due to the property of virtual elasticity of production systems and the resulting virtually increased capacity by speeding-up processes (Lemma of Increased Speed) in order to observe OTD also for increased order entry.

The finding summarized in Equation (14) and paraphrased with the Lemma of Hidden Potential Exploitation is an important revelation in the world of economic manufacturing and adds new insights of manufacturing-related properties of every production systems. This revelation is as well as important for practical reasons as for theoretical understanding of the value of speed. It helps to face the challenge of increased order entry with variable arrival rate above the nominal capacity allowing to satisfy OTD without physical capacity extension, at least for a limited timespan. Although this finding is very interesting, it has to be admitted that the virtual elasticity cannot replace consistent physical nominal extension of capacity and are limited within the scope of ΔT to observe OTD. It is ideally suited to absorb order peaks, but not structurally increased order levels.

It may emerge the question of the practical applicability, or better practical implementation, of these findings. Two related topics have to be investigated: firstly, finding manufacturing solutions to speed-up processes to increase virtual capacity, and secondly, finding scheduling algorithms allowing most of the backlogged orders to satisfy OTD. The first topic can be easily resolved by TPS-derived Lean JIT techniques, hence the importance to transform B&Q Systems into Lean JIT production systems [

· Shorten CT at the bottleneck (increase capacity without investing)

· Eliminate non value-add steps in the process such as searching (eliminate Muda)

· Reduce batch-size (in order to reduce WIP).

To generate significantly the virtual elasticity, however, introducing a SPF is mandatory, see the Main Theorem of Production Time (or SPF Dominance Theorem) and its corollary and lemma [

Lean JIT-based production systems show a superior performance versus traditional B&Q production systems and stand at the base to reduce PLT and respecting OTD [

Originally, Lean JIT production systems have been conceived for stable TR-based production regimes for the manufacturing of a deterministic product mix with adequate, but inelastic capacities. These production systems are largely representative for the automotive industry. However, in many industries, order entry is a random variable. This makes it difficult, but not impossible to implement JIT-like Lean systems. In addition, economy shows seasonal patterns requiring a variable capacity, not mentioning the variability originated by pipeline-filling effects of fragmented global supply chains during economic cycles. This complex variability of order entry may lead to costly investment in order to observe OTD resulting in production overcapacities. This does not have to be the case, at least not for limited variability and limited timespans, if the virtual elasticity of a production system is discovered, implemented, and exploited. The virtual elasticity is the result of the inherent latent production capacity resulting from the difference between EDT and PLT. This insight is an important result. Even though the findings of the paper may seem trivial, it has to be stated that this is not self-evident. Indeed, instead of investing in costly production capacities covering the upper range of order variability that is leading to overcapacities, the latent capacity intrinsic to every production systems given by increased speed of SPF and JIT production systems can be exploited especially for arrival distributions constituting a weak solution for OTD. It goes without saying that reducing PLT is not intended simply to operate the equipment at higher speeds, but by applying appropriate manufacturing and transfer principles in order to implement a high performance production system.

The authors declare no conflicts of interest regarding the publication of this paper.

Rüttimann, B.G. and Stöckli, M.T. (2020) Exploiting Virtual Elasticity of Production Systems for Respecting OTD-Part 1: Post-Optimality Conditions for Ergodic Order Arrivals in Fixed Capacity Regimes. American Journal of Ope- rations Research, 10, 321-342. https://doi.org/10.4236/ajor.2020.106018