# Patent

## Title
Method and business process for the estimation of erosion costs in assemble-to-order manufacturing operations

## Abstract

A method for forecasting the component surpluses for a target planning period is provided. To begin this method, a planner first identifies each component required to produce a product. For each component, the planner defines a planned level and an uncancelable level. The planned level for a component is the quantity at which the component is expected to be available. The uncancelable level for a component is the quantity of the component that cannot be liquidated without charge. The planner also defines a vector of connect rates for the components. After the required data has been entered, an expected surplus is computed for each component. To compute a component's expected surplus, the component is assumed to be available at its uncancelable level. The remainder of the components are assumed to be available at their respective planned levels. The mean production for the component is computed and used, along with the uncancelable level the selected component and the vector of connect rates to compute the component's expected surplus. This computation is repeated for all components in turn.

## Background
RELATED APPLICATIONS
This application claims priority to U.S. Provisional Application No. 60/229,611, entitled “Method And Business Process For Estimation Of Erosion Costs In Assemble-To-Order Manufacturing Operations.” filed Aug. 31, 2000, and naming Balazs Kralik, Michael Goldbach and Paul Dagum as inventors, which is related to U.S. Provisional Application No. 60/213,189, entitled “Method and Apparatus For The Estimation Of Mean Production For Assemble-To-Order Manufacturing Operations,” filed Jun. 21, 2000, and naming Paul Dagum, Michael Goldbach and Balazs Kralik as inventors and related to U.S. Pat. No. 6,684,193, entitled “Method and Apparatus for Optimizing Multivariate Allocation of Resources,” issued on Jan. 27, 2004, and naming Thomas Chavez and Paul Dagum as inventors and U.S. patent application Ser. No. 09/491,461, entitled “Method and Apparatus for Optimizing Multivariate Allocation of Resources,” filed on Jan. 26, 2000, and naming Thomas Chavez and Paul Dagum as inventors. This application is related to U.S. Pat. No. 7,249,049, entitled “Method and Business Process for the Estimation of Mean Production for Assemble-To-Order Manufacturing Operations,” issued on Jul. 24, 2007, and naming Paul Dagum, Michale Goldbach and Balazs Kralik as inventors. These applications are incorporated by reference herein, in their entirety and for all purposes.


TECHNICAL FIELD OF THE INVENTION
The present invention is generally related to resource-planning methods used by manufacturing companies and other organizations. More specifically, it relates to risk management in the procurement of unfinished goods by estimation of expected excess component inventory and erosion costs.
BACKGROUND OF THE INVENTION
In the field of product manufacturing, the term assemble-to-order (ATO) refers to systems where products are not produced until the demand for the products becomes known. Manufacture-to-stock (MTS) systems, on the other hand, base production on component availability. Unlike ATO manufacturing, MTS systems often produce products before the demand for those products has been determined.
Manufacturing companies typically prefer to use ATO manufacturing operations whenever possible. ATO allows manufacturers to minimize their component inventories. This reduces the amount of capital invested in components and reduces the risk that components will loose value or become obsolete before they are transformed and sold as products.
In an ATO operation, the availability of components is the key variable that controls the ability to produce. If component availability is limited (either because inventory is limited or because supplier commitments to ship the component on demand are insufficient) then fulfillment of demand either fails or is delayed. On the other hand, if excess inventory is left over, it must be liquidated or utilized later. Both outcomes are costly to the company. If the company liquidates, it may be forced to do so at substantially reduced prices. Many items lose value because of innovation. This is especially true in the especially fast moving high-tech sector where goods, such as microprocessor are continually superseded with cheaper more powerful replacements. In other industries there are shelf life considerations that make it costly to try to utilize excess inventory at a later time. In both cases there is capital tied up in extra inventory and that corresponds to incurring a cost equivalent to the firms funding rate.
The cost incurred for component inventory that exceeds the inventory required to meet demand is referred to as the erosion cost. The estimation of erosion costs is paramount in making intelligent business decisions about the proper amount of components to obtain commitments in any given period. This is particularly true when a large number of products are made from a large number of the same components.
SUMMARY OF THE INVENTION
An embodiment of the present invention includes a method for the computation of surplus components. For this method a planner first identifies each component required to produce a product. For each component, the planner defines a planned level and an uncancelable level. The planned level for a component is the quantity at which the component is expected to be available. The uncancelable level for a component is the quantity of the component that cannot be liquidated without charge. The planner also defines a vector of connect rates for the components.
After the required data has been entered, an expected surplus is computed for each component. To compute a component's expected surplus, the component is assumed to be available at its uncancelable level. The remainder of the components are assumed to be available at their respective planned levels. The mean production for the component is computed and used, along with the uncancelable level the selected component and the vector of connect rates to compute the component's expected surplus. This computation is repeated for all components in turn.
The expected component surpluses can then be used to compute the expected erosion cost per component and the total erosion cost for all components.
An assemble-to-order (ATO) manufacturing process is one in which products are manufactured from raw components only as orders are received. While ATO processes are more efficient than manufacture-to-stock (MTS) operations, in which products are produced before demand is known, they are not without significant risks. If component availability is below what is required to meet product demand, then fulfillment of orders is either delayed or fails. Conversely, if excess component inventory remains after demand is met, the excess is either returned to the supplier, usually at a loss, or held until future need. Thus any component inventory that does not exactly meet demand incurs some financial loss for the manufacturer. Being able to predict the excess component inventory and associated financial loss allows a manufacturer to make more informed business decisions.
The present invention addresses this need by providing a novel method for estimating component surplus and expected erosion costs. The erosion cost is simply the cost incurred for component inventory that exceeds the inventory required to meet demand. Estimation of erosion cost is particularly complicated when a large number of products are made from a large number of the same components, and thus the present invention is particularly advantageous under such conditions.
In addition to a business process and computer-implemented method, the present invention provides a system for computing expected component surplus and erosion costs. The system is preferably a computer system containing a memory, input/output devices, and processing means for implementing the method.
The present invention also provides a computer-readable medium of instructions executable by the computer to perform method steps for the method described above. The present invention applies equally regardless of the particular type of signal bearing media used. Examples of computer-readable medium include recordable-type media such as floppy disks and CD-ROMs and transmission-type media such as digital and analog communications links. It will be clear to one skilled in the art that the above embodiment may be altered in many ways without departing form the scope of the invention.
The foregoing is a summary and thus contains, by necessity, simplifications, generalizations and omissions of detail; consequently, those skilled in the art will appreciate that the summary is illustrative only and is not intended to be in any way limiting. As will also be apparent to one of skill in the art, the operations disclosed herein may be implemented in a number of ways, and such changes and modifications may be made without departing from this invention and its broader aspects. Other aspects, inventive features, and advantages of the present invention, as defined solely by the claims, will become apparent in the non-limiting detailed description set forth below.



BRIEF DESCRIPTION OF THE DRAWINGS
For a more complete understanding of the present invention and for further features and advantages, reference is now made to the following description taken in conjunction with the accompanying drawings, in which:
FIG. 1 is a block diagram of a computer system shown as a representative environment for deployment of the present invention.
FIG. 2 is a flowchart showing the steps associated with an embodiment of the method for computing expected surplus of the present invention.
FIG. 3 is a flowchart showing the steps associated with an embodiment of the method for computing expected erosion costs per component of the present invention.
FIG. 4 is a flowchart showing the steps associated with an embodiment of the method for computing total expected erosion cost of the present invention.



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
The preferred embodiments of the present invention and their advantages are best understood by referring to FIGS. 1 through 4 of the drawings. Like numerals are used for like and corresponding parts of the various drawings.
DEFINITIONS
Component plan: a list of quantities for each component, representing a company's component order for a given planning period.
Scenario: a set of assumptions about products and components. A scenario includes product parameters, component parameters, component consumption, component interactions, and an allocation policy.
OVERVIEW OF THE INVENTION
Problem Definition
When possible, it is beneficial for manufacturing companies to produce their outputs in an assemble-to-order (ATO) fashion. While manufacture-to-stock (MTS) manufacturers produce outputs before demand for them is revealed, ATO operations only produce outputs as the order stream arrives.
In an ATO operation, the availability of components is the key variable that controls the ability to produce. If component availability is limited (either because inventory is limited or because supplier commitment to ship the component on demand are insufficient), then fulfillment of demand either fails or is delayed.
On the other hand if excess inventory is left over, it must be liquidated or utilized later. Both outcomes are costly to the company. If the company liquidates, it can only do so at substantially reduced prices. Many items especially in the fast-moving high-tech sector lose value because of innovation (e.g. who wants last year's microprocessor when this year's is twice as fast for the same price). In other industries there are simple shelf life considerations that make it costly to try to utilize excess inventory at a later time. In addition, in both cases there is capital tied up in the extra inventory and that corresponds to incurring a cost equivalent to the firms funding rate.
The estimation of the erosion costs is paramount in making intelligent business decisions about the proper amount of components to obtain commitments for in any given period.
EXAMPLES PROBLEMS ADDRESSED BY THE PRESENT INVENTION


    
    
        1. How to compute the expected component surplus at the end of a particular planning period
        2. How to compute the value of expected erosion costs for a particular set of components in a planning period
        3. Computation of total expected erosion costs in a company for a particular planning period
    
    


The essence of this invention is the ability to solve the above problems in cases where there is a plurality of products and components.
ENVIRONMENT
In FIG. 1, a computer system 100 is shown as a representative environment for an embodiment of the present invention. Structurally, computer system 100 includes a processor, or processors 102, and a memory 104. Processor 102 can be selected from a wide range of commercially available or custom types. An input device 106 and an output device 108 are connected to processor 102 and memory 104, input device 106 and output device 108 represent all types of I/O devices such as disk drives, keyboards, modems, network adapters, printers and displays. Computer system 100 may also includes a disk drive 110 of any suitable disk drive type (equivalently, disk drive 110 may be any non-volatile mass storage system such as “flash” memory).
COMPUTING EXPECTED SURPLUS
As shown in FIG. 2, an embodiment of the present invention includes a Method 200 for the computation of surplus components (expected excess number of components remaining at the end of a period). Method 200 consists of the following steps:
In step 202 a user enters data describing one or more products and their associated components. The data entered includes all the data required for mean production computation as disclosed in the related application “Method and Business Process for the Estimation of Mean Production for Assemble-To-Order Manufacturing Operations.” In particular, this includes, for each product, information describing
1) The type of components required to produce the product,
2) The number of each component required (A),
3) Component allocations (d),
4) The mean demand for the product (μ).
5) The standard deviation of the demand for the product,
6) The correlations between the product and any other product.
This information can be input manually by the user or retrieved from a user database or other source.
For each component i, the value di (referred to above as component allocations) represents the maximum number of that is available (i.e. the expeditable amounts specified for the mean production computations). In many cases, positioned components may be returning or canceled. In these cases, only the uncancelable portion is subject to erosion. As part of step 202, the user enters (or otherwise inputs) a value di0 to represent the uncancelable portion of each component of interest. For each component i, the quantity di−di0 represents the number that can be liquidated or returned at no cost. For some components, return or liquidation of any amount has an associated cost. For each component i, of this type the value di0 is equal to di.
Steps 204 through 214 form a loop. This loop sequences each component i for which the calculation of surplus components is desired.
In step 206, each component plan di is set to its original value. This means that each component is assumed to be available in its planned quantity.
In step 208 the component plan for component i (i.e., di) is set to di0. This means that component i is assumed to be available at the uncancelable level for that component.
In step 210 the component levels specified in steps 206 and 208 (i.e., di0 for component i and di for all other components) are used to compute an expected mean production qi0 for a specified planning period. In general, there are several methods, such as Monte Carlo simulation that may be used to perform this computation. Preferably, but not necessarily, the method disclosed in the related application “Method and Business Process for the Estimation of Mean Production for Assemble-To-Order Manufacturing Operations.”
In step 212 the formula

ei=di0−ai·qi0 

is used to compute the expected number of erodible components i. In this formula di0 is the uncancelable level for component i, ai is the vector of connect rates for component i, and qi0 is the expected mean production computed in step 210. The uncancelable level for a component is the quantity of the component that cannot be liquidated without charge.

In step 216 the expected excess is reported to the user or returned to the system for further use.
COMPUTING EXPECTED EROSION COSTS PER COMPONENT
A shown in FIG. 3, an embodiment of the present invention includes a method 300 for the computing the expected erosion costs per component (the expected financial cost of having leftover inventory for a subset of components). Method 300 consists of the following steps:
In step 302 a user inputs the erosion costs per surplus component.
In step 304 the expected surplus is evaluated using Method 200. This step may be done in parallel with step 302.
In step 306 the expected surplus is multiplied by the per-unit erosion costs to obtain the expected erosion cost for each component of interest.
In step 306 the expected erosion cost for each component of interest is reported to the user.
In step 308 we report the expected erosion cost for each component of interest.
COMPUTING TOTAL EXPECTED EROSION COST
As shown in FIG. 4, an embodiment of the present invention includes a Method 400 for computing the total expected erosion cost (i.e., estimating the total erosion cost a company should expect given a certain plan for component availability). Method 400 consists of the following steps:
In step 402 the component-by-component expected erosion costs are evaluated for all components of the planning portfolio. This step may be done in parallel.
In step 404 the numbers obtained in step 402 are summed.
In step 406 the single expected erosion cost number is reported back to the user, or stored for further processing by the system.
The method of computing erosion costs can be used as the basis for computing expected cancellation fees using the following equation:




  
    〈
    
      C
      i
      
        (
        t
        )
      
    
    〉
  
  =
  
    
      〈
      
        E
        i
        
          (
          t
          )
        
      
      〉
    
    -
    
      〈
      
        E
        i
        
          (
          
            t
            -
            1
          
          )
        
      
      〉
    
  




The method of computing erosion costs can be used as the basis for computing expected expediting fees using the following equation:




  
    〈
    
      X
      i
      
        (
        t
        )
      
    
    〉
  
  =
  
    
      (
      
        
          d
          i
          
            (
            t
            )
          
        
        -
        
          d
          i
          
            (
            
              t
              -
              1
            
            )
          
        
      
      )
    
    -
    
      (
      
        
          〈
          
            E
            i
            
              (
              t
              )
            
          
          〉
        
        -
        
          〈
          
            E
            i
            
              (
              
                t
                -
                1
              
              )
            
          
          〉
        
      
      )
    
  




As will be appreciated, embodiments of the present this invention relate generally to the field of manufacturing resource planning. More specifically, the present invention relates to risk management in the procurement of unfinished goods.
Although particular embodiments of the present invention have been shown and described, it will be obvious to those skilled in the art that changes and modifications may be made without departing from the present invention in its broader aspects, and therefore, the appended claims are to encompass within their scope all such changes and modifications that fall within the true scope of the present invention.
1 REVIEW OF EROSION CUSTOM ASSEMBLE-TO-ORDER MANUFACTURING
When component stock exceeds component consumption a surplus remains. When the producer does not have the option to return unused components, the surplus is subject to value erosion, in at least three ways:

    
    
        opportunity costs of capital tied up in hard assets;
        technical erosion due to design changes and technological progress; and
        physical erosion due to spoilage.

As a result, the value of a surplus component is strictly less than its original value when it was purchased.

    
    


In some cases, the supply contract includes an option to return all or some of the component surplus subject to a cancellation fee. In the following analysis, we will assume that the cancellation fee is strictly less than the associated value erosion so that it will always be rational to exercise the cancellation option when it is available.
Similarly, a supply contract may also include an option to expedite extra components when components demand exceeds the stock. Again expediting will carry with it an expediting fee which can be thought of as a premium for rush service. Just as there is a lower limit to how many components the supplier is willing to take back, there is also an upper limit beyond which components cannot be expedited at any cost. When demand exceeds component availability the company may charge itself a “penalty” for missed (or delayed) sales.
A model is introduced of erosion and expediting based on a simplifying assumption that a distribution on component demand could be induced from product demand without regard to production policies. In this report we relax that assumption, recognizing that component consumption is not driven by product demand directly, but rather by production.
2. SIMPLE EROSION
In this section we consider the simple case of erosion when the supply contract does not allow for cancellation or expediting. As usual, we define product demand xεn as a random variable distributed normally with mean μ and covariance Σ. Given production qεn and connect rate matrix A εmxn induced component consumption is Aq. We denote the available component stock by dεm. Thus, the volume of component i subject to erosion is

ei=max(di−[Aq]i, 0)≡(di−[Aq]i)+.  (2.1)

Component supply d defines a feasible region Ω, in which component supply is sufficient to meet demand. The feasible region Ω is a defined as a polyhedral set

Ω=∩tεMHe,  (2.2)

where  denotes the set of relevant components, and He denotes the halfspace where component supply d is sufficient to meet demand, i.e.,

He={x:dt−[Ax], ≧0}.  (2.3)

In general, the feasible region can be product specific. In other words, the feasible region from the perspective of product j may be different that that of product k. We denote the feasible region with respect to product j as Ω(j), and the set of relevant components as (j) In particular, under uniform allocation (j) is the Markov Blanket associated with product j, where




  
    
      
        
          M
          ⁡
          
            (
            j
            )
          
        
        =
        
          ⋃
          
            i
            =
            D
          
        
      
    
    
      
        (
        2.4
        )
      
    
  




Under prioritized allocation feasibility The production function q:xε→qεΩ establishes a level of feasible production given product demand x and component availability d. In principle, it is desirable to produce as much as is feasible, but never more than what is demanded.

Therefore, when demand is feasible production maps demand to itself. When demand is not feasible, production maps infeasible demand to the feasible frontier

∂Ω=Ω∩(Υ,∂ll,)  (2.5)

where

∂Hi={x:di−[Ax]i=0}.

For clarity we will also denote the i-th row of A by ai. Thus production is defined as follows:





  
    
      
        
          q
          ⁡
          
            (
            
              x
              ,
              d
            
            )
          
        
        =
        
          {
          
            
              
                
                  x
                  ,
                
              
              
                
                  
                    
                      if
                      ⁢
                      
                        
                      
                      ⁢
                      x
                    
                    ∈
                    Ω
                  
                  ,
                
              
            
            
              
                
                  q
                  ∈
                  
                    ∂
                    Ω
                  
                
              
              
                
                  otherwise
                  .
                
              
            
          
        
      
    
    
      
        (
        2.6
        )
      
    
  




Expected Erosion. 
Expected Erosion is defined as




  
    
      
        
          〈
          
            E
            i
          
          〉
        
        =
        
          ∫
          
            
              
                (
                
                  
                    d
                    i
                  
                  -
                  
                    
                      [
                      Aq
                      ]
                    
                    i
                  
                
                )
              
              +
            
            ⁢
            f
            ⁢
            
              q
              )d 
              q
            
          
        
      
    
    
      
        (
        2.7
        )
      
    
  




Unfortunately, we do not have an explicit representation of the probability density function ∫(q). However, it is useful to recognize that by the definition of production in (2.6), qεΩ, and therefore

di−[Aq(x,d)], ≧0

for all iε[m]. Thus, expected erosion for component i can be computed as





  
    
      
        
          
            
              
                〈
                
                  E
                  i
                
                〉
              
              =
              
                ∫
                
                  
                    
                      (
                      
                        
                          d
                          i
                        
                        -
                        
                          
                            [
                            Aq
                            ]
                          
                          i
                        
                      
                      )
                    
                    +
                  
                  ⁢
                  
                    f
                    ⁡
                    
                      (
                      q
                      )
                    
                  
                  ⁢
                  
                    d
                    q
                  
                
              
            
          
        
        
          
            
              =
              
                ∫
                
                  
                    (
                    
                      
                        d
                        i
                      
                      -
                      
                        
                          a
                          i
                        
                        ·
                        
                          q
                          ⁡
                          
                            (
                            
                              x
                              ,
                              d
                            
                            )
                          
                        
                      
                    
                    )
                  
                  ⁢
                  
                    f
                    ⁡
                    
                      (
                      
                        
                          x
                          ;
                          μ
                        
                        ,
                        ∑
                      
                      )
                    
                  
                  ⁢
                  
                    d
                    x
                  
                
              
            
          
        
        
          
            
              =
              
                
                  d
                  i
                
                -
                
                  
                    a
                    i
                  
                  ·
                  
                    
                      〈
                      q
                      〉
                    
                    .
                  
                
              
            
          
        
      
    
    
      
        (
        2.8
        )
      
    
  




In this simple case, computation of expected erosion, reduces to computations of expected production q. Kralik and Dagum [2] present a method to compute expected production under the distributed production function

q(x,d)=x−au, where   (2.9)





  
    
      
        
          α
          ⁡
          
            (
            
              x
              ,
              d
            
            )
          
        
        =
        
          
            max
            i
          
          ⁢
          
            
              (
              
                
                  
                    
                      d
                      i
                    
                    -
                    
                      
                        a
                        i
                      
                      ·
                      x
                    
                  
                  
                    
                      a
                      i
                    
                    ·
                    u
                  
                
                ,
                0
              
              )
            
            .
          
        
      
    
    
      
        (
        2.10
        )
      
    
  




When u=x, this function satisfies the gate mapping property. Other values for u, i.e. u=μ, may not satisfy the gate mapping property over all xεn, but may represent reasonable approximations.

3 EROSION NET OF CANCELLATION
In some cases, the producer has the option to return or cancel positioned components. In this case, only the uncancelable portion of component surplus will be subject to erosion. We denote the uncancelable portion of available components by d(0). We continue to denote the full component availability as d and production q≡q(x, d), but we introduce the following additional notation

di(0)≡[di d2 . . . di(0) . . . dm]

Hi(0)≡{x:di(0)−[Ax]i≧0}

Ωi(0)≡Ω∩Hi(0)

qi(0)≡q(x,di(0))

In view of the above notation, a more precise definition of erosion is given by

Ei(0)=(di(0)−[Aq]i)+  (3.1)

Taking the expectation of Ei(0) in (3.1) presents a complication because production q is defined by the maximum component availability levels d, while erosion is computed with respect to the uncancelable levels di(0). Nevertheless, we will show that a simplification similar to (2.8) exists for a class of production functions which we refer to as gate mapping.
Gate Mapping Production
Given demand x, we define the gating component iε[m], to be the component for which the plan di represents a binding constraint on production, i.e. production q falls on the di-face of the feasible region. In order to proceed, consider the function




  
    
      
        
          
            g
            i
          
          ⁡
          
            (
            
              x
              ,
              d
            
            )
          
        
        =
        
          
            a
            i
          
          ·
          
            q
            ⁡
            
              (
              
                x
                ,
                d
              
              )
            
          
        
      
    
    
      
        (
        3.2
        )
      
    
  
  
    
      
        
          
            
              
                
              
              ⁢
              
                =
                
                  {
                  
                    
                      
                        
                          
                            a
                            i
                          
                          ·
                          x
                        
                      
                      
                        
                          
                            if
                            ⁢
                            
                              
                            
                            ⁢
                            x
                          
                          ∈
                          Ω
                        
                      
                    
                    
                      
                        
                          d
                          i
                        
                      
                      
                        
                          
                            if
                            ⁢
                            
                              
                            
                            ⁢
                            x
                          
                          ∈
                          
                            
                              Ω
                              _
                            
                            ⁢
                            
                              
                            
                            ⁢
                            and
                            ⁢
                            
                              
                            
                            ⁢
                            i
                            ⁢
                            
                              
                            
                            ⁢
                            is
                            ⁢
                            
                              
                            
                            ⁢
                            gating
                          
                        
                      
                    
                    
                      
                        
                          
                            
                              
                                g
                                _
                              
                              i
                            
                            ⁡
                            
                              (
                              
                                x
                                ,
                                d
                              
                              )
                            
                          
                          ≤
                          
                            d
                            i
                          
                        
                      
                      
                        
                          otherwise
                          .
                        
                      
                    
                  
                
              
            
          
          
            
              
            
          
        
      
    
    
      
        (
        3.3
        )
      
    
  




In the above definition, gi is some function which depends on the problem details. Its region of definition Ψi is the intersection of Ω with the complement of the region where i is gating.

For the usual local·μ production policy [2], gi has an additional property. In its region of definition, it is independent of di. This is the gate mapping property:
Definition 3.1 
Gate mapping property: an allocation policy is gate mapping if the following condition holds for every x: If for both d=(d1, . . . di, . . . dn) and d′=(d1, . . . , di, . . . , dn) the gating component is the same j≠i, then q(x, d)=q(x, d)=q(x, d′).
Lemma 3.2
When production q is gate mapping as defined above, the production mappings q and qi(0) are equivalent with respect to demand for which there is erosion associated with component i. In math: q=qi(0) for all xε{x:di(0)−[Aq],>0}.
Proof
There are two cases in which the condition di(0)>[Aq], holds:

    
    
        1. xεΩi(0)⊂Ω, or
        2. q maps xεΩ to the feasible frontier ∂Ω such that qεHs(0).
    
    


In the first case, q=qi(0)=x because under the extra constrain xεHs(0) had no effect on feasibility. In the second case, there exists a component k≠i for which the feasibility constrain xεHk is more constraining than the new constraint xεHs(0). Since production is gate mapping, infeasible demand is mapped to the part of the feasible frontier ∂Hk associated with the gating component k. Since the new constraint is not binding, it does not change production mapping.
Lemma 3.3
If there is no erosion associated with demand x when it is mapped by a gate mapping function q, there is also no erosion in component i when production is mapped by qi(0).
Proof
Let φ={x:di(0)−a,·q(x,d)), ≦0}. Because production is gate mapping, q maps xεφ to a point q εHi(0). Therefore, when the feasible region is constrained to Ωi(0) by adding the constraint qεHi(0), component i becomes gating for all xεφ and therefore qi(0) will map the demand to the face of ∂Ωi(0) defined by the hyperplane ∂Hi(0)=}q:di(0)−[Aq], =0} along which there is no erosion.
Lemma 3.4
For all xεn.di(0)−[Aq], =di(0)−[Aqi(0)]i.
Proof
The proof follows directly from Lemmas 3.2 and 3.3. From Lemma 3.2 we see that when there is erosion, the production mapping q=qi(0). From Lemma 3.3 we see that when there is no erosion under q, the new mapping qi(0) maps demand to the hyperplane ∂Hi(0) where there is no erosion
We now state the following two propositions:
Proposition 3.1
When production is gate mapping property, expected erosion for component i can be comprised as




  
    
      
        
          〈
          
            E
            i
            
              (
              0
              )
            
          
          〉
        
        =
        
          
            d
            i
            
              (
              0
              )
            
          
          -
          
            
              a
              i
            
            ·
            
              
                〈
                
                  q
                  i
                  
                    (
                    0
                    )
                  
                
                〉
              
              .
            
          
        
      
    
    
      
        (
        3.4
        )
      
    
  




Proof
This proposition follows directly from Lemma 3.4 and Equation (2.8).
Proposition 3.2
The production function defined in (2.9) and (2.10) is gate mapping.
Proof
The max function in (2.10) ensures that demand is mapped back to hyperplane associated with the gating component.
Therefore we can compute (3.4) using available methods.
4 CANCELLATION
In the previous section we demonstrated how to compute the expected volume of erodible components. In this section we consider the expected volume of cancelled components. Cancelled components are the surplus components that can be returned to the supplier and therefore are not subject to erosion.
4.1 Simple Cancellation
As before, we denote the uncancelable volume of components by d(0). We now denote the component plane by d(1), where the component plan represents the amount of component positioned. Cancellation only applies to components that have been positioned. We continue to denote available components (which may include expediting) by d.
Given a production q such that qεHi(1), the number of components i that can be cancelled is

Ci=di(1)−max(di(0), [Aq]i).  (4.1)

The cancellation is otherwise 0.
Proposition 4.1
When production is gate mapping, expected cancellation can be computed from two erosion computations:




  
    
      
        
          
            〈
            
              C
              i
            
            〉
          
          =
          
            
              〈
              
                E
                i
                
                  (
                  1
                  )
                
              
              〉
            
            -
            
              〈
              
                E
                i
                
                  (
                  0
                  )
                
              
              〉
            
          
        
        ,
      
    
    
      
        (
        4.2
        )
      
    
  
  
    
      
        where
      
    
    
      
        
      
    
  
  
    
      
        
          〈
          
            E
            i
            
              (
              t
              )
            
          
          〉
        
        ≡
        
          
            d
            i
            
              (
              t
              )
            
          
          -
          
            
              a
              i
            
            ·
            
              
                〈
                
                  q
                  i
                  
                    (
                    t
                    )
                  
                
                〉
              
              .
            
          
        
      
    
    
      
        
      
    
  




Proof
Expanding the definition in (4.1), yields




  
    C
    i
  
  =
  
    {
    
      
        
          
            
              
                d
                i
                
                  (
                  l
                  )
                
              
              -
              
                
                  [
                  Aq
                  ]
                
                i
              
              -
              
                (
                
                  
                    d
                    i
                    
                      (
                      0
                      )
                    
                  
                  -
                  
                    
                      [
                      Aq
                      ]
                    
                    i
                  
                
                )
              
            
            ,
          
        
        
          
            
              
                if
                ⁢
                
                  
                
                ⁢
                q
              
              ∈
              
                H
                i
                
                  (
                  0
                  )
                
              
            
            ,
          
        
      
      
        
          
            
              d
              i
              
                (
                l
                )
              
            
            -
            
              
                [
                Aq
                ]
              
              i
            
          
        
        
          
            
              
                if
                ⁢
                
                  
                
                ⁢
                
                  d
                  i
                  
                    (
                    0
                    )
                  
                
              
              <
              
                
                  [
                  Aq
                  ]
                
                i
              
              <
              
                d
                i
                
                  (
                  l
                  )
                
              
            
            ,
          
        
      
      
        
          0
        
        
          
            otherwise
            .
          
        
      
    
  




Taking the expected value of





  〈
  
    C
    i
  
  〉




over xεHn implies that





  
    
      
        
          〈
          
            C
            i
          
          〉
        
        =
        
          
            〈
            
              
                (
                
                  
                    d
                    i
                    
                      (
                      1
                      )
                    
                  
                  -
                  
                    
                      [
                      Aq
                      ]
                    
                    i
                  
                
                )
              
              +
            
            〉
          
          -
          
            〈
            
              
                (
                
                  
                    d
                    i
                    
                      (
                      0
                      )
                    
                  
                  -
                  
                    
                      [
                      Aq
                      ]
                    
                    i
                  
                
                )
              
              +
            
            〉
          
        
      
    
    
      
        (
        4.3
        )
      
    
  




We recognize the first term in (4.3) as simple erosion equivalent to





  
    〈
    
      E
      i
      
        (
        1
        )
      
    
    〉
  
  .




We recognize that the second term in (4.3) is equivalent to





  〈
  
    E
    i
    
      (
      0
      )
    
  
  〉




following the discussion of erosion with cancellation in §3.

4.2 MULTILEVEL CANCELLATION
In principle, there may be many levels of cancellation (or expediting). We denote a cancellation or expediting level by l, where 0≦l≦L, and the associated volumes of components by d(l), where di(0)<. . . <di(l)< . . . di(L). Thus d=d(L) represents components available for production. We can now expands the analysis to multiple levels of cancellation and introduce the following expanded notation:

d(l)≡{di(l) . . . dm(l)},

H(l)≡{x:di(l)[Ax], ≧0},

Ω(l)≡∩iε[im]He(l), and

q(l)≡q(x, d(l)).

As before

d≡d(L),

Hi≡H(L),

Ω≡Ω(L), and

q≡q(L).

Following the reasoning in §4.1, expected cancellation of level 0<l<L components is simply




  
    
      
        
          〈
          
            C
            i
            
              (
              t
              )
            
          
          〉
        
        =
        
          
            〈
            
              E
              i
              
                (
                t
                )
              
            
            〉
          
          -
          
            
              〈
              
                E
                i
                
                  (
                  
                    i
                    -
                    1
                  
                  )
                
              
              〉
            
            .
          
        
      
    
    
      
        (
        4.4
        )
      
    
  




EXPEDITING
Expediting is essentially the mirror problem to cancellation. However, instead of computing component surplus, we are computing component use. Thus while cancellation measured the nonnegative distance to an upper boundary, expediting measures then nonnegative distance to a lower boundary.
Let d(l) represent the components available at expediting level l. We now define expediting as associated with component i as




  
    
      
        
          X
          i
          
            (
            l
            )
          
        
        =
        
          {
          
            
              
                
                  
                    
                      min
                      ⁡
                      
                        (
                        
                          
                            
                              [
                              Aq
                              ]
                            
                            i
                          
                          ⁢
                          
                            d
                            i
                            
                              (
                              l
                              )
                            
                          
                        
                        )
                      
                    
                    -
                    
                      d
                      i
                      
                        (
                        
                          l
                          -
                          1
                        
                        )
                      
                    
                  
                  ,
                
              
              
                
                  
                    if
                    ⁢
                    
                      
                    
                    ⁢
                    q
                  
                  ∈
                  
                    
                      H
                      _
                    
                    i
                    
                      (
                      
                        l
                        -
                        1
                      
                      )
                    
                  
                
              
            
            
              
                
                  0
                  ,
                
              
              
                
                  otherwise
                  .
                
              
            
          
        
      
    
    
      
        (
        5.1
        )
      
    
  




Proposition 5.1
When production is gate mapping, expected expediting can be computed as




  
    
      
        
          〈
          
            X
            i
            
              (
              t
              )
            
          
          〉
        
        =
        
          
            (
            
              
                d
                i
                
                  (
                  t
                  )
                
              
              -
              
                d
                i
                
                  (
                  
                    t
                    -
                    1
                  
                  )
                
              
            
            )
          
          -
          
            
              (
              
                
                  〈
                  
                    E
                    i
                    
                      (
                      t
                      )
                    
                  
                  〉
                
                -
                
                  〈
                  
                    E
                    i
                    
                      (
                      
                        t
                        -
                        1
                      
                      )
                    
                  
                  〉
                
              
              )
            
            .
          
        
      
    
    
      
        (
        5.2
        )
      
    
  




Proof
We follow the same logic as before. Expanding the definition in (5.1), yields




  
    X
    i
    
      (
      l
      )
    
  
  =
  
    {
    
      
        
          
            
              
                (
                
                  
                    
                      [
                      Aq
                      ]
                    
                    i
                  
                  -
                  
                    d
                    i
                    
                      (
                      
                        l
                        -
                        1
                      
                      )
                    
                  
                
                )
              
              -
              
                
                  [
                  Aq
                  ]
                
                i
              
              -
              
                d
                i
                
                  (
                  l
                  )
                
              
            
            ,
          
        
        
          
            
              
                if
                ⁢
                
                  
                
                ⁢
                q
              
              ∈
              
                
                  H
                  _
                
                i
                
                  (
                  l
                  )
                
              
            
            ,
          
        
      
      
        
          
            
              
                
                  [
                  Aq
                  ]
                
                i
              
              -
              
                d
                i
                
                  (
                  
                    l
                    -
                    1
                  
                  )
                
              
            
            ,
          
        
        
          
            
              
                if
                ⁢
                
                  
                
                ⁢
                
                  d
                  i
                  
                    (
                    
                      l
                      -
                      1
                    
                    )
                  
                
              
              <
              
                
                  [
                  Aq
                  ]
                
                i
              
              <
              
                d
                i
                
                  (
                  l
                  )
                
              
            
            ,
          
        
      
      
        
          
            0
            ,
          
        
        
          
            otherwise
            .
          
        
      
    
  




Taking the expected value of Xi(l) yields




  
    〈
    
      X
      i
      
        (
        t
        )
      
    
    〉
  
  =
  
    
      〈
      
        
          (
          
            
              
                [
                Aq
                ]
              
              i
            
            -
            
              d
              i
              
                (
                
                  t
                  -
                  1
                
                )
              
            
          
          )
        
        +
      
      〉
    
    -
    
      
        〈
        
          
            (
            
              
                
                  [
                  Aq
                  ]
                
                i
              
              -
              
                d
                i
                
                  (
                  t
                  )
                
              
            
            )
          
          +
        
        〉
      
      .
    
  




With a simple transformation we can recognize these as already known quantities:




  
    
      
        
          〈
          
            
              (
              
                
                  
                    [
                    Aq
                    ]
                  
                  i
                
                -
                
                  d
                  i
                  
                    (
                    l
                    )
                  
                
              
              )
            
            +
          
          〉
        
        =
        
          
            〈
            
              
                
                  [
                  Aq
                  ]
                
                i
              
              -
              
                d
                i
                
                  (
                  l
                  )
                
              
            
            〉
          
          +
          
            〈
            
              
                (
                
                  
                    d
                    i
                    
                      (
                      l
                      )
                    
                  
                  -
                  
                    
                      [
                      Aq
                      ]
                    
                    i
                  
                
                )
              
              +
            
            〉
          
        
      
    
  
  
    
      
        =
        
          
            
              a
              i
            
            ·
            
              〈
              q
              〉
            
          
          -
          
            d
            i
            
              (
              l
              )
            
          
          +
          
            〈
            
              E
              i
              
                (
                l
                )
              
            
            〉
          
        
      
    
  




With this, we can write




  
    
      
        
          〈
          
            X
            i
            
              (
              t
              )
            
          
          〉
        
        =
        
          
            (
            
              
                d
                i
                
                  (
                  t
                  )
                
              
              -
              
                d
                i
                
                  (
                  
                    t
                    -
                    1
                  
                  )
                
              
            
            )
          
          -
          
            
              (
              
                
                  〈
                  
                    E
                    i
                    
                      (
                      t
                      )
                    
                  
                  〉
                
                -
                
                  〈
                  
                    E
                    i
                    
                      (
                      
                        t
                        -
                        1
                      
                      )
                    
                  
                  〉
                
              
              )
            
            .
          
        
      
    
    
      
        (
        5.3
        )

## Claims

1. A computer implemented method:
estimating an expected component surplus, wherein said estimating comprises
identifying a plurality of components required to produce a product;
defining a vector of connect rates for the components;
defining an uncancelable level for each of the components, each of the uncancelable levels defining a quantity of the each of the components below which the quantity of the each of the components cannot be liquidated without incurring a charge;
assuming that a selected one of the components is available at least at an uncancelable level of the uncancelable levels corresponding to the selected one of the components; and
computing the expected component surplus for the selected component using a mean production for the product, the uncancelable level and the vector of connect rates, wherein
the computing is performed by a computer.

2. The computer implemented method of claim 1, the estimating further comprising:
defining a planned level for the each of the components, each of the planned levels defining a quantity at which a corresponding one of the components is expected to be available.

3. The computer implemented method of claim 2, the estimating further comprising:
assuming that each of the remaining ones of the components is available at the planned level corresponding to the each of the remaining ones of the components, wherein the remaining ones of the components are ones of the components other than the selected one of the components.

4. The computer implemented method of claim 3, the estimating further comprising:
using the planned levels in the computing the expected component surplus.

5. The computer implemented method of claim 4, the estimating further comprising:
estimating the mean production for the product.