Skip to content
BY 4.0 license Open Access Published by De Gruyter Open Access April 11, 2019

Optimizing model and algorithm for railway freight loading problem

  • Hao Chen , Yugang Wang EMAIL logo and Raquel Martínez Lucas
From the journal Open Physics

Abstract

This paper constructed an optimizing model on multiple vehicles and multiple goods, which is used in a study on the loading problem for railway freight. We take the maximum utilization coefficient of car loading capacity, volume capacity and layout optimal degree as objective functions. Several major factors in railway transportation have been cited as constraints: center of gravity, non-overlapping freight, transportation circumscription, car marked loading capacity and volume capacity. In order to increase the car utilization coefficient and obtain flat layers in the loading scheme, an improved genetic annealing algorithm is proposed to solve this optimizing model. The proposed algorithm can efficiently obtain a satisfactory loading scheme in railway transportation. At the end, a numerical example shows that the proposed model and algorithm are better than First Fit Algorithm and Neural Network Algorithm in goods loading capacity and volume capacity, the new loading schemes are more flat on each layer.

1 Introduction

The railway freight loading problem is a common problem in the process of goods transportation. A reasonable loading scheme can reduce the unit costs and enhance the working efficiency, but also increase the safety of the transportation to a certain extent.

The loading problem means arranging multiple goods on several different cars, which aims to obtain the optimum and reasonable loading scheme. Moreover, the similar optimizing algorithm holds for highway, aviation, container and other modes of transportation [1]. In the actual transportation, there are a lot of constraints in the railway freight loading problem, so the heuristic approach is more used to obtain a satisfactory solution in solving this kind of NP-hard problem [2].

Lots of studies for loading problem have been done by domestic and international scholars. Researches abroad in this aspect have an early start. In this paper, a systematic review of the current state of research in loading problem is presented. The G&R algorithm with a concept of “layers” proposed by [3], is the foundation of current stratification algorithms. [4] Constructed a space matrix that can describe the free space and the loaded space simply and intuitively. A new approximation algorithm presented by [5] for the bin-packing problem, which has a linear running time and an absolute approximation factor: 3/2. The algorithm performs well even the number of bins is small and has a low complexity and a short computation time.

There are few reports concerning about this research in domestic, which pay more attention to the algorithm (e.g, heuristic algorithm or artificial intelligence algorithm) [6]. Took the maximum loading capacity and volume capacity as optimizing targets. An optimizing model of freight loading problem for one car with multiple goods and multiple cars with multiple goods was proposed. [7] used the approach of spatial merging and small block reservation, the freight loading optimizing algorithm was built on the basis of central skeleton and the moving strategy. According to the characteristics of the container, [8] and [9] solved the optimal layout problem of rectangular objects in three-dimensional space. Although the algorithm has a low optimizing rate, the calculation process is simple and the loading efficiency is higher [10]. Further studied the loading problem of irregular goods, performed the combinatorial optimization on L, T and U shape goods, and combined irregular shape goods into regular shapes, such that the research results are more suitable for the actual situation.

In summary, the existing researches on the freight loading problem are mostly focus on the maximum loading capacity and volume capacity,which lead to the lack of further studies on loading scheme and Layout optimal degree. Besides, how to combine the two objective functions loading capacity and volume capacity better, they have not yet involved. In this paper, an optimizing model for railway freight is proposed, we take the maximum utilization coefficient of car volume and loading capacity and layout optimal degree as objective functions, and the goods gravity central shift, non-overlapping freight and transportation circumscription are cited as constraints. Reconfiguring and optimizing a group of goods in numerical example. To obtain the optimized scheme of railway freight loading problem, an improved genetic annealing algorithm is discussed by the characteristics of the proposed model.

2 Freight loading model

We first build the model of car body and goods by the vertexes coordinate sets, and further present the freight loading model. Specific definitions are as follows.

2.1 Car coordinate system

The car coordinate system (shown in Figure 1) is denoted by V i = O 0 , D i + , O0 = (0, 0, 0) and D i + = X i , Y i , Z i are the start-point and end-point of the car coordinate set, respectively.

Figure 1 Car coordinate system
Figure 1

Car coordinate system

2.2 Freight layout parameters

We define that H = (H1, H2, · · · , Hn) as a good set. for Hj, it has a vertex set Hj = (Oj , Dj), where Oj (Oj = (xj , yj , zj)) is the closest point to the car starting point O0. By the same token, D j = x j , y j , z j . Then Hj can be denoted by the combination of xj,yj,zj,x'j,y'j,z'j. However, considering that the shape of the goods are mostly regular rectangles, we can obtain the only spatial location of the goods by the coordinates of Oj and Dj, which is recorded as (Oj , Dj).

2.3 Freight stratification theory

In order to improve the space utilization coefficient, we divide the whole car into several layers, and choose a suitable value in ( a j = x j x j , b j = y j y j , c j = z j z j ) as the layer height δ.

After we defined a certain δ, the goods will be loaded in the corresponding layer by an allowable height difference: cj (δξ , δ + ξ). ξ is determined by the layer height and the goods height. In the actual loading process, the height difference between different goods can be filled by padding wood, so as to ensure the same height of the goods layer, and make the goods placed neatly and stably. ΣAt this point, there is a set of layers k}, and(δk + ξ k)≤ Zi.

2.4 The representation of relative position for two goods

Here, Hst denote the relative position between s and t. If s and t have contact parts in space, then H s t = O t D s ¯ x O t D s ¯ y = 0. We denote an unfeasible scheme of overlapping goods by O t D s ¯ x > 0 O t D s ¯ y > 0 , where O t D s ¯ x and O t D s ¯ y are the projections of OtD onto the X and Y axis, respectively, as shown in Figure 2.Meanwhile, under the constraint of the stratification theory, good height differences are within |2ξ|. Thus, we can obtain a better scheme with a bigger ‖Hst.

Figure 2 Freight model parameters
Figure 2

Freight model parameters

2.5 Genetic annealing algorithm

Based upon past experience, several parameters Φ, pc, pm and T of the genetic algorithm are defined to solve the problems.

It is worth to note that, if the average volume of goods is small, even take a stratification approach, the number of goods loaded in car is still very large, such that with the goods number increase, the number of layout schemes also increase explosively. Since the complex calculations and slow convergence, original algorithm rarely obtains the optimal solution. Therefore, we propose an improved genetic annealing algorithm for the railway freight loading problem.

3 Railway freight loading model

3.1 Model definitions

In order to reduce the complexity of the model, some influential factors in the actual situation are defined as follows.

Definition 1 The goods can be placed in any direction, and the final placement must be parallel to the axis. Goods deformation can be ignored.

Definition 2 Goods and the outer packing are strong enough and can be stacked. If there is a suspended part below a good, we consider that a support is placed below without weight. The size of single good should be more than 0.02m3 according to the requirements of less-than-carload freight, that is, | (x jxj) (y jyj) (z jzj) | ≥ 0.02m3

Definition 3 The gravity center of goods is considered be the geometric center: Gj = ((x jxj)/2, (y jyj)/2, (z jzj)/2). Car weight and loading capacity see Regulations for Railway Freight Traffic for the details.

Definition 4 The studied goods are non-dangerous goods, or can be assorted packing, which will be sent and arrived at same station without loading and unloading.

Definition 5 Let η = (0, 1) be a determination coefficient. We define k = (xk , yk , zk) as a point, k ∈ Hs and k ∈ Ht indicate that s is in contact with t, then ηst = 1; if not, ηst = 0.

3.2 Constraints

The constraints of railway freight transportation are complex, especially in multi-type goods transportation. Consideration should be given to the overall constraints of the loaded car. The freight loading scheme, in addition to meet the Regulations for Railway Freight Traffic (hereafter this text will be abbreviated as Freight Regulations) and the Regulations for Loading and Strengthening of Railway Freight (hereafter this text will be abbreviated as Loading Regulations), but also need to meet the car restrictions. They are defined in the following equations (1)-(7) (αh and αz see Loading Regulations for equation details):

(1) max x j X i
(2) max y j Y i
(3) δ k + ξ k Z i
(4) δ k ξ k min z j max z j δ k + ξ k
(5) [ ( x j + x j ) g j / 2 ] g j X i / 2 α
(6) [ ( y j + y j ) g j / 2 ] g j Y i / 2 α
(7) ( z j + z j ) / 2 + h i g j + G i h i g j + G i 2000 m m

In addition to the above constraints, it should be satisfied with the actual situation and logical constraints, such as non-overlapping in space or overloading constraint. We obtain equations (8)-(10)

(8) g j Q i
(9) O t D s ¯ x 0 O t D s ¯ y 0
(10) x j x j y j y j z j z j 0.02 m 3

3.3 Objective functions

Generally speaking, the most important issues of railway freight transportation are the safety and efficiency. Existing literatures are mainly focus on the utilization coefficient of car volume capacity and loading capacity [11]. Based on the further analysis of the full-load situation, this paper introduces the layout optimal degree as a new objective function, and proposes a multi-object model.

3.3.1 Objective function W1

The utilization coefficient of car loading capacity. The calculation of W1 involves the type, quantity and weight of the loaded goods, we can obtain equation (11)

(11) max W 1 = ω i / m = 1 m i = 1 m μ j i g j Q i i = 1 , 2 , , m

where μ j i = 0 , 1 , if j is put on i, then μ j i = 1 , if not, μ j i = 0 .

3.3.2 Objective function W2

The utilization coefficient of car volume capacity. Similar to (11), we can obtain equation (12)

(12) max W 2 = σ i m = i = 1 m μ j i x j y j z j X i Y i Z i m

3.3.3 Objective function W3

Layout optimal degree. There are obvious pros and cons of the different arrangement of goods on the car (as shown in Figure 3), Plan A is better than Plan B and Plan C. Then we introduce layout optimal degree to denote the pros and cons of arrangement with respect to all loaded goods. We obtain equation (13)

Figure 3 Comparison of different relative positions
Figure 3

Comparison of different relative positions

(13) max W 3 = H s t = s t η s t O t D s ¯ x O t D s ¯ y η s t

Based on the analysis above, the optimizing model of railway freight loading problem is constructed as follows:

max W 1 = o m e g a i m = i = 1 m μ j i g j Q i m ;
max W 2 = σ i m = i = 1 m μ j i x j y j z j X i Y i Z i m
max W 3 = H s t = s t η s t O t D s ¯ x O t D s ¯ y η s t

such that

max x j X i ;
max y j Y i ;
δ k + ξ k Z i ;
δ k ξ k min z j max z j δ k + ξ k ;
[ ( x j + x j ) g j / 2 ] g j X i / 2 α ;
[ ( y j + y j ) g j / 2 ] g j Y i / 2 α ;
( z j + z j ) / 2 + h i g j + G i h i g j + G i 2000 m m ;
g j Q i ;
O t D s ¯ x 0 O t D s ¯ y 0 ;
x j x j y j y j z j z j 0.02 m 3 ;
η s t = 1 k H s k H t 0 e l s e ; i = 1 , 2 , , m ; j = 1 , 2 , s , t , , n .

4 Improved genetic annealing algorithm

In this paper, the improved genetic annealing algorithm (short for GA-SA) is used to solve the problem: Firstly, a group of individuals is produced by selection, crossover and mutation operations in Genetic Algorithm (short for GA). Then we can obtain the new individuals in offspring by the simulated annealing. In the above multi-objective optimizing problem, it is difficult to find a set of solutions, which can satisfy all sub goals. As a consequence, we only need a Pareto Solution of the optimizing model.

4.1 Model parameters

4.2 Combination of objective functions

Two target functions of W1 and W2 are combined by weight, which can be expressed as equation (14)

(14) max U 1 = max α W 1 + β W 2

For the weight coefficients, referring to the actual experience, we denote the weight by the volume ratio and loading ratio for goods and car [1], shown in equation (15)

(15) α β = i j μ j i g j i j μ j i x j y j z j i Q i i X i Y i Z i h i

For the objective function 3 and equation (14) there is no direct connection between the two forms, so the possibility-satisfiability method is used to transform the multi-objective programming, we obtain equation (16)

(16) λ , θ = 1 (completely satisfied with the index) λ , θ 0 (completely dissatisfied with the index)

Equation (14) is a fmax function, we obtain equation (17)

(17) λ ( r ) = 0 ( r r B ) / ( r A r B ) 1

Equation (3) is a fmin function, we obtain equation (18)

(18) θ ( r ) = 0 ( r B r ) / ( r B r A ) 1

4.3 Genetic algorithm parameters

Let the population and crossover probability denote by M and pc, respectively, H is the group of goods be loaded in same layer, then M and pc can be defined by floating value [12, 13, 14, 15, 16, 17, 18]. The size of M will change with the number of goods in H. In order to retain the excellent individuals, we define pc be inversely proportional to the individual fitness, we can obtain equation (19)

(19) p c = [ f + f c t ] / [ f + f ] , f + f r a n d b e t w e e n 0 , 1 , e l s e

4.4 Algorithm steps

4.4.1 Parents genetic coding

If cj (δξ , δ + ξ), then we put cj into H from H, (xjx j) · (yjy j) is the area occupied by cj. If Σ( xjx j) · (yjy j) > XiYi, H is loading complete, if not, then expand the range of (δξ , δ+ξ) until it meets the scope of H. At this time, the initial parent generation is defined by C0 = (0, 0, · · · , 0) and C0 = (1, 1, · · \ , 1), respectively, such that the offspring, product in genetic process, can cover the scope of the solutions as much as possible.

4.4.2 Population fitness

The values of *the objective functions are used as the population fitness Pt (t is the generation number).

4.4.3 Crossover and mutation operations

Let the crossover probability and mutation probability be denote by pc and pm, such that the crossover and mutation operations of the chromosomes for parents can be denote by C t+1Ct and C t+1Ct, respectively, then we get a new offspring.

4.4.4 Simulated annealing calculation for new offspring

Let Ct and Ct be two initial generation groups, Ct+1 and Ct+1 be two second generation groups, then carrying out simulated annealing and genetic operations by initial generation probability ρ and second generation probability 1 − ρ, we get the final individuals in new offspring. Where

we define ρ in equation (20)

(20) ρ = 1 1 + exp [ ( f t f t + 1 ) / T ]

4.4.5 Fitness calculation of new population

The same as 4.2.2, population fitness of second generation (t + 1 generation) is iterative calculated until the end condition or the end generation are satisfied.

5 Numerical example

5.1 Simulation results of freight loading

In reference to examples in [1] and [2], we set up the goods data, as shown in Table 1. Let us assume that the gravity center of goods is the geometrical center. We choose the type of C62 as the study car, and its technical parameters are as follows, tare mass: 21.7t, marked loading capacity: 60t, car body volume: 12500mm×2900mm×2000mm, car floor height 1211from rail top: 1083mm, center height of gravity for empty car: 1000mm, lateral shift for goods gravity

Table 1

Goods loading data

No. Goods code Length (mm) Width (mm) Height (mm) Weight (t) Quantity
1 1-4 2 800 2 000 1 200 4.10 4
2 5-16 1 000 1 000 1 000 0.20 12
3 17-20 3 000 800 2 000 1.10 4
4 21-34 4 800 900 1 000 3.80 14
5 35-42 2 900 800 900 0.40 8
6 43-50 1 900 1 200 900 2.10 8
7 51-60 1 250 600 800 0.30 10
8 61-70 940 400 500 0.33 10

center: 100mm, longitudinal shift for goods gravity center can be defined by loading capacity, and the height of gravity center for loaded car must be less than 2000mm.

By calculation, we can obtain the first and second layer of the first car: c 1 1 = 1000 , 1200 and c 2 1 = 800 ; the first and second layer of the second car: c 1 2 = 900 , 1000

and c 2 2 = 800 , 1000 . The specific layout is shown in Figure 4 and Table 2

Figure 4 Representation of layers
Figure 4

Representation of layers

Table 2

Usage for loaded car

Single weight First layer Second layer First layer Second layer Total
Quantity Weight Quantity Weight Quantity Weight Quantity Weight
1 4.1 4 16.40
2 0.2 2 0.40 3 0.60
3 1.1 3 3.30
4 3.8 1 3.80 6 22.80 6 22.80
5 0.4 4 1.60
6 2.1 5 10.50 2 4.20 1 2.10
7 0.3 8 2.40 2 0.60
8 0.33 4 1.32 4 1.32 2 0.66
Total 30.70 8.62 28.72 26.76 94.8

As a consequence, the improved genetic annealing algorithm obtain a superior solution, which ensures better objective functions and speeds up algorithm converging.

5.2 Discussion

We compared the First Fit Algorithm (see e.g [1]) and Neural Network Algorithm (see e.g [2]) with the proposed algorithm in a same numerical example, as shown in Table 3.

Table 3

Comparison of different algorithms

Basic technical parameters Neural Network algorithm First Fit algorithm Algorithm in this paper
First car Second car First car Second car First car Second car
Longitudinal shift tolerance (mm) 1 053.0 914.6 1 053.0 914.3 1 106.3 354.4
Longitudinal shift (mm) 81.3 278.4 81.3 278.4 7.48 132.48
Lateral shift tolerance (mm) 100.0 100.0 100.0 100.0 100.0 100.0
Lateral shift (mm) 9.6 43.6 9.6 43.6 78.1 10.33
Gravity center height for car loaded (mm) 1 654.0 1 729.0 1 654.0 1 729.0 1 575.1 1 760.9
Loading capacity (t) 42.7 49.2 44.8 49.2 39.32 55.48
Utilization coefficient of car loading capacity (%) 71.00 82.00 75.00 82.00 65.53 92.47
Total volume of the loaded goods (m3) 66.678 67.048 66.684 67.532 69.764 65.324
Utilization coefficient of car volume capacity (%) 91.97 92.48 91.98 93.15 96.23 90.10

We note that, the basic technical parameters of proposed algorithm are better than First Fit Algorithm and Neural Network Algorithm in many respects. And 3.27%, 0.64% increase in utilization coefficient of car loading capacity; 1.02%, 0.65% increase in utilization coefficient of car volume capacity for details, respectively. Moreover, the new loading schemes are more flat on each layer. By objective function equation (13), the lateral and longitudinal shift of goods gravity center were reduced by 44.69%. As a consequence, the improved genetic annealing algorithm proposed above can efficiently obtain the optimal solution of the railway freight layout problem.

6 Conclusion

In order to increase the car utilization coefficient and obtain flat layers in the loading scheme, we introduce the layout optimal degree, couple with the combination of car loading capacity and volume capacity, to construct the railway freight optimizing model. The following findings were obtained from this research.

The value of layer height δ and the acceptable horizon ξ directly affect the computational complexity and scheme results, reasonable choice of which are very important.

The weight ratio and possibility-satisfiability method for multi-objective function can express the best form of the solution. By simplifying the objective functions, compared with literature 1 and 2, the optimal solution increased by 3.27% and 0.64% in utilization coefficient of car loading capacity, respectively, and increased by 1.02% and 0.65% in utilization coefficient of car volume capacity, respectively.

The stratification theory optimizes the basic parameters of loaded car. Compared with the First Fit Algorithm, the average longitudinal shift of goods gravity center reduced from 179.85mm to 69.98mm, and the average longitudinal shift increased from 26.6mm to 44.215mm, but still in the allowed range.

On the basis of stratification theory, the upper and lower exchanges between layers can make the gravity center height lower, this method reduces the average height by 23.5mm compared with the previous algorithm.

Vector product is used to express the layout optimal degree. To some extent, the layout optimal degree describes the pros and cons of the different arrangement. But the degree is out of proportion to the position relationship, and it can only indicate the advantages and disadvantages.

Furthermore, since there are too many parameters involved in the problem, this paper makes a simplified design for the actual situation, especially the rectangular outer packing for all goods. But in fact, the shapes of the goods transported by railway have great randomness and strong irregularity, and further research is needed.

Acknowledgement

The research described in this paper was jointly supported by the Science and Technology Project of China Ministry of Railways (grant 2013X014-E, 2012S14071 and 2009X012-I) and the Fundamental Research Funds for the Central Universities (grant SWJTU09BR135). These supports are gratefully acknowledged.

List of notation

1 − ρ is the annealing probability of second
generation
α and β are the weight of W1 and W2;
δ is the layer height;
η and µ j i are two determination coefficients;
λ and θ are the possibility-satisfiability of equations
(3) and (14);
Φ is the final evolutionary generation;
ρ is the annealing probability of initial generation
ξ is the acceptable horizon of good heights
when we load goods into a layer;
D i + is the end-point of the car body coordinate
set;
f (ct) is the fitness value with respect to an individual
ct;
H is a good set;
H is the group of goods be loaded in a same
layer;
h i is the plane height of car body;
i is the car sequence number;
M is the population;
r is the specific value of target function;
T is the simulated annealing temperature in
initial state;
t is the generation number;
αh is the lateral shift for gravity center of
goods;
αz is the longitudinal shift for gravity center
of goods;
δk is the layers set;
ωi is the coefficient of utilization for car loading
capacity;
σi is the utilization coefficient of volume capacity
for i;
x and y are the projections of OtD onto the X and
Y axis;
aj ,bj,cjare the length, width and height of
H j;
C0 and C 0 are the initial parent generations;
Ct and C t are two initial generation groups
C t + 1 Ct+1 and are two second generation groups
Dj is the loading end-point of H j with respect
to the loaded car;
f + and f are the maximum and minimum value of
individual fitness;
ft and ft+1 are the fitness of Ct+1 and C t + 1 .
Gj is the gravity center;
gj is the weight of j ;
hi is the gravity central height of empty car;
Hj is the good in H, (Hj ∈ H);
Hst is the relative position value between s
and t;
O0 is the start-point of the car body coordinate
set;
Oj is the loading start-point of H j with respect
to the loaded car;
OsOt and DsDt are start-point and end-point coordinates
with respect to s and t;
pc is the crossover probability;
pm is the mutation probability;
Pt + is the population fitness;
Qi is the maximum loading capacity of car i;
rA and rB are the ideal index values that are completely
satisfactory and dissatisfactory;
Vi is the car coordinate system;
Xi , Yi , Zi are the total loading limits of the goods
on the X, Y and Z axis with respect to the
loaded car;

References

[1] Guo Y., Chen Z., Tang B., Optimization model and algorithms for loading and packing multi-category freights into multi-vehicle, J. China Railway Sci., 2011, 32(03), 111-116.Search in Google Scholar

[2] Xu Q., Research on multiple vehicle and multiple freight loading problem based on hopfield network, Logist. Technol., 2010, 33(12), 101-104.Search in Google Scholar

[3] George J.A., Robinson D.F., A heuristic for packing boxes into a Container, Comput. Oper. Res., 1980, 7, 147-156.10.1016/0305-0548(80)90001-5Search in Google Scholar

[4] Ngoi B.K.A., Tay M.L., Chua E.S., Applying spatial representation techniques to the container packing problem, Int. J. Prod. Res., 1994, 32, 111-123.10.1080/00207549408956919Search in Google Scholar

[5] Berghammer R., Reuter F., A linear approximation algorithm for bin packing with absolute approximation factor 3/2, Sci. Comput. Prog., 2003, 48, 67-80.10.1016/S0167-6423(03)00011-XSearch in Google Scholar

[6] Lei D., Tang B., Zhang Y.Optimizingmodel and algorithms of loading and packing multi-piece freights into one car. J. Chin. Railway Soc., 2011, 33(9), 1-9.Search in Google Scholar

[7] Lei D., Onoptimizing theory an application of goods loading and carrying, 2005, Doctoral thesis Cent. S. Univ.Search in Google Scholar

[8] Jiang Y., Cha J., He D., Research on the packing of loading rectangular freight into a container, J. China Railway Soc., 2000, 22(6), 13-18.Search in Google Scholar

[9] Pu R., Wu L., Li G., Simple and rapid optimization of packing method for single goods unconstrained three-dimensional packing problem, Technol. Innov. Manag., 2017, 38(2), 132-135.Search in Google Scholar

[10] Zhang G., The study of multi-shape irregular cargo loading optimization problem. Master’s thesis Shandong Univ., 2007.Search in Google Scholar

[11] Wang H., Li Z., Model and algorithm of optimal loading for multi-category goods and multi-vehicles, J. Chin. Railway Soc., 2010, 32(1), 93-97.Search in Google Scholar

[12] Chen Y., Tang M., Tong R., Packing of polygons using genetic simulated annealing algorithm, J. Comp. Aid, Design Comp. Graph., 2003, 15(5), 598-903.Search in Google Scholar

[13] Reddy M.V., Hayashi S., Choi D., Cho H., Chang Y., Short chain and mediumchain fatty acids production using food waste under non-augmented and bio-augmented conditions, J. Clean. Prod., 2018, 176, 645-653.10.1016/j.jclepro.2017.12.166Search in Google Scholar

[14] Awati V., Jyoti M., Homotopy analysis method for the solution of lubrication of a long porous slider. Appl.Math. Nonlinear Sci., 2016, 1(2), 507-516.10.21042/AMNS.2016.2.00040Search in Google Scholar

[15] Gao W.,WangW., New Isolated toughness condition for fractional (G, F, N) - critical graph, Colloq. Math., 2017, 147(1), 55-65.10.4064/cm6713-8-2016Search in Google Scholar

[16] Peng W., Liu Z., Motahari-Nezhad M., Banisaeed M., Shahraki S., Beheshti M., A detailed study of oxy-fuel combustion of biomass in a circulating fluidized bed (Cfl) combustor: evaluation of catalytic performance of metal aanoparticles (Al, Ni) for combustion eflciency improvement, Energy, 2016, 109, 1139-1147.10.1016/j.energy.2016.04.130Search in Google Scholar

[17] Baig A.Q., Naeem M., Gao W., Revan and hyper-revan indices of octahedral and icosahedral networks, Appl. Math. Nonlin. Sci., 2018, 3(1), 33-40.10.21042/AMNS.2018.1.00004Search in Google Scholar

[18] Halim A.A., Hanafiah M.M., Khairi A. Ammonia removal from sewage wastewater using chemically modified sand, Appl.Search in Google Scholar

Received: 2018-10-28
Accepted: 2019-01-28
Published Online: 2019-04-11

© 2019 H. Chen et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Downloaded on 29.3.2024 from https://www.degruyter.com/document/doi/10.1515/phys-2019-0014/html
Scroll to top button