# The Closed-Loop Supply Chain Network Equilibrium with Products Lifetime and Carbon Emission Constraints in Multiperiod Planning Horizon

Authors: Guitao Zhang; Hao Sun; Jinsong Hu; Gengxin Dai
Journal: Discrete Dynamics in Nature and Society
 (2014)

---

## Abstract

This paper studies a closed-loop supply chain network equilibrium problem in multiperiod planning horizons with consideration of product lifetime and carbon emission constraints. The closed-loop supply chain network consists of suppliers tier, manufacturer tier, retailers tier, and demand markets tier, in which the manufacturers collect used products from the demand markets directly. Product lifetime is introduced to denote the maximum times of manufacturing and remanufacturing, and the relation between adjacent periods is described by inventory transfer. By variational inequalities and complementary theory, the optimal behaviors of all the players are modeled, and, in turn, the governing closed-loop supply chain network equilibrium model is established. The model is solved by modified project contraction algorithm with fixed step. Optimal equilibrium results are computed and analyzed through numerical examples. The impacts of collection rate, remanufacturing conversion rate, product lifetime, and carbon emission cap on equilibrium states are analyzed. Finally, several managerial insights are given to provide decision support for entrepreneurs and government official along with some inspirations for future research.


---

## Body

## 1. Introduction

Closed-loop supply chain (CLSC) has received increasing attention coming from theorists and entrepreneurs in recent years. From the viewpoint of overall lifecycle, it compromises the forward activities of manufacturing, distribution, and retailing and backward activities of collection and remanufacturing. Remanufacturing, as the significant part of CLSC, is a complicated process in which the core components of used products are disassembled, repaired, and reused [1]. Compared with production from virgin materials, remanufacturing tends to be energy saving, less material consuming and often has a lower impact on the environment. Related studies have shown that the remanufacturing cost of a core component is typically 40–60% less than those of brand-new production by raw materials [2]. Many famous enterprises, such as Xerox, Robert Bosch Tool, Black and Decker, and Hewlett-Packard, have implemented remanufacturing and CLSC strategies successfully [3, 4].In CLSC operations, the dynamic production, remanufacturing, and pricing strategies in multiperiod planning horizons are of great importance for the firms, which is directly related to their long-term profits. But for the decision-makers, it is much more complicated to make the optimal strategies in several planning horizons than in a single period, because they have to weigh the profit of the current period against those of the subsequent periods. One can easily deduce that if the price is high in the current period, the firm’s profit in this period might increase whereas the quantity of the products sold decreases; in turn, in the subsequent periods, the number of remanufacturable products available decreases, therefore reducing the subsequent periods’ potential profits. On the other hand, if the price is low in the current period, the profits in this period might decrease, but the manufacturers have more end-of-life products to remanufacture in future periods [5]. Then, in this context, how should the decision-makers make the best choice? Aiming at this issue, the scholars have done a lot of research from the viewpoint of a monopolistic manufacturer, duopoly manufacturers, and single manufacturer-single retailer/collector bilateral monopoly. Their research provides several meaningful management insights for the readers and the practitioners. However, as we know, the actual CLSC is likely to have a wider variety of channel members and complex cooperative and competitive relations. For example, a complete CLSC may have several mutual competitive suppliers, competitive manufacturers, competitive retailers, competitive collectors, and demand markets, in which each kind of players is sorted in the same tier, so the combination of any a player from each tier can form a cooperative chain. From another perspective, some players in the network such as a typical manufacturer or a typical retailer may be included in different CLSC. Therefore, the CLSC usually exists in the form of network structure. According to Hammond and Beullens [6], if and only if all the players in the CLSC network agree to the transaction prices and the transaction amounts of products, the equilibrium state of the CLSC network is achieved. Then, the following questions come up: for such a complex CLSC network, what are the players’ optimal decisions in dynamic multiperiod planning horizons? What is the dynamic equilibrium state of the CLSC network?In addition, the factors of product lifetime are usually neglected in the previous dynamic multiperiod CLSC modeling, and almost all of the literature assumed the used products could be remanufactured only once. But in reality, the core components may undergo multiple cycles of remanufacturing. For example, an end-of-life toner cartridge can be retreaded up to four times before landfill [7]. A single-use camera can be remanufactured approximately 6 times [8]. Many others including the maximum lives of car rites and computer chips are 3 and 4, respectively [9]. Therefore, it is of great importance to take the factors of finite product lifetime into account in dynamic CLSC modeling.Last but not least, as huge consumptions of fossil energies and accumulative carbon emissions of greenhouse gases have aggravated weather deterioration and environment pollution most recently, more and more attention on how to reduce carbon emissions has been received worldwide. Accordingly, a series of regulation policies are introduced by governments. Kyoto protocol is one of the most famous protocols signed at the third conference of the parties of United Nations Framework Convention on Climate Change held in Kyoto, Japan, in October, 1997. After Kyoto protocol, other policies, for example, EU-ETS, RGGI, are published subsequently [10].  Predictably, in the future the influence on enterprise operations will be more and more coming from carbon emission polices and gradually will be one of the prerequisites in the production/operations management for the firms. Therefore, for the players in CLSC network, how should they arrange their logistic activities when faced with the carbon emission policies?Motivated by all the above analysis, this paper establishes a multiperiod CLSC network equilibrium model with consideration of product lifetime and carbon emission constraints comprehensively, analyzes the impact of a variety of parameters on CLSC network equilibrium, and obtains several managerial insights. We argue that this study is one of the first studies to explore CLSC network equilibrium in dynamic multiperiod planning horizons, not to mention that the other two important factors are considered simultaneously.The rest of this paper is organized as follows. In Section2, we review the literatures related to our research and highlight the contributions of our paper. In Section 3, we give the model assumptions and variable notations. In Section 4, we present our model in which the optimal behavior of various players in the multiperiod CLSC network is described and the governing equilibrium conditions are formulated as a finite-dimension variational inequality problem. In Section 5, we provide a solving algorithm for the model. In addition, by a numerical example, we analyze the equilibrium results and reveal meaningful insights related to product lifetime and carbon emission constraints.


## 2. Literature Review

Our work is related to three streams of research, and the literature in each research will be reviewed subsequently. We will also point out how our study differs from the existing literature.The first stream of research related to our work is on the dynamic production, pricing policies, and coordination mechanisms of CLSC in two-period, multiperiod, and infinite planning horizon settings [5, 11–18]. With the consideration of multiple product lifetimes, Geyer et al. [9] investigated the economy of remanufacturing and demonstrated the need to coordinate production cost structure, collection rate, constraints of limited component durability, and finite product life cycles in the process of production operations. Although the above-mentioned literatures lay a solid foundation for future research, they are confined merely to the case of a monopoly manufacturer, duopoly manufacturers, or single manufacturer-single retailer/collector bilateral monopoly. Our study differs from the above-mentioned literature in that we research into dynamic production and pricing decisions in a complex CLSC network.The second stream of research related to our work is on network equilibrium of CLSC and multiperiod traditional/forward supply chain network equilibrium. Hammond and Beullens [6] construct an oligopolistic CLSC network model including manufacturers and demand markets under WEEE legislation. Yang et al. [19] use variational inequality method to model the CLSC network including suppliers, manufacturers who are involved in the production of a homogeneous product from virgin materials and collected materials, retailers, and recovery centers that collect used products from demand markets and can obtain subsidies from official organizations. Qiang et al. [20] establish a CLSC network model considering the competition, distribution channel investment, and uncertainties. The literature mentioned above deal with static or single-period CLSC network equilibrium problems.Recently, a few authors explore multiperiod traditional supply chain network equilibrium problems. Cruz and Wakolbinger [21] develop a framework for the analysis of the optimal levels of corporate social responsibility (CSR) activities in a multiperiod supply chain network consisting of manufacturers, retailers, and consumers and describe the problem of carbon emissions. Hamdouch [22] establishes a three-tier equilibrium model with capacity constraints and retailers’ purchase strategy from a multiperiod perspective. Cruz and Liu [23] analyze the effects of levels of social relationship on a multiperiod supply chain network with multiple decision-makers associated with different tiers. But to our knowledge, there is no research linking the dynamic multiperiod decisions and CLSC network equilibrium up to now. One intention of our work is to fill the gap.The third stream of research related to our work is on the production planning with carbon cap schemes and carbon trade schemes. Emission cap scheme, as a traditional regulation, restricts the maximum carbon emission volumes of the firm at a certain time. In contrast, emissions trading schemes provide pollutant emitters with flexibility in how they comply with the regulations. It is not discussed in detail in this paper. We refer readers to the work Gong and Zhou [24] to reach a deeper understanding about production planning with carbon emission schemes. Herein, we will give a brief review on the literature that studies supply chain network equilibrium and CLSC supply chain operations problems with carbon emission constraints. Wu et al. [25] propose an electric power supply chain network equilibrium model with carbon taxes. Based on this model, Nagurney et al. [26] develop a modeling and computational framework that allowed for the determination of optimal carbon taxes applied to electric power plants in electric power supply chain networks. More recently, it has been incorporated into reverse logistics and CLSC modeling. Krikke [27] gives a decision framework for optimizing the combined disposition and location-transport decisions on carbon footprints in CLSC and then applies it to a copier case. Kannan et al. [28] establish a biobjective 0-1 mixed-integer linear programming model for a carbon footprint based reverse logistics network design. The difference between our study and the existing literature lies in that we combine the CLSC network equilibrium with carbon emission constraints.


## 3. Model Assumptions and Notations

### 3.1. Model Assumptions

Considering a CLSC network consisting of multiple suppliers, multiple manufacturers, multiple retailers, and multiple demand markets, in which the manufacturers make homogenous products and simultaneously collected used products from demand markets at the end of periodt for remanufacturing in the period t+1, the criterion of each player in the network is its total profit maximum in multiperiod planning horizontal. The two-period CLSC network comprises two suppliers, two manufacturers, two retailers, and two demand markets which can be described as in Figure 1. We denote a typical raw material supplier by s, a typical manufacturer by m, a typical retailer market by n, and a typical demand market by k. In particular, m1(2) represents the 1st manufacturer in 2nd period, and we can interpret the other notations in the same way. In the top of Figure 1, the dotted lines between the same manufacturers in different period represent the transfer of inventory between different planning periods. The real lines express the forward transactions between suppliers and manufacturers and manufacturers and retailers, and the thin dotted lines illustrate the reverse transactions between manufacturers and demand markets.Figure 1
An illustration of a CLSC network with two periods.Due to the complexity of the problem, we make the following assumptions obeying the economic theories.Assumption 1.
The remanufactured products are the same as new products made by raw materials, and there is no difference between their selling prices in demand markets in each period. This assumption comes from the example of single-use cameras of Eastman Kodak Company. The used cameras are typically upgraded to the quality of new ones, and both products can be perfectly substituted by each other [29].Assumption 2.
Compared with the production by raw materials, the manufacturers firstly choose to remanufacture the collected products because of cost saving. The assumption is commonly used in CLSC management literature to ensure the economy of remanufacturing.Assumption 3.
It is assumed that the manufacturers have carbon emission constraints during manufacturing and remanufacturing in each period. It is different from the case of traditional/forward supply chain, in which the carbon constraints are usually applied for multiperiod of time, for example, one year. However, in CLSC, it will take a relatively long time from product design, production, retailing, collection, and remanufacturing. The average return lags alone are often at least half a year [30].Assumption 4.
The carbon emission volume of remanufacturing a used product is lower than that of manufacturing a new product [27].Assumption 5.
The manufacturers and demand markets are separate in space and the place, and the distance between different tiers is included in the transaction cost functions.Assumption 6.
The number of planning period is larger than product lifetime to illustrate the impact of product lifetime.Assumption 7.
One retailer only copes with one demand market.Assumption 8.
The production functions and transaction cost functions in the model are continuously differentiable and convex. This function property is required to be strictly satisfied in CLSC network equilibrium research [6, 19].


### 3.2. Variables and Notations

Consider the following:T: the number of total planning periods;Z: the product lifetime (Z≥2), which means that the product can be manufactured once and remanufactured Z-1 times;t: a typical planning period, t=1,2,…,T;s: a typical supplier, s=1,2,…,S;m: a typical manufacturer, m=1,2,…,M;n: a typical retailer, n=1,2,…,N;k: a typical demand market, k=1,2,…,K;α: the collection rate of used products from demand market k; in this paper, we assume that α is a constant parameter in different periods;βr: the average conversion rate of raw materials;βiu: the average remanufacturing conversion rate of used products in period i, i=1,…,min⁡{t-1,Z-1}, which in general satisfies βiu≤βi-1u, i=2,…,min⁡{t-1,Z-1}. Moreover, in order to simplify the formulation expressions hereinafter, we introduce β0u=1;qsmr(t): the amount of raw materials provided by supplier s to manufacturer m in period t; group all of qsmr(t) into a SMT-dimension column vector Q1;qsr(t): the amount of raw materials provided by supplier s to all manufacturers in period t; group all of qsr(t) into a ST-dimension column vector q1r and qsr(t)=∑m=1Mqsmr(t);qmr(t): the amount of raw materials that manufacturer m uses for production in period t; group all of qmr(t) into a MT-dimension column vector qr;qmu(t): the amount of collected products used in remanufacturing process in period t(t≥2) and qmu(t)=∑k=1Kqkmu(t-1)=∑n=1N∑i=1min⁡{t-1,Z-1}βi-1u⋯β0uαiqmnr(t-i); group all of qmu(t) into a M(T-1)-dimension column vector qu;qmn(t): the transaction amount of new products between manufacturer m and retailer n in period t group all of qmn(t) into a MNT-dimension column vector Q2;qmnr(t): the transaction amount of the new products made from raw materials between manufacturer m and retailer n in period t; group all of qmnr(t) into a MNT-dimension column vector Q3;qkmu(t): the transaction amount of used products between manufacturer m and demand market n in period t; group all of qkmu(t) into a NM-dimension column vector Q4(t), and group all of qkmu(t) into a NMT-dimension column vector Q4, and then we have qkmu(t-1)=∑i=1min⁡{t-1,Z-1}βi-1u⋯β0uαiqmnr(t-i);qnk(t): the transaction amount of products between manufacturer m and retailer n in period t; group all of qnk(t) into a NKT-dimension column vector Q5;fsr(q1r(t)): production cost of raw material by supplier s in period t;fmr(t)=fmr(βr,qr(t)): production cost of product by manufacturer m in period t;fmu(t)=fmu(βiu,qmnr(t-i)): remanufacturing cost of manufacturer m in period t, i=1,…,min⁡{t-1,Z-1};ckmu(t): transportation cost of used products undertaken by manufacturer m in period t(t≥2) from demand market k to manufacturer m;cm(qmu(t)): disposal cost of manufacturer m dealing with the collected products qmu(t) in period t(t≥2);Im(t): the inventory level of manufacturer m in period t; group all of Im(t) into a MT-dimension column vector IM;Hm(t)=Hm(Im(t)): the inventory cost function of manufacturer m in period t;ρsmr(t): the transaction price between supplier s and manufacturer m in period t, an endogenous variable;ρmn(t): the transaction price between manufacturer m and retailer n in period t, an endogenous variable;ρ-: unit disposal cost of used products by manufacturers, a constant;cmd(t)=cmd(βiu,qmr(t)): transportation cost generated from moving the wastes to landfill by manufacturer m in period t(t≥2);ρkmu(t): the transaction price of used products between manufacturer m and demand market k in period t(t≥2), an endogenous variable;csmr(t)=csmr(qsmr(t)): the transaction cost undertaken by the supplier between manufacturer m and supplier s in period t;cmn(t)=cmn(qmn(t)): the transaction cost undertaken by the manufacturer between manufacturer m and retailer n in period t;cn(qmn(t)): handling cost of retailer n, mainly referring to package and exhibition cost;cnkN(t)=cnkN(qnk(t)): the transaction cost between retailer n and demand market k in period t at retailer n;cnkK(t)=cnkK(qnk(t)): the transaction cost between retailer n and demand market k in period t at demand market k;αn(Q4(t)): the disutility function of consumer in demand market k when he sells the used products to manufacturers;ρnk(t): the price of unit product associated with retailer n and demand market k in period t, and all of ρnk(t) are the elements of the NT-dimension column vector ρ3;ρ4k(t): the purchasing price of unit product associated with demand market k in period t; all of ρ4k(t) in period t are the elements of the K-dimension column vector ρ4(t), and all of ρ4k(t) belong to KT-dimension column vector ρ4;dk(ρ4(t)): the demand volume of the products with the demand price ρ4k(t) at retailer n in period t, which is a monotonic decreasing function depending on ρ4k(t);ex: the carbon emissions when manufacturing and remanufacturing, in which x=r denotes the case of production by unit raw materials and x=u denotes the case of remanufacturing by used products; in general, eu<er is satisfied;Bm(t): carbon emission cap of manufacturer m in period t; group all of the Bm(t) into a MT-dimension column vector B;ssm(t): the carbon tax paid by manufacturer m in period t, an endogenous variable; group all of the ssm(t) into a MT-dimension column vector ss.


## 3.1. Model Assumptions

Considering a CLSC network consisting of multiple suppliers, multiple manufacturers, multiple retailers, and multiple demand markets, in which the manufacturers make homogenous products and simultaneously collected used products from demand markets at the end of periodt for remanufacturing in the period t+1, the criterion of each player in the network is its total profit maximum in multiperiod planning horizontal. The two-period CLSC network comprises two suppliers, two manufacturers, two retailers, and two demand markets which can be described as in Figure 1. We denote a typical raw material supplier by s, a typical manufacturer by m, a typical retailer market by n, and a typical demand market by k. In particular, m1(2) represents the 1st manufacturer in 2nd period, and we can interpret the other notations in the same way. In the top of Figure 1, the dotted lines between the same manufacturers in different period represent the transfer of inventory between different planning periods. The real lines express the forward transactions between suppliers and manufacturers and manufacturers and retailers, and the thin dotted lines illustrate the reverse transactions between manufacturers and demand markets.Figure 1
An illustration of a CLSC network with two periods.Due to the complexity of the problem, we make the following assumptions obeying the economic theories.Assumption 1.
The remanufactured products are the same as new products made by raw materials, and there is no difference between their selling prices in demand markets in each period. This assumption comes from the example of single-use cameras of Eastman Kodak Company. The used cameras are typically upgraded to the quality of new ones, and both products can be perfectly substituted by each other [29].Assumption 2.
Compared with the production by raw materials, the manufacturers firstly choose to remanufacture the collected products because of cost saving. The assumption is commonly used in CLSC management literature to ensure the economy of remanufacturing.Assumption 3.
It is assumed that the manufacturers have carbon emission constraints during manufacturing and remanufacturing in each period. It is different from the case of traditional/forward supply chain, in which the carbon constraints are usually applied for multiperiod of time, for example, one year. However, in CLSC, it will take a relatively long time from product design, production, retailing, collection, and remanufacturing. The average return lags alone are often at least half a year [30].Assumption 4.
The carbon emission volume of remanufacturing a used product is lower than that of manufacturing a new product [27].Assumption 5.
The manufacturers and demand markets are separate in space and the place, and the distance between different tiers is included in the transaction cost functions.Assumption 6.
The number of planning period is larger than product lifetime to illustrate the impact of product lifetime.Assumption 7.
One retailer only copes with one demand market.Assumption 8.
The production functions and transaction cost functions in the model are continuously differentiable and convex. This function property is required to be strictly satisfied in CLSC network equilibrium research [6, 19].


## 3.2. Variables and Notations

Consider the following:T: the number of total planning periods;Z: the product lifetime (Z≥2), which means that the product can be manufactured once and remanufactured Z-1 times;t: a typical planning period, t=1,2,…,T;s: a typical supplier, s=1,2,…,S;m: a typical manufacturer, m=1,2,…,M;n: a typical retailer, n=1,2,…,N;k: a typical demand market, k=1,2,…,K;α: the collection rate of used products from demand market k; in this paper, we assume that α is a constant parameter in different periods;βr: the average conversion rate of raw materials;βiu: the average remanufacturing conversion rate of used products in period i, i=1,…,min⁡{t-1,Z-1}, which in general satisfies βiu≤βi-1u, i=2,…,min⁡{t-1,Z-1}. Moreover, in order to simplify the formulation expressions hereinafter, we introduce β0u=1;qsmr(t): the amount of raw materials provided by supplier s to manufacturer m in period t; group all of qsmr(t) into a SMT-dimension column vector Q1;qsr(t): the amount of raw materials provided by supplier s to all manufacturers in period t; group all of qsr(t) into a ST-dimension column vector q1r and qsr(t)=∑m=1Mqsmr(t);qmr(t): the amount of raw materials that manufacturer m uses for production in period t; group all of qmr(t) into a MT-dimension column vector qr;qmu(t): the amount of collected products used in remanufacturing process in period t(t≥2) and qmu(t)=∑k=1Kqkmu(t-1)=∑n=1N∑i=1min⁡{t-1,Z-1}βi-1u⋯β0uαiqmnr(t-i); group all of qmu(t) into a M(T-1)-dimension column vector qu;qmn(t): the transaction amount of new products between manufacturer m and retailer n in period t group all of qmn(t) into a MNT-dimension column vector Q2;qmnr(t): the transaction amount of the new products made from raw materials between manufacturer m and retailer n in period t; group all of qmnr(t) into a MNT-dimension column vector Q3;qkmu(t): the transaction amount of used products between manufacturer m and demand market n in period t; group all of qkmu(t) into a NM-dimension column vector Q4(t), and group all of qkmu(t) into a NMT-dimension column vector Q4, and then we have qkmu(t-1)=∑i=1min⁡{t-1,Z-1}βi-1u⋯β0uαiqmnr(t-i);qnk(t): the transaction amount of products between manufacturer m and retailer n in period t; group all of qnk(t) into a NKT-dimension column vector Q5;fsr(q1r(t)): production cost of raw material by supplier s in period t;fmr(t)=fmr(βr,qr(t)): production cost of product by manufacturer m in period t;fmu(t)=fmu(βiu,qmnr(t-i)): remanufacturing cost of manufacturer m in period t, i=1,…,min⁡{t-1,Z-1};ckmu(t): transportation cost of used products undertaken by manufacturer m in period t(t≥2) from demand market k to manufacturer m;cm(qmu(t)): disposal cost of manufacturer m dealing with the collected products qmu(t) in period t(t≥2);Im(t): the inventory level of manufacturer m in period t; group all of Im(t) into a MT-dimension column vector IM;Hm(t)=Hm(Im(t)): the inventory cost function of manufacturer m in period t;ρsmr(t): the transaction price between supplier s and manufacturer m in period t, an endogenous variable;ρmn(t): the transaction price between manufacturer m and retailer n in period t, an endogenous variable;ρ-: unit disposal cost of used products by manufacturers, a constant;cmd(t)=cmd(βiu,qmr(t)): transportation cost generated from moving the wastes to landfill by manufacturer m in period t(t≥2);ρkmu(t): the transaction price of used products between manufacturer m and demand market k in period t(t≥2), an endogenous variable;csmr(t)=csmr(qsmr(t)): the transaction cost undertaken by the supplier between manufacturer m and supplier s in period t;cmn(t)=cmn(qmn(t)): the transaction cost undertaken by the manufacturer between manufacturer m and retailer n in period t;cn(qmn(t)): handling cost of retailer n, mainly referring to package and exhibition cost;cnkN(t)=cnkN(qnk(t)): the transaction cost between retailer n and demand market k in period t at retailer n;cnkK(t)=cnkK(qnk(t)): the transaction cost between retailer n and demand market k in period t at demand market k;αn(Q4(t)): the disutility function of consumer in demand market k when he sells the used products to manufacturers;ρnk(t): the price of unit product associated with retailer n and demand market k in period t, and all of ρnk(t) are the elements of the NT-dimension column vector ρ3;ρ4k(t): the purchasing price of unit product associated with demand market k in period t; all of ρ4k(t) in period t are the elements of the K-dimension column vector ρ4(t), and all of ρ4k(t) belong to KT-dimension column vector ρ4;dk(ρ4(t)): the demand volume of the products with the demand price ρ4k(t) at retailer n in period t, which is a monotonic decreasing function depending on ρ4k(t);ex: the carbon emissions when manufacturing and remanufacturing, in which x=r denotes the case of production by unit raw materials and x=u denotes the case of remanufacturing by used products; in general, eu<er is satisfied;Bm(t): carbon emission cap of manufacturer m in period t; group all of the Bm(t) into a MT-dimension column vector B;ssm(t): the carbon tax paid by manufacturer m in period t, an endogenous variable; group all of the ssm(t) into a MT-dimension column vector ss.


## 4. The Multiperiod CLSC Network Model

### 4.1. The Optimal Behavior of Suppliers

The suppliers supplies raw material to the manufactures in each period and seeks for the profit maximization in the planning horizons. Based on the notations described above, we can express the criterion for raw material supplier s as follows:
(1)πs(qsmr*(t),qsr*(t))=max⁡{∑t=1T∑m=1Mρsmr(t)qsmr(t)-∑t=1T∑m=1Mcsmr(t)hhhhhhh-∑t=1Tfsr(q1r(t))}(2)s.t.∑m=1Mqsmr(t)≤qsr(t)(3)(qsmr(t),qsr(t))∈R+SM+S,∀s,m.The constraint (2) expresses that production amount of raw materials by supplier s is higher than that of the transaction between the supplier s and all the manufacturers.It is assumed that all of the suppliers compete in a noncooperative fashion. Therefore, we can simultaneously express the optimal conditions of the suppliers as the variational inequality: determine(Q1*,q1r*,ηs)∈R+SMT+2ST, satisfying
(4)∑t=1T∑s=1S∑m=1M[∂csmr(t)∂qsmr(t)-ρsmr(t)+ηs(t)]hhhhhhh×[qsmr(t)-qsmr*(t)]+∑t=1T∑s=1S[∂fsr(qsr(t))∂qsr(t)-ηs(t)]hhhhhhh×[qsr(t)-qsr*(t)]+∑t=1T∑s=1S[qsr(t)-∑m=1Mqsmr(t)]×[ηs(t)-ηs*(t)]≥0,hhhhhhhhhhhhhhhhhhhhh∀(Q1,q1r,ηs)∈R+SMT+2ST.In (4), ηs(t) is the Lagrange multiplier corresponding to constraint (2) and ηs is the ST-dimension column vector with the elements of ηs(t).


### 4.2. The Optimal Behavior of Manufacturers

The manufacturerm purchases the raw material from suppliers to make products, sells his new products to the retailers in every period, and manages inventory between different periods according to market demand; in the same time, the manufacturer m collects the used products consumed by demand markets at the end of period t, which will be remanufactured in the next period t+1.In the first planning period, the manufacturers just make products with virgin materials; from the second period, the manufacturers can make products with virgin materials and collected products simultaneously. With the consideration of the lifetimeZ, the collected products can be remanufactured at most Z-1 times. When their remanufactured times arrive Z-1, they have no value anymore and must be disposed in an appropriate way.For convenience, letβi=βiuβi-1u⋯β1u, i=1,…,min⁡{t-1,Z-1}. The profit function of manufacturer m can be stated as follows:
(5)πm(qsmr(t),qmr(t),qmu(t),qmn(t),qmnr(t),qnmu(t),Im(t))=max⁡{∑t=1T∑n=1Nρmn(t)qmn(t)-∑t=1Tfmr(t)-∑t=2Tfmu(t)hhhhhhh-∑t=2Tcm(qmu(t))-∑t=1T∑k=1Kcnmu(t)-∑t=1THm(t)hhhhhhh-∑t=1T∑m=1Mρsmr(t)qsmr(t)hhhhhhh-∑t=1T∑k=1Kρnmu(t)qnmu(t)-cmn(t)hhhhhhh-ρ-∑t=2T∑i=1min⁡⁡{t-1,Z-1}βi-1(1-βiu)αiqmnr(t-i)hhhhhhh-∑t=2T∑i=1min⁡⁡{t-1,Z-1}cmd(t)hhhhhhh-er∑t=1Tsm(t)qmr(t)-eu∑t=2Tsm(t)qmu(t)}(6)s.t.Im(t)+∑n=1Nqmn(t)=∑n=1N∑i=1min⁡⁡(t-1,Z-1)βiαiqmnr(t-i)hhhhhhhhhhhhhhhhhhhhhhhhhhhh+βrqmr+Im(t-1)(7)qmu(t)=∑n=1Nqnmu(t-1),∀t≥2(8)∑n=1Nqnmu(t-1)=∑n=1N∑i=1min⁡{t-1,Z-1}βi-1αiqmnr(t-i)∀t≥2;(9)∑n=1Nqmn(t)=∑n=1Nqmnr(t)+∑n=1N∑i=1min⁡{t-1,Z-1}βiαiqmnr(t-i)(10)qmr(t)≤∑s=1Sqsmr(t)(11)qsmr(t),qmr(t),qmu(t),qmn(t),qmnr(t),qnmu(t),Im(t)≥0,hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhihh∀m,s,n.Constraint (6) describes that the current inventory plus the products amount selling to all retailers is equal to the sum of the products amount made from used products, virgin materials, and the inventory from the prior period. Constraint (7) expresses that the manufacturer m collects the used products in period t-1 and remanufactures them in period t. The constraint (8) means that all the remanufacturable used products (in other words, their remanufacturing times are less than Z-1) are collected from the consumers by the manufacturer m at the end of period t-1. Constraint (9) expresses that, apart from remanufacturing the used products, the manufacturers still need to produce with raw materials in each period.Similar to [6, 20], all the manufactures compete in a noncooperative fashion, and each manufacturer seeks to maximize his profit based on the other manufacturers’ optimal decisions, that is, deciding the amount of raw materials purchased from the suppliers, the used products collected from the consumers, the products sold to retailers, and the inventory transferring to next period in every period. Therefore, the optimal conditions of all the manufacturers can be described by the following variational inequality: determine (qr*,qu*,Q1*,Q2*,Q3*,Q4*,IM*,γ1*,γ2*,γ3*,γ4*,γ5*)∈ΩM, satisfying
(12)∑t=1T∑m=1M∑s=1S[ρsmr*(t)-γ5m*(t)]×[qsmr(t)-qsmr*(t)]+∑t=1T∑m=1M[∂fmr(t)∂qmr(t)-βrγ1m*(t)+γ5m*(t)+ersm*(t)]hhhhhhhh×[qmr(t)-qmr*(t)]+∑t=2T∑m=1M[∂cm(qmu*(t))∂qmu(t)+γ2m*(t-1)+eusm*(t)]hhhhhhhh×[qmu(t)-qmu*(t)]+∑t=1T∑m=1M∑n=1N[∂cmn(t)∂qmn(t)+γ1m*(t)+γ4m*(t)-ρmn*(t)]hhhhhhhhhh×[qmn(t)-qmn*(t)]+∑t=1T∑n=1N∑m=1M[∂cnmu*(t)∂qnmu(t)+ρnmu*(t)-γ2m*(t)-γ3m*(t)]hhhhhhhhhh×[qnmu(t)-qnmu*(t)]+∑t=1T∑n=1N∑m=1M[∑i=1min⁡⁡{t,N-1}[∂fmu*(t+i-1)∂qmnr(t)+cmd(t-i)∂qmnr(t)]hhhhhhhhhhhhh-∑i=0min⁡⁡{t+1,Z-1}βiαiγ4m*(t+i)hhhhhhhhhhhhh-∑i=1min⁡⁡{t+1,Z-1}βiαiγ1m*(t+i)+ρ-hhhhhhhhhhhhh×∑i=1min⁡⁡{t+1,Z-1}βi-1(1-βiu)αihhhhhhhhhhhhh+∑i=1min⁡⁡{t+1,Z-1}βi-1αiγ3m*(t+i-1)]hhhhhhhhhhh×[qmnr(t)-qmnr*(t)]+∑t=1T∑m=1M[∂Hm*(t)∂Im(t)+γ1m*(t)-γ1m*(t+1)]hhhhhhhh×[Im(t)-Im*(t)]+∑t=1T∑m=1M[∑n=1N∑i=1min⁡⁡{t-1,Z-1}βiαiqmnr*(t-i)hhhhhhhhhhhhhhhhhhh+βrqmr*(t)+Im*(t-1)hhhhhhhhhhhhhhhhhhh-∑n=1Nqmn*(t)-Im*(t)]hhhhhhhh×[γ1m(t)-γ1m*(t)]+∑t=1T∑m=1M[∑n=1Nqnmu*(t)-qmu*(t+1)]hhhhhhhh×[γ2m(t)-γ2m*(t)]+∑t=2T∑m=1M[∑n=1Nqnmu*(t-1)hhhhhhhhhh-∑n=1N∑i=1min⁡⁡{t-1,Z-1}βi-1αiqmnr*(t-i)]hhhhhhhh×[γ3m(t)-γ3m*(t)]+∑t=1T∑m=1M[∑n=1Nqmnr*(t)+∑n=1N∑i=1min⁡⁡{t-1,Z-1}βiαiqmnr*(t-i)hhhhhhhhhh-∑n=1Nqmn*(t)]×[γ4m(t)-γ4m*(t)]+∑t=1T∑m=1M[∑s=1Sqsmr*(t)-qmr*(t)]×[γ5m(t)-γ5m*(t)]≥0∀(qr,qu,Q1,Q2,Q3,Q4,IM,γ1,γ2,γ3,γ4,γ5)∈ΩM,
where ΩM={(qr,qu,Q1,Q2,Q3,Q4,IM,γ1,γ2,γ3,γ4,γ5)∈R+MT+M(T-1)+MST+2MNT+KMT+MT×R4MT×R+MT}.In (12), γ1m(t), γ2m(t), γ3m(t), γ4m(t), and γ5m(t) are the Lagrange multipliers corresponding to constraints (6), (7), (8), (9), and (10), and γ1, γ2, γ3, γ4, and γ5 are the MT-dimension column vector with the elements of γ1m(t), γ2m(t), γ3m(t), γ4m(t), and γ5m(t), respectively.From the forth term of (12), we can find that when the transaction amount qmn*(t)>0, ρmn*(t)=(∂cmn*(t)/∂qmn(t))+γ1m*(t)+γ4m*(t).In the model, we consider the manufacturers’ carbon emission constraints from the government in each period, and then the manufacturers’ behaviors can be characterized by the following equation:(13)Bm(t)-(erqmr(t)+euqmu(t)){≥0,ssm(t)=0,=0,ssm(t)>0.Equation (13) describes that if the manufactures’ carbon emission amount is lower than Bm(t) in period t, the carbon tax is zero; on the other hand, if the manufactures’ carbon emission amount is equal to Bm(t) in period t, the manufacturers must pay taxes.The above equilibrium condition corresponds to the following variational inequality: determiness*∈R+MT, satisfying
(14)∑t=1T∑m=1M[Bm*(t)-(erqmr*(t)+euqmu*(t))]×[ssm(t)-ssm*(t)]≥0,∀ss∈R+MT.In a word, the manufacturers’ optimal behaviors must meet (10) and (14) simultaneously.


### 4.3. The Optimal Behavior of Retailers

The retailers need to decide how many products to purchase from every manufacturer and to sell to the consumers in demand markets to meet the random demand in markets.For retailern, the profits can be defined as follows:
(15)πn(Q2,Q5)=max⁡(∑t=1Tρnk(t)qnk(t)hhhhhhhhh-∑t=1Tcn(t)-∑t=1T∑m=1Mρmn(t)qmn(t)hhhhhhhhh-∑t=1TcnkN(t))(16)s.t.qnk(t)≤∑m=1Mqmn(t)(17)qmn(t),qnk(t)≥0,∀m,n,k.We also assume that the retailers compete in a noncooperative manner. The optimal behavior of each retailer is to maximize his profit when the other retailers’ behaviors are given. So, the optimal conditions of all retailers can be described by variational inequality: determine(Q2*,Q5*,μ*)∈ΩN, satisfying
(18)∑t=1T∑n=1N[cnkN*(t)∂qnk(t)-ρnk*(t)+μn*(t)]×[qnk(t)-qnk*(t)]+∑t=1T∑n=1N∑m=1M[ρmn*(t)+∂cn*(t)∂qmn-μn*(t)]hhhhhhhhhhh×[qmn(t)-qmn*(t)]+∑t=1T∑n=1N[∑m=1Mqmn(t)-qnk(t)]×[μn(t)-μn*(t)]≥0hhhhhhhhhhhhhhhhhhhhhhhhhhh∀(Q2,Q5,μ)∈ΩN,
where ΩN={(Q2,Q5,μ)∈R+MNT+NKT+NT}.In (18), μn(t) is the Lagrange multiplier corresponding to constraints (16) and μ is the NT-dimension column vector with the elements of μn(t).From the first term of (18), we find that when the transaction volume qnk(t)>0, ρnk*(t)=(cnkN*(t)/∂qnk(t))+μn*(t).


### 4.4. The Optimal Behavior of Demand Markets

In the forward supply chain, the consumers of each demand market decide to purchase the products from the retailer or not based on the appreciation level of products and also decide the purchasing price and amount.According to literature [19], the transaction between the demand market and the retailer should satisfy a complement condition; that is, for some demand market k,
(19)dk(ρ4(t)){=∑n=1Nqnk(t),ρ4k(t)>0,≤∑n=1Nqnk(t),ρ4k(t)=0.Equation (19) illustrates that if the product amount supplied by each retailer is equal to that of demand in markets, then the transaction can proceed and the transaction price is positive. The demand amount depends on the transaction price paid by consumers in demand market k.And also similar to [19], in the reverse supply chain, the demand markets’ behavior should be characterized by the following equation:
(20)αk(Q4*(t)){=ρkmu*(t),qkmu(t)>0,≥ρkmu*(t),qkmu(t)=0,
where αn(Q4(t)) is the disutility of consumers in the demand market n at period t. Equation (20) expresses that the higher collection price paid by manufacturers, the more used products the manufacturers collect. Moreover, the collected amount is lower than or equal to that of the products sold to the demand market, so the following constraint should be satisfied:
(21)∑m=1Mqkmu(t)≤∑n=1Nqnk(t).Letλk(t) denote the Lagrange multiplier corresponding to (21), group all of λk(t) into a KT-dimension column vector λ.Besides the purchasing price paying to the retailers, the consumers must pay the transaction cost to obtain the products; hence, we have(22)ρnk*(t)+cnkK*(t){=ρ4k(t),qnk*(t)>0,≥ρ4k(t),qnk*(t)=0.Equations (19), (20), (21), and (22) can be described by the following variational inequality: determine (Q4*,Q5*,ρ4*,λ*)∈ΩK, satisfying
(23)∑t=1T∑n=1N∑m=1M[αn(Q4*(t))-ρkmu*(t)+λk*(t)]hhhhhh×[qkmu(t)-qkmu*(t)]+∑t=1T∑n=1N∑k=1K[ρnk*(t)+cnkK*(t)-λk*(t)]hhhhhhhhh×[qnk(t)-qnk*(t)]+∑t=1T∑k=1K[∑n=1Nqnk*(t)-dk(ρ4*(t))]hhhhhhh×[ρ4k(t)-ρ4k*(t)]+∑t=1T∑k=1K[∑n=1Nqnk*(t)-∑m=1Mqkmu*(t)]hhhhhhh×[λk(t)-λk*(t)]≥0,hhhhhhhh∀(Q4,Q5,ρ,λ)∈ΩK,
where ΩK={(Q4,Q5,ρ4,λ)∈R+NMT+NKT+2NT}.From the first term of inequality (23), we can conclude that when the transaction amount of collected products qnmu(t)≥0, the transaction price of collected products between manufacturers and demand markets ρnmu*(t)=αn(Q4*(t))+λn*(t). In fact, since ∑m=1Mqmn(t)>∑m=1Mqnmu(t) and λn*(t)=0, we can get ρnmu*(t)=αn(Q4*(t)).


### 4.5. The Equilibrium Conditions of Multiperiod CLSC Network

Theorem 1.
In equilibrium of the multiperiod CLSC network, the suppliers’ optimal behaviors, as expressed by variational inequality (4), the manufacturers’ optimal behaviors, as expressed by variational inequality (12) and (14), the retailers’ optimal behaviors, as expressed by variational inequality (18), and the consumers’ optimal behaviors, as expressed by variational inequality (23), must be satisfied simultaneously. One can state it in the following variational inequality: determine (q1r*,qr*,qu*,Q1*,Q2*,Q3*,Q4*,IM*,ss*,ρ*,ηs*,γ1*,γ2*,γ3*,γ4*,γ5*,μ*,λ*)∈Ω, satisfying
(24)∑t=1T∑s=1S∑m=1M[∂csmr*(t)∂qsmr(t)+ηs*(t)-γ5m*(t)]hhhhhh×[qsmr(t)-qsmr*(t)]+∑t=1T∑s=1S[∂fsr*(qsr(t))∂qsr(t)-ηs*(t)]hhhhhhh×[qsr(t)-qsr*(t)]+∑t=1T∑m=1M[∂fmr*(t)∂qmr(t)-βrγ1m*(t)+ersm*(t)+γ5m*(t)]hhhhhhhh×[qmr(t)-qmr*(t)]+∑t=2T∑m=1M[∂cm(qmu*(t))∂qmu(t)+γ2m*(t-1)+eusm*(t)]hhhhhhhh×[qmu(t)-qmu*(t)]+∑t=1T∑m=1M∑n=1N[∂cmn*(t)qmn(t)+∂cn*(t)∂qmnhhhhhhhhihhhhh+γ1m*(t)+γ4m*(t)-λk*(t)-μn*(t)∂cmn*(t)qmn(t)+∂cn*(t)∂qmn]hhhhhhhhihh×[qmn(t)-qmn*(t)]+∑t=1T∑k=1K∑m=1M[∂ckmu*(t)∂qkmu(t)-γ2m*(t)-γ3m*(t)hhhhhhhhhhhhh+αn(Q4*(t))+λk*(t)∂ckmu*(t)∂qkmu(t)]hhhhhhhhihh×[qkmu(t)-qkmu*(t)]+∑t=1T∑n=1N∑m=1M[∑i=1min⁡⁡{t,N-1}[∂fmu*(t+i-1)∂qmnr(t)+cmd(t-i)∂qmnr(t)]hhhhhhhhhhhhhhh-∑i=0min⁡⁡{t+1,Z-1}βiαiγ4m*(t+i)hhhhhhhhhhhhhhh-∑i=1min⁡⁡{t+1,Z-1}βiαiγ1m*(t+i)hhhhhhhhhhhhhhh+ρ-∑i=1min⁡⁡{t+1,Z-1}βi-1(1-βiu)αihhhhhhhhhhhhhhh+∑i=1min⁡⁡{t+1,Z-1}βi-1αiγ3m*(t+i-1)]hhhhhhhhihhh×[qmnr(t)-qmnr*(t)]+∑t=1T∑m=1M[∂Hm*(t)∂Im(t)+γ1m*(t)-γ1m*(t+1)]hhhhhhhh×[Im(t)-Im*(t)]+∑t=1T∑m=1M[Bm(t)-(erqmr*(t)+euqmu*(t))]hhhhhhhh×[ssm(t)-ssm*(t)]+∑t=1T∑n=1N[cnkN*(t)∂qnk(t)+cnkK*(t)-λk*(t)+μn*(t)]hhhhh×[qnk(t)-qnk*(t)]+∑t=1T∑k=1K[∑n=1Nqnk*(t)-dk(ρ4*(t))]hhhhhhhh×[ρ4k(t)-ρ4k*(t)]+∑t=1T∑s=1S[qsr*(t)-∑m=1Mqsmr*(t)]hhhhhhhh×[ηs(t)-ηs*(t)]+∑t=1T∑m=1M[∑n=1N∑i=1min⁡⁡{t-1,Z-1}βiαiqmnr*(t-i)+βrqmr*(t)hhhhhhhhhhh+Im*(t-1)-∑n=1Nqmn*(t)-Im*(t)]hhhhhhhh×[γ1m(t)-γ1m*(t)]+∑t=1T∑m=1M[∑n=1Nqnmu*(t)-qmu*(t+1)]hhhhhhhh×[γ2m(t)-γ2m*(t)]+∑t=2T∑m=1M[∑n=1Nqnmu*(t-1)hhhhhhhhhh-∑n=1N∑i=1min⁡⁡{t-1,Z-1}βi-1αiqmnr*(t-i)]hhhhhhhh×[γ3m(t)-γ3m*(t)]+∑t=1T∑m=1M[∑n=1Nqmnr*(t)hhhhhhhhhhh+∑n=1N∑i=1min⁡⁡{t-1,Z-1}βiαiqmnr*(t-i)-∑n=1Nqmn*(t)]hhhhhhhh×[γ4m(t)-γ4m*(t)]+∑t=1T∑m=1M[∑s=1Sqsmr*(t)-qmr*(t)]hhhhhhhh×[γ5m(t)-γ5m*(t)]+∑t=1T∑n=1N[∑m=1Mqmn(t)-qnk(t)]hhhhhhh×[μn(t)-μn*(t)]+∑t=1T∑k=1K[qnk*(t)-∑m=1Mqkmu*(t)]hhhhhhh×[λk(t)-λk*(t)]≥0,∀(q1r,qr,qu,Q1,Q2,Q3,Q4,IM,ss,hhv,ρ,ηs,γ1,γ2,γ3,γ4,γ5,μ,λ)∈Ω,
where Ω=ΩS×ΩM×ΩN×ΩK.


### 4.6. Solving Algorithm and Parameter Setting

For easy formulation hereinafter, we group the terms of the multiplication signs in inequality (24) into a column vector vector F(y)={Fsmt1,Fst2,Fmt3,Fm(t-1)4,Fmnt5,Fmnt6,Fkmt7,Fmt8,Fmt9,Fnt10,Fnt11,Fst12,Fmt13,Fmt14,Fmt15,Fmt16,Fmt17,Fnt18,Fnt19}∀s,m,n,k,t and introduce y=(q1r,qr,qu,Q1,Q2,Q3,Q4,IM,ss,v,ρ,ηs,γ1,γ2,γ3,γ4,γ5,μ,λ)∈Ω; thus, we can rewrite the variational inequality problem (24) in the standard form: determine y*=(q1r*,qr*,qu*,Q1*,Q2*,Q3*,Q4*,IM*,ss*,v*,ρ*,ηs*,γ1*,γ2*,γ3*,γ4*,γ5*,μ*,λ*), satisfying
(25)F(y)T·(y-y*)≥0,∀y∈Ω.In order to solve the variational inequality (25), we have several algorithms to select, for example, the modified project contraction algorithm [31], logarithmic-quadratic proximal prediction-correction method [32]. In this paper, we employ the modified project contraction algorithm to solve the problem for two reasons. Firstly, the modified project contraction algorithm provides an iterative solution framework that can deal with variational inequalities defined on polyhedral sets and get the optimal solutions and Lagrangian multipliers simultaneously. So, it satisfies all the requirements of a solution method for our model. Secondly, the design of its iterative steps is much simpler than that of logarithmic-quadratic proximal prediction-correction method. Set the related parameters of the modified project contraction algorithm as follows: the original value of decision variables and Lagrange multipliers is set to 1; the convergence criterion, for example, the absolute value of difference of decision variables and Lagrange multipliers between two steps, is lower than or equal to 10−8.


## 4.1. The Optimal Behavior of Suppliers

The suppliers supplies raw material to the manufactures in each period and seeks for the profit maximization in the planning horizons. Based on the notations described above, we can express the criterion for raw material supplier s as follows:
(1)πs(qsmr*(t),qsr*(t))=max⁡{∑t=1T∑m=1Mρsmr(t)qsmr(t)-∑t=1T∑m=1Mcsmr(t)hhhhhhh-∑t=1Tfsr(q1r(t))}(2)s.t.∑m=1Mqsmr(t)≤qsr(t)(3)(qsmr(t),qsr(t))∈R+SM+S,∀s,m.The constraint (2) expresses that production amount of raw materials by supplier s is higher than that of the transaction between the supplier s and all the manufacturers.It is assumed that all of the suppliers compete in a noncooperative fashion. Therefore, we can simultaneously express the optimal conditions of the suppliers as the variational inequality: determine(Q1*,q1r*,ηs)∈R+SMT+2ST, satisfying
(4)∑t=1T∑s=1S∑m=1M[∂csmr(t)∂qsmr(t)-ρsmr(t)+ηs(t)]hhhhhhh×[qsmr(t)-qsmr*(t)]+∑t=1T∑s=1S[∂fsr(qsr(t))∂qsr(t)-ηs(t)]hhhhhhh×[qsr(t)-qsr*(t)]+∑t=1T∑s=1S[qsr(t)-∑m=1Mqsmr(t)]×[ηs(t)-ηs*(t)]≥0,hhhhhhhhhhhhhhhhhhhhh∀(Q1,q1r,ηs)∈R+SMT+2ST.In (4), ηs(t) is the Lagrange multiplier corresponding to constraint (2) and ηs is the ST-dimension column vector with the elements of ηs(t).


## 4.2. The Optimal Behavior of Manufacturers

The manufacturerm purchases the raw material from suppliers to make products, sells his new products to the retailers in every period, and manages inventory between different periods according to market demand; in the same time, the manufacturer m collects the used products consumed by demand markets at the end of period t, which will be remanufactured in the next period t+1.In the first planning period, the manufacturers just make products with virgin materials; from the second period, the manufacturers can make products with virgin materials and collected products simultaneously. With the consideration of the lifetimeZ, the collected products can be remanufactured at most Z-1 times. When their remanufactured times arrive Z-1, they have no value anymore and must be disposed in an appropriate way.For convenience, letβi=βiuβi-1u⋯β1u, i=1,…,min⁡{t-1,Z-1}. The profit function of manufacturer m can be stated as follows:
(5)πm(qsmr(t),qmr(t),qmu(t),qmn(t),qmnr(t),qnmu(t),Im(t))=max⁡{∑t=1T∑n=1Nρmn(t)qmn(t)-∑t=1Tfmr(t)-∑t=2Tfmu(t)hhhhhhh-∑t=2Tcm(qmu(t))-∑t=1T∑k=1Kcnmu(t)-∑t=1THm(t)hhhhhhh-∑t=1T∑m=1Mρsmr(t)qsmr(t)hhhhhhh-∑t=1T∑k=1Kρnmu(t)qnmu(t)-cmn(t)hhhhhhh-ρ-∑t=2T∑i=1min⁡⁡{t-1,Z-1}βi-1(1-βiu)αiqmnr(t-i)hhhhhhh-∑t=2T∑i=1min⁡⁡{t-1,Z-1}cmd(t)hhhhhhh-er∑t=1Tsm(t)qmr(t)-eu∑t=2Tsm(t)qmu(t)}(6)s.t.Im(t)+∑n=1Nqmn(t)=∑n=1N∑i=1min⁡⁡(t-1,Z-1)βiαiqmnr(t-i)hhhhhhhhhhhhhhhhhhhhhhhhhhhh+βrqmr+Im(t-1)(7)qmu(t)=∑n=1Nqnmu(t-1),∀t≥2(8)∑n=1Nqnmu(t-1)=∑n=1N∑i=1min⁡{t-1,Z-1}βi-1αiqmnr(t-i)∀t≥2;(9)∑n=1Nqmn(t)=∑n=1Nqmnr(t)+∑n=1N∑i=1min⁡{t-1,Z-1}βiαiqmnr(t-i)(10)qmr(t)≤∑s=1Sqsmr(t)(11)qsmr(t),qmr(t),qmu(t),qmn(t),qmnr(t),qnmu(t),Im(t)≥0,hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhihh∀m,s,n.Constraint (6) describes that the current inventory plus the products amount selling to all retailers is equal to the sum of the products amount made from used products, virgin materials, and the inventory from the prior period. Constraint (7) expresses that the manufacturer m collects the used products in period t-1 and remanufactures them in period t. The constraint (8) means that all the remanufacturable used products (in other words, their remanufacturing times are less than Z-1) are collected from the consumers by the manufacturer m at the end of period t-1. Constraint (9) expresses that, apart from remanufacturing the used products, the manufacturers still need to produce with raw materials in each period.Similar to [6, 20], all the manufactures compete in a noncooperative fashion, and each manufacturer seeks to maximize his profit based on the other manufacturers’ optimal decisions, that is, deciding the amount of raw materials purchased from the suppliers, the used products collected from the consumers, the products sold to retailers, and the inventory transferring to next period in every period. Therefore, the optimal conditions of all the manufacturers can be described by the following variational inequality: determine (qr*,qu*,Q1*,Q2*,Q3*,Q4*,IM*,γ1*,γ2*,γ3*,γ4*,γ5*)∈ΩM, satisfying
(12)∑t=1T∑m=1M∑s=1S[ρsmr*(t)-γ5m*(t)]×[qsmr(t)-qsmr*(t)]+∑t=1T∑m=1M[∂fmr(t)∂qmr(t)-βrγ1m*(t)+γ5m*(t)+ersm*(t)]hhhhhhhh×[qmr(t)-qmr*(t)]+∑t=2T∑m=1M[∂cm(qmu*(t))∂qmu(t)+γ2m*(t-1)+eusm*(t)]hhhhhhhh×[qmu(t)-qmu*(t)]+∑t=1T∑m=1M∑n=1N[∂cmn(t)∂qmn(t)+γ1m*(t)+γ4m*(t)-ρmn*(t)]hhhhhhhhhh×[qmn(t)-qmn*(t)]+∑t=1T∑n=1N∑m=1M[∂cnmu*(t)∂qnmu(t)+ρnmu*(t)-γ2m*(t)-γ3m*(t)]hhhhhhhhhh×[qnmu(t)-qnmu*(t)]+∑t=1T∑n=1N∑m=1M[∑i=1min⁡⁡{t,N-1}[∂fmu*(t+i-1)∂qmnr(t)+cmd(t-i)∂qmnr(t)]hhhhhhhhhhhhh-∑i=0min⁡⁡{t+1,Z-1}βiαiγ4m*(t+i)hhhhhhhhhhhhh-∑i=1min⁡⁡{t+1,Z-1}βiαiγ1m*(t+i)+ρ-hhhhhhhhhhhhh×∑i=1min⁡⁡{t+1,Z-1}βi-1(1-βiu)αihhhhhhhhhhhhh+∑i=1min⁡⁡{t+1,Z-1}βi-1αiγ3m*(t+i-1)]hhhhhhhhhhh×[qmnr(t)-qmnr*(t)]+∑t=1T∑m=1M[∂Hm*(t)∂Im(t)+γ1m*(t)-γ1m*(t+1)]hhhhhhhh×[Im(t)-Im*(t)]+∑t=1T∑m=1M[∑n=1N∑i=1min⁡⁡{t-1,Z-1}βiαiqmnr*(t-i)hhhhhhhhhhhhhhhhhhh+βrqmr*(t)+Im*(t-1)hhhhhhhhhhhhhhhhhhh-∑n=1Nqmn*(t)-Im*(t)]hhhhhhhh×[γ1m(t)-γ1m*(t)]+∑t=1T∑m=1M[∑n=1Nqnmu*(t)-qmu*(t+1)]hhhhhhhh×[γ2m(t)-γ2m*(t)]+∑t=2T∑m=1M[∑n=1Nqnmu*(t-1)hhhhhhhhhh-∑n=1N∑i=1min⁡⁡{t-1,Z-1}βi-1αiqmnr*(t-i)]hhhhhhhh×[γ3m(t)-γ3m*(t)]+∑t=1T∑m=1M[∑n=1Nqmnr*(t)+∑n=1N∑i=1min⁡⁡{t-1,Z-1}βiαiqmnr*(t-i)hhhhhhhhhh-∑n=1Nqmn*(t)]×[γ4m(t)-γ4m*(t)]+∑t=1T∑m=1M[∑s=1Sqsmr*(t)-qmr*(t)]×[γ5m(t)-γ5m*(t)]≥0∀(qr,qu,Q1,Q2,Q3,Q4,IM,γ1,γ2,γ3,γ4,γ5)∈ΩM,
where ΩM={(qr,qu,Q1,Q2,Q3,Q4,IM,γ1,γ2,γ3,γ4,γ5)∈R+MT+M(T-1)+MST+2MNT+KMT+MT×R4MT×R+MT}.In (12), γ1m(t), γ2m(t), γ3m(t), γ4m(t), and γ5m(t) are the Lagrange multipliers corresponding to constraints (6), (7), (8), (9), and (10), and γ1, γ2, γ3, γ4, and γ5 are the MT-dimension column vector with the elements of γ1m(t), γ2m(t), γ3m(t), γ4m(t), and γ5m(t), respectively.From the forth term of (12), we can find that when the transaction amount qmn*(t)>0, ρmn*(t)=(∂cmn*(t)/∂qmn(t))+γ1m*(t)+γ4m*(t).In the model, we consider the manufacturers’ carbon emission constraints from the government in each period, and then the manufacturers’ behaviors can be characterized by the following equation:(13)Bm(t)-(erqmr(t)+euqmu(t)){≥0,ssm(t)=0,=0,ssm(t)>0.Equation (13) describes that if the manufactures’ carbon emission amount is lower than Bm(t) in period t, the carbon tax is zero; on the other hand, if the manufactures’ carbon emission amount is equal to Bm(t) in period t, the manufacturers must pay taxes.The above equilibrium condition corresponds to the following variational inequality: determiness*∈R+MT, satisfying
(14)∑t=1T∑m=1M[Bm*(t)-(erqmr*(t)+euqmu*(t))]×[ssm(t)-ssm*(t)]≥0,∀ss∈R+MT.In a word, the manufacturers’ optimal behaviors must meet (10) and (14) simultaneously.


## 4.3. The Optimal Behavior of Retailers

The retailers need to decide how many products to purchase from every manufacturer and to sell to the consumers in demand markets to meet the random demand in markets.For retailern, the profits can be defined as follows:
(15)πn(Q2,Q5)=max⁡(∑t=1Tρnk(t)qnk(t)hhhhhhhhh-∑t=1Tcn(t)-∑t=1T∑m=1Mρmn(t)qmn(t)hhhhhhhhh-∑t=1TcnkN(t))(16)s.t.qnk(t)≤∑m=1Mqmn(t)(17)qmn(t),qnk(t)≥0,∀m,n,k.We also assume that the retailers compete in a noncooperative manner. The optimal behavior of each retailer is to maximize his profit when the other retailers’ behaviors are given. So, the optimal conditions of all retailers can be described by variational inequality: determine(Q2*,Q5*,μ*)∈ΩN, satisfying
(18)∑t=1T∑n=1N[cnkN*(t)∂qnk(t)-ρnk*(t)+μn*(t)]×[qnk(t)-qnk*(t)]+∑t=1T∑n=1N∑m=1M[ρmn*(t)+∂cn*(t)∂qmn-μn*(t)]hhhhhhhhhhh×[qmn(t)-qmn*(t)]+∑t=1T∑n=1N[∑m=1Mqmn(t)-qnk(t)]×[μn(t)-μn*(t)]≥0hhhhhhhhhhhhhhhhhhhhhhhhhhh∀(Q2,Q5,μ)∈ΩN,
where ΩN={(Q2,Q5,μ)∈R+MNT+NKT+NT}.In (18), μn(t) is the Lagrange multiplier corresponding to constraints (16) and μ is the NT-dimension column vector with the elements of μn(t).From the first term of (18), we find that when the transaction volume qnk(t)>0, ρnk*(t)=(cnkN*(t)/∂qnk(t))+μn*(t).


## 4.4. The Optimal Behavior of Demand Markets

In the forward supply chain, the consumers of each demand market decide to purchase the products from the retailer or not based on the appreciation level of products and also decide the purchasing price and amount.According to literature [19], the transaction between the demand market and the retailer should satisfy a complement condition; that is, for some demand market k,
(19)dk(ρ4(t)){=∑n=1Nqnk(t),ρ4k(t)>0,≤∑n=1Nqnk(t),ρ4k(t)=0.Equation (19) illustrates that if the product amount supplied by each retailer is equal to that of demand in markets, then the transaction can proceed and the transaction price is positive. The demand amount depends on the transaction price paid by consumers in demand market k.And also similar to [19], in the reverse supply chain, the demand markets’ behavior should be characterized by the following equation:
(20)αk(Q4*(t)){=ρkmu*(t),qkmu(t)>0,≥ρkmu*(t),qkmu(t)=0,
where αn(Q4(t)) is the disutility of consumers in the demand market n at period t. Equation (20) expresses that the higher collection price paid by manufacturers, the more used products the manufacturers collect. Moreover, the collected amount is lower than or equal to that of the products sold to the demand market, so the following constraint should be satisfied:
(21)∑m=1Mqkmu(t)≤∑n=1Nqnk(t).Letλk(t) denote the Lagrange multiplier corresponding to (21), group all of λk(t) into a KT-dimension column vector λ.Besides the purchasing price paying to the retailers, the consumers must pay the transaction cost to obtain the products; hence, we have(22)ρnk*(t)+cnkK*(t){=ρ4k(t),qnk*(t)>0,≥ρ4k(t),qnk*(t)=0.Equations (19), (20), (21), and (22) can be described by the following variational inequality: determine (Q4*,Q5*,ρ4*,λ*)∈ΩK, satisfying
(23)∑t=1T∑n=1N∑m=1M[αn(Q4*(t))-ρkmu*(t)+λk*(t)]hhhhhh×[qkmu(t)-qkmu*(t)]+∑t=1T∑n=1N∑k=1K[ρnk*(t)+cnkK*(t)-λk*(t)]hhhhhhhhh×[qnk(t)-qnk*(t)]+∑t=1T∑k=1K[∑n=1Nqnk*(t)-dk(ρ4*(t))]hhhhhhh×[ρ4k(t)-ρ4k*(t)]+∑t=1T∑k=1K[∑n=1Nqnk*(t)-∑m=1Mqkmu*(t)]hhhhhhh×[λk(t)-λk*(t)]≥0,hhhhhhhh∀(Q4,Q5,ρ,λ)∈ΩK,
where ΩK={(Q4,Q5,ρ4,λ)∈R+NMT+NKT+2NT}.From the first term of inequality (23), we can conclude that when the transaction amount of collected products qnmu(t)≥0, the transaction price of collected products between manufacturers and demand markets ρnmu*(t)=αn(Q4*(t))+λn*(t). In fact, since ∑m=1Mqmn(t)>∑m=1Mqnmu(t) and λn*(t)=0, we can get ρnmu*(t)=αn(Q4*(t)).


## 4.5. The Equilibrium Conditions of Multiperiod CLSC Network

Theorem 1.
In equilibrium of the multiperiod CLSC network, the suppliers’ optimal behaviors, as expressed by variational inequality (4), the manufacturers’ optimal behaviors, as expressed by variational inequality (12) and (14), the retailers’ optimal behaviors, as expressed by variational inequality (18), and the consumers’ optimal behaviors, as expressed by variational inequality (23), must be satisfied simultaneously. One can state it in the following variational inequality: determine (q1r*,qr*,qu*,Q1*,Q2*,Q3*,Q4*,IM*,ss*,ρ*,ηs*,γ1*,γ2*,γ3*,γ4*,γ5*,μ*,λ*)∈Ω, satisfying
(24)∑t=1T∑s=1S∑m=1M[∂csmr*(t)∂qsmr(t)+ηs*(t)-γ5m*(t)]hhhhhh×[qsmr(t)-qsmr*(t)]+∑t=1T∑s=1S[∂fsr*(qsr(t))∂qsr(t)-ηs*(t)]hhhhhhh×[qsr(t)-qsr*(t)]+∑t=1T∑m=1M[∂fmr*(t)∂qmr(t)-βrγ1m*(t)+ersm*(t)+γ5m*(t)]hhhhhhhh×[qmr(t)-qmr*(t)]+∑t=2T∑m=1M[∂cm(qmu*(t))∂qmu(t)+γ2m*(t-1)+eusm*(t)]hhhhhhhh×[qmu(t)-qmu*(t)]+∑t=1T∑m=1M∑n=1N[∂cmn*(t)qmn(t)+∂cn*(t)∂qmnhhhhhhhhihhhhh+γ1m*(t)+γ4m*(t)-λk*(t)-μn*(t)∂cmn*(t)qmn(t)+∂cn*(t)∂qmn]hhhhhhhhihh×[qmn(t)-qmn*(t)]+∑t=1T∑k=1K∑m=1M[∂ckmu*(t)∂qkmu(t)-γ2m*(t)-γ3m*(t)hhhhhhhhhhhhh+αn(Q4*(t))+λk*(t)∂ckmu*(t)∂qkmu(t)]hhhhhhhhihh×[qkmu(t)-qkmu*(t)]+∑t=1T∑n=1N∑m=1M[∑i=1min⁡⁡{t,N-1}[∂fmu*(t+i-1)∂qmnr(t)+cmd(t-i)∂qmnr(t)]hhhhhhhhhhhhhhh-∑i=0min⁡⁡{t+1,Z-1}βiαiγ4m*(t+i)hhhhhhhhhhhhhhh-∑i=1min⁡⁡{t+1,Z-1}βiαiγ1m*(t+i)hhhhhhhhhhhhhhh+ρ-∑i=1min⁡⁡{t+1,Z-1}βi-1(1-βiu)αihhhhhhhhhhhhhhh+∑i=1min⁡⁡{t+1,Z-1}βi-1αiγ3m*(t+i-1)]hhhhhhhhihhh×[qmnr(t)-qmnr*(t)]+∑t=1T∑m=1M[∂Hm*(t)∂Im(t)+γ1m*(t)-γ1m*(t+1)]hhhhhhhh×[Im(t)-Im*(t)]+∑t=1T∑m=1M[Bm(t)-(erqmr*(t)+euqmu*(t))]hhhhhhhh×[ssm(t)-ssm*(t)]+∑t=1T∑n=1N[cnkN*(t)∂qnk(t)+cnkK*(t)-λk*(t)+μn*(t)]hhhhh×[qnk(t)-qnk*(t)]+∑t=1T∑k=1K[∑n=1Nqnk*(t)-dk(ρ4*(t))]hhhhhhhh×[ρ4k(t)-ρ4k*(t)]+∑t=1T∑s=1S[qsr*(t)-∑m=1Mqsmr*(t)]hhhhhhhh×[ηs(t)-ηs*(t)]+∑t=1T∑m=1M[∑n=1N∑i=1min⁡⁡{t-1,Z-1}βiαiqmnr*(t-i)+βrqmr*(t)hhhhhhhhhhh+Im*(t-1)-∑n=1Nqmn*(t)-Im*(t)]hhhhhhhh×[γ1m(t)-γ1m*(t)]+∑t=1T∑m=1M[∑n=1Nqnmu*(t)-qmu*(t+1)]hhhhhhhh×[γ2m(t)-γ2m*(t)]+∑t=2T∑m=1M[∑n=1Nqnmu*(t-1)hhhhhhhhhh-∑n=1N∑i=1min⁡⁡{t-1,Z-1}βi-1αiqmnr*(t-i)]hhhhhhhh×[γ3m(t)-γ3m*(t)]+∑t=1T∑m=1M[∑n=1Nqmnr*(t)hhhhhhhhhhh+∑n=1N∑i=1min⁡⁡{t-1,Z-1}βiαiqmnr*(t-i)-∑n=1Nqmn*(t)]hhhhhhhh×[γ4m(t)-γ4m*(t)]+∑t=1T∑m=1M[∑s=1Sqsmr*(t)-qmr*(t)]hhhhhhhh×[γ5m(t)-γ5m*(t)]+∑t=1T∑n=1N[∑m=1Mqmn(t)-qnk(t)]hhhhhhh×[μn(t)-μn*(t)]+∑t=1T∑k=1K[qnk*(t)-∑m=1Mqkmu*(t)]hhhhhhh×[λk(t)-λk*(t)]≥0,∀(q1r,qr,qu,Q1,Q2,Q3,Q4,IM,ss,hhv,ρ,ηs,γ1,γ2,γ3,γ4,γ5,μ,λ)∈Ω,
where Ω=ΩS×ΩM×ΩN×ΩK.


## 4.6. Solving Algorithm and Parameter Setting

For easy formulation hereinafter, we group the terms of the multiplication signs in inequality (24) into a column vector vector F(y)={Fsmt1,Fst2,Fmt3,Fm(t-1)4,Fmnt5,Fmnt6,Fkmt7,Fmt8,Fmt9,Fnt10,Fnt11,Fst12,Fmt13,Fmt14,Fmt15,Fmt16,Fmt17,Fnt18,Fnt19}∀s,m,n,k,t and introduce y=(q1r,qr,qu,Q1,Q2,Q3,Q4,IM,ss,v,ρ,ηs,γ1,γ2,γ3,γ4,γ5,μ,λ)∈Ω; thus, we can rewrite the variational inequality problem (24) in the standard form: determine y*=(q1r*,qr*,qu*,Q1*,Q2*,Q3*,Q4*,IM*,ss*,v*,ρ*,ηs*,γ1*,γ2*,γ3*,γ4*,γ5*,μ*,λ*), satisfying
(25)F(y)T·(y-y*)≥0,∀y∈Ω.In order to solve the variational inequality (25), we have several algorithms to select, for example, the modified project contraction algorithm [31], logarithmic-quadratic proximal prediction-correction method [32]. In this paper, we employ the modified project contraction algorithm to solve the problem for two reasons. Firstly, the modified project contraction algorithm provides an iterative solution framework that can deal with variational inequalities defined on polyhedral sets and get the optimal solutions and Lagrangian multipliers simultaneously. So, it satisfies all the requirements of a solution method for our model. Secondly, the design of its iterative steps is much simpler than that of logarithmic-quadratic proximal prediction-correction method. Set the related parameters of the modified project contraction algorithm as follows: the original value of decision variables and Lagrange multipliers is set to 1; the convergence criterion, for example, the absolute value of difference of decision variables and Lagrange multipliers between two steps, is lower than or equal to 10−8.


## 5. Numerical Examples

Now, we consider a CLSC network with two suppliers, two manufacturers, two retailers, and two demand markets, the lifetime of products is two or three, and the planning horizon is five. In this construction, we give some numerical examples to illustrate the efficiency of our model and analyze the influence of some parameters to obtain managerial insights. The related cost functions are given as follows.Suppliers’ raw material producing cost functions:(26)fsr(q1r(t))=(1+4t)qsr(t)2+qsr(t)+1,hhhhhhhhhhis=1,2;t=1,2,3,4,5;manufacturers’ production cost functions using raw material:(27)fmr(t)=(0.5+1.5t)(βrqmr(t))2+βrqmr(t)βrq3-mr(t)+2βrqmr(t),hhhhhhm=1,2;t=1,2,3,4,5;manufacturers’ remanufacturing cost functions using collected products:(28)fmu(2)=4.5(β1uα∑n=12qmnr(1))2+1.5(β1uα∑n=12qmnr(1)),hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhm=1,2;fmu(t)=(4.5-0.5t)×(β1uα∑n=12qmnr(t-1)+2δβ2uβ1uα2∑n=12qmnr(t-2))2+1.5(β1uα∑n=12qmnr(t-1)+2δβ2uβ1uα2∑n=12qmnr(t-2)),hhhhhhhhhhhhhhhhhhhhhhhihhhm=1,2;t=3,4,5;the manufacturers’ inventory cost functions:(29)Hm(t)=(2+(-1)t0.5)Im(t),hhhhm=1,2;t=1,2,3,4,5;the manufacturers’ disposal cost functions dealing with the collected products:(30)cm(qmu(t))=2.5qmu(t)2+qmu(t)+2,hhhhhhhhiim=1,2;t=2,3,4,5;the transportation cost functions from demand markets to manufacturers:(31)cnmu(t)=0.5qnmu(t)2+3qnmu(t)+1,m=1,2;n=1,2;t=2,3,4,5;the manufacturers’ transportation cost functions to move the waste of disposal products to landfill:(32)cmd(2)=0.5[(1-β1u)α∑n=12qmnr(1)]2+3.5(1-β1u)α∑n=12qmnr(1),cmd(t)=0.5[δ(1-β2u)β1uα2∑n=12qmnr(t-2)+(1-β1u)α∑n=12qmnr(t-1)]2+3.5[δ(1-β2u)β1uα2∑n=12qmnr(t-2)+(1-β1u)αhhhhh×∑n=12qmnr(t-1)]fort=3,4,5;m=1,2;the consumers’ disutility functions in demand markets:(33)αn(Q4(t))=0.5∑m=12∑n=12qnmu(t),hhforn=1,2;t=1,2,3,4,5;the transaction cost functions undertaken by the supplier between suppliers and manufacturers:(34)csmr(t)=0.5(qsmr(t))2+qsmr(t),fors=1,2;m=1,2;t=1,2,3,4,5;the transaction cost functions undertaken by the manufacturer between manufacturers and retailers:(35)cmn(t)=(qmn(t))2+3qmn(t)+0.5,form=1,2;n=1,2;t=1,2,3,4,5;the disposal functions of retailers:(36)cn(Q2)=0.25(∑m=1Mqmn(t))2,forn=1,2;t=1,2,3,4,5;the transaction cost functions undertaken by retailern to demand market k:
(37)cnkN(t)=0.5(qnk(t))2+qnk(t)+1,n=1,2;k=1,2;t=1,2,3,4,5;the transaction cost functions undertaken by demand markets purchasing products from retailers:(38)cnkK(t)=qnk(t)+0.5,n=1,2;t=1,2,3,4,5;the demand functions of the consumer markets:(39)dk(ρtD)=-(2.2-0.4t)ρ3-kD(t)-2ρkD(t)+120,hhhhhhhhhhhhk=1,2,t=1,2,3,4,5.In order to illustrate the impact of lifetime on equilibrium results, we introduceδ, a Boolean parameter, to represent two cases of lifetime values: δ=0 denotes that the product lifetime is two, and δ=1 denotes that the product lifetime is three. The total planning horizon is five.


### 5.1. Numerical Example 1

Whenδ=0, for example, the product lifetime is 2, ρ-=1, βr=1, β1u=0.9, and β2u=0.6, we analyze the effect of the collection rate α on the equilibrium results given in Table 1, in which the data in each column are the equilibrium solutions corresponding to a value of α and the data in each row are the optimal values of main variables.Table 1
δ
=
0, the CLSC network equilibrium results when α increases.
α
0.05
0.1
0.15
0.2
0.25
0.3
q
s
r
,
q
m
r  s=1,2;m=1,2
3.09271.97561.55401.27531.1344
3.19292.03571.59701.30871.1617
3.28622.09171.63691.33981.1871
3.37232.14341.67381.36851.2106
3.45102.19061.70761.39481.2321
3.52192.23321.73801.41841.2514
q
m
u  m=1,2
0.01620.04600.07370.1339
0.04630.10600.15950.2871
0.09010.17810.25370.4638
0.14690.26120.35170.6695
0.21610.35520.44800.9106
0.29650.46140.53581.1957
q
m
n
r  m=1,2;n=1,2
0.16200.45990.73741.33931.8174
0.23170.52990.79771.43531.6534
0.30020.59370.84551.54601.4854
0.36720.65300.87921.67371.3112
0.43220.71040.89611.82121.1281
0.49410.76900.89291.99280.9325
q
m
n  m=1,2;n=1,2
0.16200.46720.75811.37251.8776
0.23170.55070.84541.50711.7826
0.30020.63420.92571.66021.6942
0.36720.71910.99671.83191.6125
0.43220.80761.05592.02291.5379
0.49410.90241.10062.23391.4706
I
m  m=1,2
2.76873.82453.90372.50030
2.72953.70553.70702.14510
2.68583.59013.53591.78370
2.63783.47513.39061.41180
2.58663.35643.27181.02410
2.53373.22883.18090.61360
ρ
n  n=1,2
30.669733.585037.478341.853048.5833
30.493633.368737.245641.480348.8597
30.320533.152637.031541.056549.1170
30.151232.932836.842140.580949.3546
29.987032.703736.684240.052149.5717
29.830732.458436.565139.467549.7674
ρ
k
m
u  k=1,2
0.01620.04600.07370.13390
0.04630.10600.15950.28710
0.09010.17810.25370.46380
0.14690.26120.35170.66950
0.21610.35520.44800.91060
0.29650.46140.53581.19570
π
s
,
s
=
1,2
32.2546
34.7208
37.0795
39.3119
41.3978
43.3149
π
m
,
m
=
1,2
75.9138
77.0719
77.1677
75.8899
72.8294
67.4213
π
n
,
n
=
1,2
18.6856
19.5621
20.9267
22.8044
25.2330
28.2688
Total profit
126.8549
131.3548
135.1739
138.0062
139.4602
139.0050From Table1, we can find that when α increases, the volume of raw material qsr (equal to qmr) supplied to manufacturers, the used products qmu collected by manufacturers, and the transaction qmn between manufacturers and retailers increase too, so the increase of α stimulates the transaction of the network. Moreover, in order to obtain and remanufacture more used products, the manufacturers and retailers should sell more products and hence collect more used products at the end of each period, so that the transaction price ρn between retailers and demand markets is decreasing all the time except the last period. This point characterizes the dynamic network fully in which the players are foresighted rather than shortsighted.From Table1, we also find that Im(t) decreases when α increases. It is because of the transaction volume increasing, and the manufacturers have no motivation to transfer more inventories to next period, and Im(5)=0 means all the surplus products are sold out in the last period.In the following, we analyze the impact ofα increasing on the profits of players in the CLSC network. The suppliers’ and retailers’ profits increase because of their transaction volume with manufacturers increasing. The profits of manufacturers increase at the beginning but then decrease, and the maximum values appear at about the point 0.15 of α, so we must explain this interesting phenomena. When α increases, the amount of collected products also increases; therefore, the collection price increases; that is, the manufacturers must pay more to demand markets for obtaining extra used products. When the collection price reaches a threshold, the cost saving via remanufacturing cannot make up the increased collection cost, so the manufacturers’ profits reach the maximum and then decrease when the price continues to increase.


### 5.2. Numerical Example 2

Whenδ=1, for example, the product lifetime is 3, ρ-=1, βr=1, β1u=0.9, and β2u=0.6, we analyze the impact of the collection rate α on the equilibrium results. The equilibrium results are shown in Table 2.Table 2
δ
=
1, the CLSC network equilibrium states when α increases.
α
0.05
0.1
0.15
0.2
0.25
0.3
q
s
r
,
q
m
r s=1,2; m=1,2
3.10731.98441.56031.28021.1384
3.22812.05681.61201.32051.1713
3.34702.12821.66301.36011.2037
3.46272.19761.71261.39871.2353
3.57342.26411.76001.43561.2655
3.67722.32631.80451.47021.2938
q
m
u m=1,2
0.01770.04830.07780.1358
0.05300.11750.18020.2946
0.10700.20870.31050.4802
0.18050.32220.47210.6971
0.27380.45720.66770.9499
0.38700.61160.89961.2440
q
m
n
r m=1,2;n=1,2
0.17660.47480.75641.32431.8031
0.26520.56370.85001.39651.6190
0.35680.64750.94761.47291.4261
0.45120.72421.04981.55371.2244
0.54760.79121.15731.63931.0139
0.64500.84521.27111.73020.7944
q
m
n m=1,2;n=1,2
0.17660.48280.77801.35901.8637
0.26520.58760.90221.47601.7492
0.35680.69571.03941.60871.6365
0.45120.80541.18991.75831.5268
0.54760.91441.35381.92641.4218
0.64501.01941.53062.11451.3234
I
m m=1,2
2.75413.78883.83632.46790
2.69773.62723.53922.06660
2.63333.46653.23431.64860
2.56033.30942.92241.21360
2.47833.15982.60530.76220
2.38713.02302.28530.29510
π
s
,
s
=
1,2
32.6103
35.6032
38.6502
41.7106
44.7290
47.6346
π
m
,
m
=
1,2
75.7360
76.3586
75.2496
71.6741
64.7021
53.2055
π
n
,
n
=
1,2
18.5700
19.4041
20.8727
23.0729
26.1039
30.0609
Total profit
126.9163
131.3659
134.7725
136.4576
135.5350
130.9010In Table2, the changing trends of some variables such as qsr, qmr, qmu, qmn, Im, and profits of suppliers, manufacturers, and retailers are similar to those in Table 1 depending on α, and, in turn, we will pay more attention to the differences between the two tables.In respect to the profits of suppliers and retailers, there are changes but not large in the course ofα increasing. But the profits of manufacturers change dramatically. The maximum value of manufacturers’ profit is 77.1677 in Table 1, which is at about the point 0.1 of α and larger than the corresponding value 76.3586 in Table 2. This phenomenon shows that the smaller the product lifetime is, the more profits the manufacturers can earn, which means that the environmental and resource targets are different from those of the enterprise to seek maximized profits.For CLSC, the maximum total profits are 139.4602 and 136.4576 at the points 0.25 and 0.2 ofα when δ=0 and δ=1, respectively. The maximum profits of the whole CLSC and the manufacturers appear in different points of α, so the objects of the manufactures and the whole CLSC are also different. If the manufacturers can give up a part of interests, then the performance of CLSC will be better.Moreover, we compare the CLSC profits in two cases ofδ=0 and δ=1 when α is set to be a same value. It can be seen that when α≤0.1, the CLSC profit with Z=3 is larger than that with Z=2. However, when α≥0.15, the quantitative relationship between the two is opposite. The phenomenon shows clearly that more product lifetime is not always a benefit to the CLSC. Too high collection rate will result in the fact that the collection burden of the manufacturers aggravates and the operational efficiency of the CLSC reduces, although it contributes to the cost saving by remanufacturing. Therefore, the government should set the collection rate in a suitable region based on the manufacturers’ collection cost structure and benefits from remanufacturing.


### 5.3. Numerical Example 3

WhenBm(1)=3, δ=0, ρ-=1, βr=1, β1u=0.9, and β2u=0.6, we analyze the impact of the collection rate α on the equilibrium results. The profits of various players in CLSC are shown in Figure 2.Profits of various players in CLSC.
(a)
Profit of suppliers
(b)
Profit of manufacturers
(c)
Profit of retailersFrom Figure2, we can make a conclusion that as α increases, the profits of suppliers and retailers increase while the profit of manufacturers increases at the beginning and then decreases. This is the same as the above analysis in Table 1. The maximum value of the profit of manufacturer differs at another point of α; that is, the maximum value point in Table 1 is at α=0.15 but at α=0.1 in Figure 2. This is because of the existence of carbon emission cap in the 1st period; the profits of all players are lower than that without carbon emission cap by comparing Figure 2 with Table 1.


### 5.4. Numerical Example 4

Whenα=0.25 and δ=0, ρ-=1, βr=1, β1u=0.9, and β2u=0.6, we analyze the effect of different carbon emission caps of the first 3 periods on the equilibrium results shown in Table 3.Table 3
δ
=
0 and α=0.25, the CLSC network equilibrium results of various emission caps.
B
m
(
1
)
,
m
=
1,2
3.2
3
3.2
3.2
3.2
3.2
B
m
(
2
)
,
m
=
1,2
∞
∞
2.3
2
2.3
2.3
B
m
(
3
)
,
m
=
1,2
∞
∞
∞
∞
1.8
1.6
q
s
r
,
q
m
r s=1,2;m=1,2
3.22.20791.72001.40441.2400
32.22181.72981.41211.2463
3.22.19751.72071.40501.2405
3.21.89061.74301.42231.2546
3.22.20251.62211.41081.2452
3.22.21311.41081.42341.2555
q
m
u m=1,2
0.20450.34620.43830.9018
0.19520.33900.43060.8947
0.20500.34540.43780.9012
0.21870.32140.42210.8849
0.19510.35590.42820.8967
0.17390.37840.40780.8870
s
m m=1,2
2.099300
3.772100
2.11260.14360
2.50254.36940
2.24410.21321.8580
2.52560.36225.8363
π
s
,
s
=
1,2
39.8379
38.7104
39.7563
37.5784
38.8658
37.1742
π
n
,
n
=
1,2
70.3260
68.5201
70.1648
65.7157
68.0905
63.8557
π
m
,
m
=
1,2
24.3986
23.7460
24.3539
23.0658
23.9874
23.2293
Total profit
134.5625
130.9765
134.2750
126.3599
130.9437
124.2592Table3 presents main variables associated with different carbon caps. As mentioned above, sm are the Lagrange multiplier associated with the carbon cap, and sm>0 means that the carbon emission reaches the cap. We compare the first 2 columns in Table 3 and note that when the carbon emission cap of the 1st period becomes smaller, Lagrange multiplier associated with the 1st period is larger, while the amounts of raw material supplied to the manufacturers are larger except in the 1st period due to the carbon emission constraints. The amount of collected products of the manufacturers is smaller in all planning periods.The profits of all players in the network decrease when the cap decreases. So, the carbon emission constraint is harmful to the CLSC in the economic aspect without considering the environmental target.Now, we compare the 1st column with the 3rd and 4th columns, in which the carbon emission cap in the 1st period is fixed, and in the 2nd period it appears and then becomes smaller. The amount of raw materials supplied to the manufacturers is larger than the 3rd period because there are no constraints since this period. The amount of collected products by the manufacturers is larger in the 1st period and smaller in the other periods. The profits of all players in the network decrease when the cap of the 2nd period decreases.Finally, we compare the 3rd column with the 5th and 6th columns. The amount of raw materials supplied to the manufacturers is larger than the 4th period with no constraints since this period. The amount of collected products by the manufacturers is larger in the 2nd period and smaller in other periods. The profits of all players in the network also decrease.From the analysis above, we can make the following conclusion: when setting carbon emission caps in some periods, the supply amounts of raw materials in these periods become lower, while the amounts of collected products in the previous period are larger with the caps; in other period, the collected used products of the manufacturers are smaller. The profits of all players in the CLSC network decrease as the caps become smaller. With the carbon emission caps appearing in Table3, the decision-makers consider the whole planning horizon comprehensively to obtain the optimal strategy, but not merely a single period.


## 5.1. Numerical Example 1

Whenδ=0, for example, the product lifetime is 2, ρ-=1, βr=1, β1u=0.9, and β2u=0.6, we analyze the effect of the collection rate α on the equilibrium results given in Table 1, in which the data in each column are the equilibrium solutions corresponding to a value of α and the data in each row are the optimal values of main variables.Table 1
δ
=
0, the CLSC network equilibrium results when α increases.
α
0.05
0.1
0.15
0.2
0.25
0.3
q
s
r
,
q
m
r  s=1,2;m=1,2
3.09271.97561.55401.27531.1344
3.19292.03571.59701.30871.1617
3.28622.09171.63691.33981.1871
3.37232.14341.67381.36851.2106
3.45102.19061.70761.39481.2321
3.52192.23321.73801.41841.2514
q
m
u  m=1,2
0.01620.04600.07370.1339
0.04630.10600.15950.2871
0.09010.17810.25370.4638
0.14690.26120.35170.6695
0.21610.35520.44800.9106
0.29650.46140.53581.1957
q
m
n
r  m=1,2;n=1,2
0.16200.45990.73741.33931.8174
0.23170.52990.79771.43531.6534
0.30020.59370.84551.54601.4854
0.36720.65300.87921.67371.3112
0.43220.71040.89611.82121.1281
0.49410.76900.89291.99280.9325
q
m
n  m=1,2;n=1,2
0.16200.46720.75811.37251.8776
0.23170.55070.84541.50711.7826
0.30020.63420.92571.66021.6942
0.36720.71910.99671.83191.6125
0.43220.80761.05592.02291.5379
0.49410.90241.10062.23391.4706
I
m  m=1,2
2.76873.82453.90372.50030
2.72953.70553.70702.14510
2.68583.59013.53591.78370
2.63783.47513.39061.41180
2.58663.35643.27181.02410
2.53373.22883.18090.61360
ρ
n  n=1,2
30.669733.585037.478341.853048.5833
30.493633.368737.245641.480348.8597
30.320533.152637.031541.056549.1170
30.151232.932836.842140.580949.3546
29.987032.703736.684240.052149.5717
29.830732.458436.565139.467549.7674
ρ
k
m
u  k=1,2
0.01620.04600.07370.13390
0.04630.10600.15950.28710
0.09010.17810.25370.46380
0.14690.26120.35170.66950
0.21610.35520.44800.91060
0.29650.46140.53581.19570
π
s
,
s
=
1,2
32.2546
34.7208
37.0795
39.3119
41.3978
43.3149
π
m
,
m
=
1,2
75.9138
77.0719
77.1677
75.8899
72.8294
67.4213
π
n
,
n
=
1,2
18.6856
19.5621
20.9267
22.8044
25.2330
28.2688
Total profit
126.8549
131.3548
135.1739
138.0062
139.4602
139.0050From Table1, we can find that when α increases, the volume of raw material qsr (equal to qmr) supplied to manufacturers, the used products qmu collected by manufacturers, and the transaction qmn between manufacturers and retailers increase too, so the increase of α stimulates the transaction of the network. Moreover, in order to obtain and remanufacture more used products, the manufacturers and retailers should sell more products and hence collect more used products at the end of each period, so that the transaction price ρn between retailers and demand markets is decreasing all the time except the last period. This point characterizes the dynamic network fully in which the players are foresighted rather than shortsighted.From Table1, we also find that Im(t) decreases when α increases. It is because of the transaction volume increasing, and the manufacturers have no motivation to transfer more inventories to next period, and Im(5)=0 means all the surplus products are sold out in the last period.In the following, we analyze the impact ofα increasing on the profits of players in the CLSC network. The suppliers’ and retailers’ profits increase because of their transaction volume with manufacturers increasing. The profits of manufacturers increase at the beginning but then decrease, and the maximum values appear at about the point 0.15 of α, so we must explain this interesting phenomena. When α increases, the amount of collected products also increases; therefore, the collection price increases; that is, the manufacturers must pay more to demand markets for obtaining extra used products. When the collection price reaches a threshold, the cost saving via remanufacturing cannot make up the increased collection cost, so the manufacturers’ profits reach the maximum and then decrease when the price continues to increase.


## 5.2. Numerical Example 2

Whenδ=1, for example, the product lifetime is 3, ρ-=1, βr=1, β1u=0.9, and β2u=0.6, we analyze the impact of the collection rate α on the equilibrium results. The equilibrium results are shown in Table 2.Table 2
δ
=
1, the CLSC network equilibrium states when α increases.
α
0.05
0.1
0.15
0.2
0.25
0.3
q
s
r
,
q
m
r s=1,2; m=1,2
3.10731.98441.56031.28021.1384
3.22812.05681.61201.32051.1713
3.34702.12821.66301.36011.2037
3.46272.19761.71261.39871.2353
3.57342.26411.76001.43561.2655
3.67722.32631.80451.47021.2938
q
m
u m=1,2
0.01770.04830.07780.1358
0.05300.11750.18020.2946
0.10700.20870.31050.4802
0.18050.32220.47210.6971
0.27380.45720.66770.9499
0.38700.61160.89961.2440
q
m
n
r m=1,2;n=1,2
0.17660.47480.75641.32431.8031
0.26520.56370.85001.39651.6190
0.35680.64750.94761.47291.4261
0.45120.72421.04981.55371.2244
0.54760.79121.15731.63931.0139
0.64500.84521.27111.73020.7944
q
m
n m=1,2;n=1,2
0.17660.48280.77801.35901.8637
0.26520.58760.90221.47601.7492
0.35680.69571.03941.60871.6365
0.45120.80541.18991.75831.5268
0.54760.91441.35381.92641.4218
0.64501.01941.53062.11451.3234
I
m m=1,2
2.75413.78883.83632.46790
2.69773.62723.53922.06660
2.63333.46653.23431.64860
2.56033.30942.92241.21360
2.47833.15982.60530.76220
2.38713.02302.28530.29510
π
s
,
s
=
1,2
32.6103
35.6032
38.6502
41.7106
44.7290
47.6346
π
m
,
m
=
1,2
75.7360
76.3586
75.2496
71.6741
64.7021
53.2055
π
n
,
n
=
1,2
18.5700
19.4041
20.8727
23.0729
26.1039
30.0609
Total profit
126.9163
131.3659
134.7725
136.4576
135.5350
130.9010In Table2, the changing trends of some variables such as qsr, qmr, qmu, qmn, Im, and profits of suppliers, manufacturers, and retailers are similar to those in Table 1 depending on α, and, in turn, we will pay more attention to the differences between the two tables.In respect to the profits of suppliers and retailers, there are changes but not large in the course ofα increasing. But the profits of manufacturers change dramatically. The maximum value of manufacturers’ profit is 77.1677 in Table 1, which is at about the point 0.1 of α and larger than the corresponding value 76.3586 in Table 2. This phenomenon shows that the smaller the product lifetime is, the more profits the manufacturers can earn, which means that the environmental and resource targets are different from those of the enterprise to seek maximized profits.For CLSC, the maximum total profits are 139.4602 and 136.4576 at the points 0.25 and 0.2 ofα when δ=0 and δ=1, respectively. The maximum profits of the whole CLSC and the manufacturers appear in different points of α, so the objects of the manufactures and the whole CLSC are also different. If the manufacturers can give up a part of interests, then the performance of CLSC will be better.Moreover, we compare the CLSC profits in two cases ofδ=0 and δ=1 when α is set to be a same value. It can be seen that when α≤0.1, the CLSC profit with Z=3 is larger than that with Z=2. However, when α≥0.15, the quantitative relationship between the two is opposite. The phenomenon shows clearly that more product lifetime is not always a benefit to the CLSC. Too high collection rate will result in the fact that the collection burden of the manufacturers aggravates and the operational efficiency of the CLSC reduces, although it contributes to the cost saving by remanufacturing. Therefore, the government should set the collection rate in a suitable region based on the manufacturers’ collection cost structure and benefits from remanufacturing.


## 5.3. Numerical Example 3

WhenBm(1)=3, δ=0, ρ-=1, βr=1, β1u=0.9, and β2u=0.6, we analyze the impact of the collection rate α on the equilibrium results. The profits of various players in CLSC are shown in Figure 2.Profits of various players in CLSC.
(a)
Profit of suppliers
(b)
Profit of manufacturers
(c)
Profit of retailersFrom Figure2, we can make a conclusion that as α increases, the profits of suppliers and retailers increase while the profit of manufacturers increases at the beginning and then decreases. This is the same as the above analysis in Table 1. The maximum value of the profit of manufacturer differs at another point of α; that is, the maximum value point in Table 1 is at α=0.15 but at α=0.1 in Figure 2. This is because of the existence of carbon emission cap in the 1st period; the profits of all players are lower than that without carbon emission cap by comparing Figure 2 with Table 1.


## 5.4. Numerical Example 4

Whenα=0.25 and δ=0, ρ-=1, βr=1, β1u=0.9, and β2u=0.6, we analyze the effect of different carbon emission caps of the first 3 periods on the equilibrium results shown in Table 3.Table 3
δ
=
0 and α=0.25, the CLSC network equilibrium results of various emission caps.
B
m
(
1
)
,
m
=
1,2
3.2
3
3.2
3.2
3.2
3.2
B
m
(
2
)
,
m
=
1,2
∞
∞
2.3
2
2.3
2.3
B
m
(
3
)
,
m
=
1,2
∞
∞
∞
∞
1.8
1.6
q
s
r
,
q
m
r s=1,2;m=1,2
3.22.20791.72001.40441.2400
32.22181.72981.41211.2463
3.22.19751.72071.40501.2405
3.21.89061.74301.42231.2546
3.22.20251.62211.41081.2452
3.22.21311.41081.42341.2555
q
m
u m=1,2
0.20450.34620.43830.9018
0.19520.33900.43060.8947
0.20500.34540.43780.9012
0.21870.32140.42210.8849
0.19510.35590.42820.8967
0.17390.37840.40780.8870
s
m m=1,2
2.099300
3.772100
2.11260.14360
2.50254.36940
2.24410.21321.8580
2.52560.36225.8363
π
s
,
s
=
1,2
39.8379
38.7104
39.7563
37.5784
38.8658
37.1742
π
n
,
n
=
1,2
70.3260
68.5201
70.1648
65.7157
68.0905
63.8557
π
m
,
m
=
1,2
24.3986
23.7460
24.3539
23.0658
23.9874
23.2293
Total profit
134.5625
130.9765
134.2750
126.3599
130.9437
124.2592Table3 presents main variables associated with different carbon caps. As mentioned above, sm are the Lagrange multiplier associated with the carbon cap, and sm>0 means that the carbon emission reaches the cap. We compare the first 2 columns in Table 3 and note that when the carbon emission cap of the 1st period becomes smaller, Lagrange multiplier associated with the 1st period is larger, while the amounts of raw material supplied to the manufacturers are larger except in the 1st period due to the carbon emission constraints. The amount of collected products of the manufacturers is smaller in all planning periods.The profits of all players in the network decrease when the cap decreases. So, the carbon emission constraint is harmful to the CLSC in the economic aspect without considering the environmental target.Now, we compare the 1st column with the 3rd and 4th columns, in which the carbon emission cap in the 1st period is fixed, and in the 2nd period it appears and then becomes smaller. The amount of raw materials supplied to the manufacturers is larger than the 3rd period because there are no constraints since this period. The amount of collected products by the manufacturers is larger in the 1st period and smaller in the other periods. The profits of all players in the network decrease when the cap of the 2nd period decreases.Finally, we compare the 3rd column with the 5th and 6th columns. The amount of raw materials supplied to the manufacturers is larger than the 4th period with no constraints since this period. The amount of collected products by the manufacturers is larger in the 2nd period and smaller in other periods. The profits of all players in the network also decrease.From the analysis above, we can make the following conclusion: when setting carbon emission caps in some periods, the supply amounts of raw materials in these periods become lower, while the amounts of collected products in the previous period are larger with the caps; in other period, the collected used products of the manufacturers are smaller. The profits of all players in the CLSC network decrease as the caps become smaller. With the carbon emission caps appearing in Table3, the decision-makers consider the whole planning horizon comprehensively to obtain the optimal strategy, but not merely a single period.


## 6. Conclusion

In this paper, we studied the multiperiod CLSC network including manufacturers, retailers, and demand markets taking the product lifetime into account. By using the variational inequalities and complementary theory, we established the CLSC network equilibrium model and solved the model by LQP P-C algorithm. In the end, we provided four numerical examples to illustrate the validity of the modeling framework with or without carbon emission cap, and then we explored the managerial insights in practice. Through the analysis in the numerical examples, we make conclusions as follows.(1)
The increase of product lifetime will enhance the burden of manufacturers and then decrease their profits; the other players in the CLSC are almost not affected; the objects of the manufactures and the whole CLSC are different in general.(2)
The increase ofα stimulates the transaction of the network, but it is not always a benefit to the CLSC. The government should set the collection rate in a suitable region based on the manufacturers’ collection cost structure and benefits from remanufacturing.(3)
The existence of carbon emission cap will lead to the fact that the profits of all players are lower than those with no carbon emission cap.(4)
When setting carbon emission caps in some periods, the supply amounts of raw materials become lower in these periods, while the amount of collected products is larger in the previous period with the caps; the endogenous carbon tax is larger when the cap of carbon emission is lower.We hope that the insights and corresponding results obtained in this paper may be useful and helpful in theory and practice and enhance the operating quality of enterprises in CLSC.



---

*Source: 784637-2014-06-16.xml*