Abstract
Different material supplyrelated decisions intensively affect the efficiency of manufacturers. To obtain a suitable supplyrelated protocol, this study proposes a supply selection model which considers both manufacturers’ development orientation and material ordering. In contrast to traditional approaches that rely on expert opinions, the proposed approach in this study allows the time series analysis (ARIMA) to forecast the trend in manufacturers’ development during the execution of the plan. Based on the predicted trend, taking the minimum of total material management cost as the objective function, the control function and optimization conditions are constructed to select the appropriate protocol. The dynamic prediction protocol is obtained by considering the variation in production and material costs by an evolutionary algorithm. The model enables users to determine material supply protocol in continuous time and autonomously adapt to changes in the manufacturers’ production goals within a lower convergence time.
1. Introduction
The material ordering management solution needs to determine the quantities and ordering paths for product raw materials in a timevarying environment. In some papers, different relative operators were investigated in the last years. Inventory issues and the selection of materialcontrol solutions have been important parts of the overall supply chain coordination process and have been the focus of attention in the field of supply chain management research. The final object of the economical activities is to increase the people’s happiness [1], and the optimal program can save resources then thus contributing to promoting the people’s satisfaction. The optimal program can save resources and thus is to contribute to people’s happiness. And the model can also be used to optimize other domain issues, such as some transportation schedules. In this study, ordering planning is based on a manufacturer’s perspective consideration of several questions from a theoretical point of view, but also from a practical one: What factors influence ordering planning? How will changes within the company orient the ordering planning of raw materials? How should the planning model respond to this orientation? These questions not only suggest how the corresponding model should be built but also determine the situation in which the model needs to be adapted. The study improves the adaptability and generalization of the prediction model and optimizes the information interaction and response speed of the ordering and distribution network.
This paper is outlined as follows. In Sections 2 and 3, (1) we investigate the research background of this problem, (2) we propose dynamic multiobjective optimization based on development prediction, and (3) we propose an evolutionary algorithm combined with dynamic strategy (RS/AS) to solve dynamic multiobjective optimization problem. In Section 4, we investigate two cases to clearly demonstrate the experimental process and to prove the effectiveness of the method. In Section 5, a brief conclusion is provided to prove insights from above researches.
2. Literature Reviews
Recent years have witnessed remarkable progress in material ordering and transportation, with a variety of approaches proposed in some papers and applied in the real application. The method of finding the optimal ordering strategy that minimizes the total inventory by an EPQ model has been extensively studied in the past literature. Wagner and Whitin [2] proposed a forward dynamic algorithm for minimum dynamic inventory for the assessment of the manufacturer’s raw material demand situation and gave the ordering strategy in this case. Mirmohammadi et al. [3] proposed a material demand planning based on a limited period and constructed an optimal sequential tree to solve this problem, but its computational time may increase in exponential form with the sharp growing number of data. The rapid development of the computer and information sciences has witnessed that the big financial data used broadly in the investment field [4]. And actually, the text data by crawling can also be applied to the optimal decision of the supply chain. Qiao and Zhang [5] applied the dynamic multiobjective algorithm to wastewater treatment and showed that the parameters of the Pareto front undergo some changes during the multiobjective optimal control process, but there is no need to restart the optimization process due to the similarity of two consecutive changes. However, the above methods still suffer from limited accuracy and premature convergences. By combining the above literatures, this study establishes an adaptive dynamic multiobjective optimization model for the dynamically changing process of material ordering control to analyze the ordering and storage strategies under different goal orientations and obtain the optimal results using multiple search algorithms which save computational resources and increase rate of convergence by exploiting the similarity of dynamic scenarios and proposes precise suggestions for different goal orientations according to the differences like policy orientations and material properties.
3. Methods
3.1. Future Demand Forecast
According to the need of constructing the model, it is necessary to have a general understanding of the development cycle of the manufacturer’s business. It is more interpretable to humans to quantify the impact of both environment and the manufacture itself. The special events such as COVID19 pandemic have brought great to the price volatility [6], and further, the price volatility will influence the supply chain. Time series forecasting reveals future trends through the historical data of the time series, i.e., the seasonal cycle can be obtained through time series forecasting thus facilitating the development of the length of time for which the ordering strategy applies similarly to the method which Benvenuto et al. used [7]. The outliers in this prediction are clustered and denoised to obtain a data series that is easy to deal with [8]. After data processing and cleaning, a general observation of the data shows that the series contains both a longterm trend and a seasonal effect with a 24week cycle, and the sum of seasonal factor is approximately zero, so the seasonal model could apply to the forecast.
From the above method, we get the maximum amount of products that can be produced by the manufacturer in the current week in the case of sufficient raw materials compared with the production capacity. For different material suppliers, except for some extreme cases where the supply is always large and close to 0, the quantity of supplied materials mostly obeys the production law of the manufacturer, i.e., there is a cyclical supply law with a 24week cycle. We demonstrate these advantages through using supplier S007 as an example:
To determine and separate the seasonal components, seasonal indices need to be calculated to determine the seasonal effects in the time series and assume weekly supply of suppliers be the forecast object indicator . We linearly combine the forecast object indicator as follows: where is the seasonal term, is the linear fit term, and is the error term. Using approximately six months (24 weeks) as a period, the 240 weeks of data were divided into 10 groups. We define the average value of each group as set of features from class :
Hence, we know that is , the value of week , group I, and is the average value of group .
We define the season term of that series is as follows from class .
After obtaining the seasonal term index, the series with the seasonal term removed can be obtained which is , and after sorting this decomposed series , a regression fit is performed on this array .
3.2. Model Definition
We define these variables which are required in the following models: (i). The supply of the supplier in the week j(ii). The zeroone variable of whether the supplier in week supplies(iii). The original material storage at the beginning of week
Next, a nested weighted recurrent network [9] can be used to investigate supply chain network coherence. Based on the results of the time series analysis, the weekly raw material ordering forecast for each seasonal cycle can be obtained, and based on this forecast and the manufacturer’s production capacity, the weekly material stock quantity state transfer equation can be established in (4). (i)(ii) is the material consumption of week ()(iii)
The optimization strategy uses a combination of continuoustime model and multiobjective optimization model to measure the uncertainty constraint and to ensure the reliability of the optimal design through the hybrid optimization strategy [10] which is an important step of our pipeline.
3.3. Establishment of Flexible Constraints for Multiobjective Planning
In this section, for a complex optimal supply problem, the complexity of the working conditions and uncertainty quantification as well as the relationship between objectives and constraints are fully considered. In addition, since multiple objectives need to be balanced, multiobjective optimization is required to make decisions by balancing the objectives when the constraints are satisfied.
Objective 1. Excessive suppliers may lead to cumbersome handover, so it is necessary to ensure that as few suppliers as possible are selected to meet production requirements. The price ratio of A, B, and C raw materials in the data set is 0.72 : 0.726 : 0.72. Take the example of group A suppliers who supply A raw materials in the data set, group A contains 146 suppliers, and the plan will cover the next 24 weeks. The zeroone variable measures whether the supply is provided in (5).
Objective 2. To minimize the ordering cost, the quantities of raw materials A, B, and C ordered are set to , , and . Since in practice, the unit prices of raw materials A and B are 20% and 10% higher than those of raw materials C, respectively, and the price of C is set to be the unit price, this optimization objective can be expressed as (6).
Objective 3. To minimize the relative difference between the products that can be produced with the order quantity and the production capacity, the difference is defined as when the products that can be produced with the weekly intake are greater than the weekly production capacity, and as when the products that can be produced with the weekly intake are less than the weekly production capacity. The objective of keeping this relative difference as small as possible is established:
3.4. Linear Weighted Transformation for Multiobjective Planning
The three objective functions are dimensionless and normalized using the linear weighting method:
3.5. Solution for Dynamic Multiobjective Planning
The difficulty of multiobjective optimization is affected by the change of Paretooptimal solutions (PSs) and different ways of its moving. We commonly use population evolutionary algorithm to find the PS. By the end of evolutionary, the diversity of population will loss and then make a decline in the adaption competency. In this study, we use the evolutionary algorithm combining with dynamic strategy (including the RS and AS) to significantly improve the convergence rate.
3.5.1. Reset Strategy (RS)
When the environment changes, we use a small amount of new environment information earlier to predict the possible movement direction of PS in the new environment [11]. Then, we also reinitialize the population by using the estimated direction and local search, making it close to the Paretooptimal solution set in the new environment.
Half of the crowd is generated by the direction of PS movement estimated by the new environmental information, and the other half is generated by a local search of the current crowd. is the noninferior solution obtained by the algorithm in the th environment, and the moving step between two environmental changes is the Euclidean distance between the centroids of two noninferior solution sets.
This distance reflects the movement direction of PS in the new environment. Through the domain search of individuals in , we obtain the predicted distribution set and the noninferior order depending on individual adaptation. Then, we notice that we can obtain the predicted movement direction by calculating the centroid of the set.
We design the corresponding moving step by integrating historical environmental moving distance and current environmental information. A new evolutionary population can be obtained by screening half of from the crowding distance.
3.5.2. Adjust Strategy (AS)
We investigate the PS with higher speed of convergence by adjust strategy (AS). After obtaining more information of the new environment, we adjust and update the current population to make more individuals close to the Paretooptimal solution set of the new environment in the current population and to accelerate the speed of convergence.
When the new environment information increases, we obtain the PS movement direction by the current group’s center of mass .
Similar to the RS strategy, we train to generate estimations of new members of the population and replace underperforming individuals by crowding distance. is a complement to the diversity of population in the new environment during AS strategy.
is random normal function, is , and .
By combining traditional evolutionary algorithms with adjustment strategies, we obtain solutions approaching to PS in the new environment. Since we predict the position of Paretooptimal front (PF) more accurately, we can also obtain the PS of the set more accurately, which content the requirements of algorithm design.
4. Result
4.1. Practical Ordering Example
4.1.1. Time Series Forecasting
In this section, we investigate the feasibility of each module in the proposal method. Considering the class A supplier S007 in the dataset to confirm our thoughts, the time series analysis is performed according to this pipeline, and the final fitted function is shown in Figure 1. After removing the effect of the seasonal term, the data is fitted as a gently rising straight line, from which the trend of supply and demand can be known.
The data series can be further verified to be a seasonal series based on the comparison between the series with the seasonal term removed and the original series.
As shown in Figure 2, the original data are more consistent with the predicted values, so the prediction of the data combines the strength both accurate and reliable. Therefore, it is feasible to transform the original series into a smooth series through the loworder difference and then use the ARMA model to fit the smooth series, which can provide information on the seasonal development pattern required for the model building. After removing the seasonal factor, since the lag coefficient of ARMA tends to 0 at this point, it is known that the data series is smooth.
The trend predicted by ARMA gives a rough limitation of the availability of suppliers, which is further used as a constraint in the multiobjective optimization process, and the model can be solved as described in Methods.
The trend in time series shows a significant increase in supply availability in weeks 1 and 5 but also shows that in most weeks, the supply is less than the supply limitation, and there is a hidden danger of oversupply for the producer. Therefore, in order to ensure that inventory levels are always no less than the number of raw materials stocked to meet the twoweek production demand and to even out the transportation pressure across weeks, more stock can be predetermined to be stocked in weeks 1 and 5.
4.1.2. Rigid Requirements for Multiobjective Planning
The supply is replenished at the peak 1 and 5 weeks of the forecasted incoming stock; the weekly availability of the three raw materials is defined; one of the offset variables and must be 0 is affirmed. Inspired by these constrains, we develop a new algorithm to solve such multiobjective planning problem.
4.1.3. Result of Multiobjective Planning
In this section, we address the problem arising from the tradeoff between objectives and constraints. The result of multiobjective planning was solved by linear weighted transformation. Through the multiple search algorithm, it is found that the best results are obtained at 7 : 2 : 1 by using plane separation [12]. Based on the allocation of the supply selection status for that week, the supply efficiency is calculated using the total supply forecast table as follows in (10). where is the total amount of materials ordered for the th supplier in week and is the total amount of materials provided by the th supplier in week . The resulting utilization rates for each state are shown in Table 1 which presents the manufacturer’s utilization rate for raw materials for the next 24 weeks, which shows that the utilization rate for materials purchased according to the ordering plan is basically greater than 93%.
4.1.4. Feasibility Verification of Result
The industry has a monopoly or dominant position in the larger supply of suppliers. The topranked suppliers are selected by CRICTIC importance assessment and compared with the planning results obtained from the multiobjective optimization strategy and provided that the two can be verified with each other, and it means that the selected strategy has feasibility. To verify that, the following results are proposed to confirm our proposal.
A matrix of size is constructed, with representing the number of suppliers and representing the number of evaluation indicators. By dimensionless processing, all indicators are normalized to positive indicators. We investigate this method to consider the standard deviation between indicators to measure their difference fluctuation and the correlation between indicators to solve the endogeneity problem. We complement our method with considering the size of the information entropy of indicators for weight assignment, and get the weight score ranking by program. Then, by comparison, it can be seen that the manufacturer material sources are mainly the top 50 suppliers in this ranking such as in Table 2.
4.2. Dynamic Multiobjective Change in Practical Ordering Example
Considering the well solved on prior constraints, we further use the more challenging changing situation to evaluate the proposal method. Assume that the ordering scenario has a new dynamic due to the larger volume of material C, which requires more transshipment resources. The emergence of new objectives and changes in the original objectives lead to a consequent change in the optimization, which is the result of the combined conflicting objectives of the optimization algorithm [13]. To save transshipment resource plan to purchase as much as possible of material A and as little as possible of material C, set the quantity of material A ordered to , the quantity of material B ordered to , and the quantity of material C ordered to . Optimization objective one is replaced by the following:
The approximated POF is searched based on similar scenarios by a multiple search algorithm which referring some similar multiobjective optimization methods [5, 14, 15].

The optimal allocation ratio can be solved by this program as the algorithm in Algorithm 1, in this case, changes to 0.80 : 0.15 : 0.05.
We investigate our method with contrastive learning to further combat the results of the two optimizations in Table 3 which compares the different order quantity between the second ordering scenario and the first ordering scenario after changing the optimization objective, reflecting the consideration of dynamic factors for this scenario. Solution for multiobjective optimization is the fundamental problems in this study. Our algorithm has proven to be highly effective and efficient in dynamic multiobjective ordering strategy because of the saved storage costs. The empirical results show that the method satisfies the manufacturer’s need for planning for ordering materials, satisfies reasonable strategy changes under target differences and can recommend for the planning of material ordering problems.
4.3. Practical Transport Example
4.3.1. Model Definition
The process of transshipment is similar despite the different optimization objectives and can be modeled analogously: (i). The amount of material available from supplier in th week(ii). Whether supplier provides material in week (01 variable)(iii). Whether supplier uses forwarder to deliver material (01 variable)(iv) and are decision variables
Optimization objective
Requirements for multiobjective planning
Due to the constraints mentioned in this section, the solution of this model requires a tradeoff between the objectives and the limitations, and the same algorithm as above can be used to solve for the optimal transportation method. Considering its variable correlation, the yield can be expressed by a topological network of recursive form [16]. On this basis, the maximum capacity of yield can be accommodated is inverted by setting the weekly capacity to and the weekly production to as follows:
4.3.2. Result of Multiobjective Planning
The specific transporters corresponding to each supplier in regular situation are allocated according to Tables 4 and 5, which allows for the most rational allocation plan and a stable partnership between the suppliers and the transporters. Tables 4 and 5 show the weekly volume of each transporter where 0 means no shipment and the corresponding value represents the shipment of the corresponding value. We notate the utilization efficiency and quantify it by the data in Tables 4 and 5.
This equation characterizes the ratio between the average and maximum transshipment volume of the th transporter, where refers to the week of the th transporter in the equation, and represents the total number of weeks carried by the transporter. We obtain that the efficiency of each transporter is greater than 60%. In particular, the efficiency of transporter 2 is as high as 91.31%, which shows the high efficiency of the model and the improvement over the current stateoftheart.
5. Conclusion
The study assisted in the establishment of the ordering and transportation strategy, and a time series prediction model was established to forecast the corresponding material supply based on the historical data of manufacturer and supplier capacity and supply, analyzing its trend and seasonality. At the same time, considering the number of suppliers and material procurement costs, dynamic multiobjective programming model is established. Then, multiple search algorithm using a minority new environment information and the old environment information to adjust the population with candidate solutions is used to solve this problem thus improving the convergence speed of the algorithm to meet market demand and business needs when a change in the environment is detected. The experience in the paper considers the minimum inventory strategy and empirically proves the good adaptability of the model through the establishment of a hybrid strategy model. Through the comparison of ordering strategies under dynamic target changes, it is found that the method provides decision reference and theoretical support to the customization of raw material ordering strategies for actual projects.
The empirical results of the ordering strategy show that, first, the positive effects of material storage space and seasonal factors on the supply chain system should be increased to facilitate centralized decisionmaking for greater benefits; second, the dynamic ordering strategy model responds to the digital representation of demand and has a positive effect on the overall efficiency improvement of the supply chain. The prediction process in our study can be similar to the prediction of coke price [17] combining with the deep learning model to improve the consideration of internal and external variables. Moreover, this paper does not consider the impact of nonimmediately decaying materials on ordering, storage, and preservation, which will be further considered in future studies.
Data Availability
The original data of this paper comes from the data of competition question C for CUMCM2021(China University mathematical modeling competition), which can be downloaded publicly from the official website.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Authors’ Contributions
Kejia He contributed to conceptualization, data curation, writing—original draft, and writing—review and editing. HongYu Cheng contributed to software and visualization. Yuchen Zhou contributed to methodology and project administration. CuihuaXie contributed to funding acquisition, supervision, and validation.
Acknowledgments
This study was supported by a project grant from Zhejiang Philosophy and Social Sciences Foundation (Grand No. 19NDJC145YB).