#### Abstract

Owing to the transparency in supply chains, enhancing competitiveness of industries becomes a vital factor. Therefore, many developing countries look for a possible method to save costs. In this point of view, this study deals with the complicated liberalization policies in the global supply chain management system and proposes a mathematical model via the flow-control constraints, which are utilized to cope with the bonded warehouses for obtaining maximal profits. Numerical experiments illustrate that the proposed model can be effectively solved to obtain the optimal profits in the global supply chain environment.

#### 1. Introduction

Traditionally, supply chain management (SCM) mainly offers different ways to reduce the production and transportation costs such that either the total expenditures in a supply chain can be minimized [1, 2] or the profits can be maximized [3, 4] to enhance the industrial competitiveness. These concepts above have been formed as SCM mathematical models in the last few decades such as supplier’s pricing policy in a just-in-time environment [5], pricing strategy for deteriorating items using quantity discount when customer demand is sensitive [6], and an optimization approach for supply chain management models with quantity discount policy [7].

On the other hand, a part of studies is focused on demand forecasting [8–11] to decrease the impact of the bullwhip effect [12, 13]. Additionally, some studies utilize the neural network method and regression analysis to improve the accuracy for forecasting customer demand [14–16], and furthermore, these connection weights in neural networks are assigned weighting based on fuzzy analytic hierarchy process methods without any tunings [14, 17]. By the previous mentions, we cannot guarantee that the solution are optimal. Based on the reason, this study proposes a deterministic method, an approximate method for linearizing nonlinear time series analysis model, to forecast customer demand in Section 3.2.

Many researchers have suggested that information sharing is a key influence on SCM environments [11, 18] and that it impacts the SCM performance in terms of both total costs and service levels [11, 19, 20]. However, the development of a global SCM model must share information as much as possible. We then suppose that the information in our experimental tests is shared to simplify the SCM situations [21–23].

In order to enhance the competitiveness of industries, a liberalization policy discussed in recent years becomes an important factor. Many developing countries have implemented a liberalization policy with bonded warehouses in the global supply chain environments. The liberalization policy with bonded warehouses in the different countries around the world, makes the goods be stored in the bonded warehouse, the liability of exporter is temporally cancelled, and the importer and warehouse proprietor incur liability. A bonded warehouse is one building or another secured area, where dutiable goods may be stored or processed without payment of duty. It may be managed by the state or by a private enterprise. The liberalization policy with bonded warehouses not only enhances the industries competitiveness but also saves inventory, production and transportation costs. Of course, the objective is to find more bonded warehouses, and more costs can be saved. In this study, the policy has proved the advantage for the exporter obtaining more profit in our numerical examples [24–28].

A global SCM is seen as a network that includes upstream, midstream and downstream sectors linkage, where we consider vendors as upstream, facilities and warehouses as midstream, and retailers (or customers) as downstream. The liberalization policy focuses on bonded warehouses of midstream in different countries around the world is not considered in the current SCM model [7, 29]. In this study, we then propose a global SCM model to be suitable for the international liberalization policy. The objective is to maximize the profits associated with inventory, production, transpiration, and distribution in the period of times. The global SCM model is formed as a linear mixed-integer program and solved by commerce software (i.e., [30]). Then, a flow chart for the proposed procedure is depicted in Figure 1.

This paper is organized as follows. Section 2 defines the notations in global SCM and forms a basic model of global SCM. Section 3 describes a procedure of global SCM. Section 4 shows the system architecture of a global SCM. Finally, numerical examples demonstrate that the proposed global SCM method can enhance profit ratio in Section 5.

#### 2. A Global SCM Model

There are four stages in the global SCM environment discussed: the vendor, facility, warehouse, and retailer. Each stage may be located in different nations around the world. The general structure of a global supply chain network has been displayed in Figure 2. For simplifying presentation, both income and cost statements occur in the dashed rectangle as follows.(i)Income statement(a)sales of goods in warehouse,(b)transfer price of goods in facility. (ii)Cost statement(a)procurement cost from vendor,(b)transportation cost (i.e., vendor to facility, facility to warehouse, and warehouse, to retailer),(c)inventory cost (i.e., warehouse).

Notations are described in Section 2.1, and a basic model of global SCM is formulated in Section 2.2.

##### 2.1. Notations

There are seven entity sets : The set of nations : The set of vendors : The set of facilities : The set of warehouses : The set of retailers : The set of times : The set of parts (or goods) for presenting each entity associated with notations as shown above.

All parameters are composite style, which are composed of uppercase English strings and suffixes with one or more lowercase alphabet characters, and each suffix single lowercase alphabet character presents the element of relative entity sets. These parameters are : The sales price of part in warehouse . : The transfer price of part in facility . : The procurement cost of part from vendor . : The production cost of part in facility . : The basis of transportation costs from vendor to facility . : The basis of transportation costs from facility to facility . : The basis of transportation costs from facility to warehouse . : The basis of transportation costs from warehouse to retailer . : The basis of inventory costs in facility . : The basis of inventory costs in warehouse . : The transportation time from vendor to facility . : The transportation time from facility to facility . : The transportation time from facility to warehouse . : The transportation time from warehouse to retailer . : The production time of part in facility . : The amount for part reproducing as part . : The weighting of transportation costs for part. : The weighting of inventory costs for part. : The required quantities of part for retailer at time . : The tax in nation . : The import tax of part in nation . * *: The vendor in nation . : The facility in nation . : The warehouse in nation . : The retailer in nation . : The supplied upper level of part in the facility . : The supplied upper level of part from vendor to facility . : The supplied upper level of part from facility to facility . : The supplied upper level of part from facility to warehouse . : The supplied upper level of part from warehouse to retailer . : The inventoried upper level of part in facility . : The inventoried upper level of part in warehouse .

All decision variables in global SCM are also compositestyle with lowercase suffix alphabet characters. The details of decision variables (integer variables) are described below: : The amount of part transported from vendor to facility at time : The amount of part produced in facility at time : The amount of part transported from facility to facility at time : The amount of part transported from facility to warehouse at time : The amount of part transported from warehouse to retailer at time : The inventory of part in facility at time : The inventory of part in warehouse at time : The amount of bonded part transported from vendor to facility at time : The amount of bonded part produced in facility at time : The amount of bonded part transported from facility to facility at time : The amount of bonded part transported from facility to warehouse at time : The inventory of bonded part in facility at time .

##### 2.2. A Proposed Model

We separately use , , and as the symbols of sales profits, losses, and taxes in each nation at different times. The objective of the proposed model is to maximize the sales profits associated with inventory, production, transpiration, and distribution in each nation at a period of times. The linear global SCM model is formulated as follows.

*A Global SCM Model*
where all decision variables (i.e., , etc.) belong to integers.

The objective is to maximize the profits (i.e., expression (1)) in the global SCM model. The revenue of sales (or losses) are obtained from the following statements.

*Income Statement* (2): sales of part in warehouse minus transportation costs from warehouse to retailer , (7): transfer price of part from facility to warehouses (i.e., transfer price), (8): transfer price of part from facility to facility (i.e., transfer price),

*Cost Statement* (3): transportation and ordering costs from facility to warehouse , (4): import duty of part from facility to warehouse , (5): inventory costs of part in warehouse , (6): fixed cost in warehouse , (9): transportation and ordering costs from facility to warehouse , (10): import duty of part from facility to facility, (11): transportation and ordering costs from vendor to facility , (12): import duty of parts from vendor to facility, (13): production cost in facility, (14): inventory cost in warehouse , (15): fixed costs in facilities.

Continuously, there are two kinds of constraints in the proposed model. One is flow conservation and the other is upper-lower bound constraints, which are described below.

*Flow Conservation*

We have These flow conservation constraints are described as follows: (16): the limitation of the flow conservation in each facility, (17): the limitation of the flow conservation in each warehouse, (18): the limited bonded parts of the flow conservation in each facility, (19): the limited bonded parts of the flow conservation in each warehouse, (20): all of the supplied goods satisfy with retailer demand.

*Upper-Lower Bound*

We have for ,, and ,,, and .

The upper-lower bound constraints are noted as follows: (21): capacities in which facilities produce goods, (22): inventory of parts in facilities, (23): inventory of parts in warehouses, (24): quantities of raw parts transported from vendors to facilities, (25): quantities of parts transported from facilities to facilities, (26): quantities of parts transported from facilities to warehouses, (27): quantities of parts transported from the warehouse to retailers, (28): quantities of bonded parts transported from vendors to facilities, (29): quantities of bonded parts transported from facilities to other facilities, (30): quantities of bonded parts transported from facilities to warehouses.

#### 3. A Procedure of Global SCM

There are different procedures in a global supply chain environment such as optimal transportation in logistic management, demand forecasting in preprocessing, bill of materials’ (*BOM*) table, and considering the liberalization policies with bonded warehouses. These procedures are discussed below.

##### 3.1. Logistic Management

In this study, we assume that a global SCM goes through four stages such as procurement, production, warehousing, and distribution (i.e.,, where the semifinished goods will be transferred to another facility for reproduction (i.e., ) to enhance the transfer price of goods. The course of the logistic management needs to go through a phase or a cycle (i.e., ). The processing flow of global SCM is presented in Figure 3.

In order to avoid delaying the demand of retailers in the future, in the beginning, we let the quantities of inventories equal zero in each stage in the first phase. Then, we will add times for utilizing demand forecast method. The period of times will be calculated in this global. Occurring in the period from to times, all of the goods in transportation and processing will be restored to the inventories of the adjacent warehouses.

##### 3.2. Demand Forecast Method—Time Series Analysis

As mentioned above, we need to forecast the quantity of retailers’ (i.e., customer) demand. Time series forecasting is a well-known method [31, 32] to forecast future demands based on known past events.

By referring to Li’s method [33], we then have the following nonlinear time series analysis’ model.

*A Nonlinear Demand Forecasting Model (NDF Model)*
where is sum of seasons (usually ) and , , and denotes the sales, quarter numbers, and seasonal indexes, respectively. For example, if , then represent the first quarter, represent the second quarter, represent the third quarter, and represent the fourth quarter.

A NDF model is nonlinear because of existing absolute functions and nonlinear terms. Following linear systems can linearize them. The linear approximation of the function is a continuous function , which is equal to (i.e., at each grid point ) and a linear function of its variables on the interior or edges of each of the indicated triangles. To describe the linear approximation of the function , we then have the following linear system (referring to the models in Babayev [34]).

Let where and .

Expressions (36)~(38) define the point () as the convex linear combination of the vertices of the chosen triangle. Expressions (39)~(41) imply that if , then and *;* or if , then and .

The mentions above give a basic formulation satisfying the piecewise linear approximation of . Similarly, can be also linearized by the same way (for the proof, refer to Babayev [34]).

On the other hand, we use additional binary variables to treat the absolute function. A piecewisely linear demand forecasting model can be reformulated as follows.

*A Piecewisely Linear Demand Forecasting Model (PLDF Model)*
where is a big enough constant, , , , and are piecewisely linear variables of and , respectively, referred to(35)–(41).

Considering the following situations we know that the nonlinear objective function can be linearized by linear system, that is, (44) and (45).(i)if , then , which is forced by (44) because of ,(ii)if , then , which is forced by (45) because of .

##### 3.3. Bill of Materials (*BOM*)

A bill of materials’ (*BOM*) table is expressed as the “parts list” of components. In order to understand the *BOM* for this global SCM model, the symbols ofandcould be expressed as a material or goods (i.e., parameter ). Hence, we can utilize a two-dimensional matrix to express the relationship between *partial *and *goods*. For example, a with a can produce a , and a with two can produce a . The *BOM* table of the sparse matrix is expressed in Figure 4. The aim is to design a friendly interface for our system design.

##### 3.4. Liberalization Policy with Bonded Warehouses

This study proposes a global SCM model which is capable of treating liberalization policy for bonded warehouses in the different countries around the world [24–28]. A bonded warehouse is one building or any other secured area, where dutiable goods may be stored or processed without payment of duty. It may be managed by the state or by private enterprise.

While the goods are stored in the bonded warehouse, the liability of exporter is temporally cancelled, and the importer and warehouse proprietor incur liability. These goods are in the bonded warehouse, and they may be processed by assembling, sorting, repacking, or others. After processing and within the warehousing period, the goods may be exported without the payment of duty, or they may be withdrawn for consumption upon payment of duty at the rate applicable to the goods in their manipulated condition at the time of withdrawal. Bonded warehouses provide specialized storage services such as deep freeze, bulk liquid storage, commodity processing, and coordination with transportation and are an integral part of the global supply chain.

The liberalization policy with bonded warehouses not only enhances the industries competitiveness but also saves inventory, production, and transportation costs. Of course, the objective aims to find more bonded warehouses, and more costs could be saved.

#### 4. System Development

Currently, global supply chain management (global SCM) system always needs experts offering the qualities of goods in bonded warehouses to calculate the production capacity since there is a lack of kernel technique (i.e., deterministic mathematical model) for developing global SCM system. At the same time, considering the policy of trade liberalization with bonded warehouses by manual is a difficult job. Therefore, this study develops a global SCM system with optimization technique to cope with the bonded warehouses for obtaining the maximal profits, automatically.

##### 4.1. System Design

This study designs a system to apply for a global SCM system based on a proposed mathematical model. The architecture of the global SCM system includes an application server in the middle tier to distribute all of the processing jobs. Establishing a Broker Server can not only distribute the jobs into different application servers, but also the jobs work much more efficient. The basic architecture of the global SCM system is shown in Figure 5.

Furthermore, the system architecture includes three tiers which are composed of various modules. In the first tier, GUI layer, an intuitive and friendly interface for users is built. In the second tier, the application (AP) server layer, is a kernel in the global SCM system whose base is based on the Java programming language and an optimal software component (i.e., LINGO’s Solver). In the last tier, database layer stores all of the data such as the username, password, environment setting, and logistics information into database system. Moreover, the data must be backed up to the other database to protect it from hacking or crashing. Then, the system architecture is shown in Figure 6.

#### 5. Numerical Examples

In our numerical examples, we assume that the situation occurred in three different overseas locations. There are three vendors (), four facilities (), four warehouses (), three retailers (). Under the liberalization policies, these limited quantities of bonded warehouses are listed in the Appendix. There are three parts () supplied from the vendor to the facility, and the parts be produced another two parts (,) in different facilities. The period of time includes two demands: one is actual customer demand (i.e., ), and the other is forecasting demand (i.e., ), so that the total period of time is 17 (i.e., to 17). The other parameters are expressed in Tables 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, and 26 of the Appendix. The numerical example was solved via a global SCM system; these logistic results at different times are depicted in Figures 7, 8, 9, 10 and 11.

Solving the example, the objective value is $603,010 with consideration of the bonded warehouses. We then calculate them without bonded warehouses again, and the objective value is 514,680. Comparing both of them, the profit ratio increases to 14.6% as shown in Table 1.

In order to observe the advantages of the proposed method, this study experiments different numerical examples by generating parameters below.(i)The total locations are randomly generated from a uniform distribution according to 2 ≤ vendors ≤ 20, 3 ≤ facilities ≤ 20, 4 ≤ warehouses ≤ 20, and 2 ≤ retailers ≤ 40.(ii)The period of time is randomly generated from the uniform distribution over the range of [2, 16].(iii)The capacities of bonded warehouses are according to the capacities of original warehouses based on uniform distribution from 50% to 90%.(iv)In order to balance between the supply and demand, all of the coefficients follows feasible sets for the global SCM models.According to the computational results followed rules (i) to (iv), we observe that the liberalization policy with bonded warehouses directly affects the activities of the global SCMs. When we fixed the number of bonded warehouses to 5, the profit ratio approaches 15%. When we take the number of bonded warehouses being 40 as an example, the profit ratio enhances to 20%. Finally, the profit ratio is close to 25% at 50 bonded warehouses. That is, considering that more nodes have bounded constraints, there is a greater cost saving. The tendency of the other 10 tests is depicted in Figure 12.

The proposed model considers the liberalization policy with bonded warehouses as flow-control constraints to obtain optimal profits. Additionally, numerical experiments also illustrate that the proposed model can enhance the profit ratio when the number of bonded warehouses is large.

#### 6. Conclusions

This study proposes a global SCM model which is capable of treating liberalization policy with bonded warehouses. All of the constraints about bonded warehouses can be added to the proposed global SCM model. Additionally, the proposed model is embedded into a global SCM system we developed. The system can automatically calculate the profit ratio to further enhance the competitiveness. Numerical experiments’ reports illustrate that the proposed model can be effectively solved to obtain the optimal profit.

Future research can consider different quantity discount policies, artificial neural network method for solving the uncertainty demand forecast, heuristic methods for accelerating solving speed, and linearization technique for resolving the proposed demand forecasting model.

#### Appendix

For more details, see Tables 2–26.

#### Acknowledgments

The authors would like to thank the editor and anonymous referees for providing most valuable comments for us to improve the quality of this paper. This research was supported by the project granted by ROC NSC 101-2811-E-009-011.