In order to simplify the complex calculation and solve the difficult solution problems of neutrosophic number optimization models (NNOMs) in the practical production process, this paper presents two methods to solve NNOMs, where Matlab built-in function “fmincon()” and neutrosophic number operations (NNOs) are used in indeterminate environments. Next, the two methods are applied to linear and nonlinear programming problems with neutrosophic number information to obtain the optimal solution of the maximum/minimum objective function under the constrained conditions of practical productions by neutrosophic number optimization programming (NNOP) examples. Finally, under indeterminate environments, the fit optimal solutions of the examples can also be achieved by using some specified indeterminate scales to fulfill some specified actual requirements. The NNOP methods can obtain the feasible and flexible optimal solutions and indicate the advantage of simple calculations in practical applications.
1. Introduction
Traditional inventory models [1–4] and production planning models [5–7] involve deterministic constrained functions and/or objective functions in deterministic environments. Nevertheless, uncertainty is nearly universal in real world. Therefore, many uncertain optimization methods were proposed for optimization problems with uncertain variables, interval numbers, stochastic, and fuzzy logics [8–15]. In many applied fields, such as management, engineering, and design problems, uncertain programming has been broadly carried out so far. In order to obtain the optimal crisp values of the objective function and the optimal feasible crisp solutions of the decision variables, the constrained functions and/or objective functions are usually changed into some crisp or deterministic programming problems in existing uncertain programming approaches. So, the aforementioned transformed methods are not really meaningful indeterminate approaches because the real indeterminate optimization problems can only indicate indeterminate solutions rather than optimal crisp solutions in indeterminate environments. Nevertheless, indeterminate programming problems imply the corresponding indeterminate optimal values of the objective function and indeterminate optimal solutions for the decision variables under indeterminate environments. So, it is necessary to find some fit optimization approaches for dealing with indeterminate programming problems with indeterminate solutions.
Smarandache [16–18] is a pioneer of indeterminacy theories which provide the new minds to solve indeterminacy problems. He adopted the imaginary value denoted by I and then introduced a neutrosophic number (NN) z = x + yI for x, y ∈ R (R: the set of all real numbers) composed of the determinate part x and indeterminate part yI. As for describing indeterminate and incomplete information, obviously, NNs in the indeterminacy theories are a useful mathematical tool. With the development of indeterminacy theories, NNs were also applied to fault diagnosis [19, 20] and decision making [21, 22] under indeterminate environments.
Further, thick function or interval function named neutrosophic function, neutrosophic precalculus, and neutrosophic calculus were provided by Smarandache [23] in 2015, where thick function e: S ⟶ E(S) (E(S) is the set of all interval functions) as the form of an interval function e(x) = [e1(x), e2(x)]. The indeterminate function was applied in engineering problems successfully. For example, Ye et al. [24, 25] and Chen et al. [26, 27] proposed expressions of neutrosophic function and applied NNs in analyzing the joint roughness coefficient. Later, Ye [28] used neutrosophic linear equations of NNs to solve traffic flow problems.
At present, neutrosophic linguistic numbers, hesitant neutrosophic linguistic numbers, and their aggregation operators were applied to multiattribute decision making [29–31].
But in real situations, affected by each kind subjective and objective reasons, such as absences of precise information judged by decision makers or experts, loss of data, and measurement errors, there exist some indeterminate problems. As for the concepts of NNs, NN functions containing indeterminacy I can represent the indeterminate problems with partial certainty and partial uncertainty under indeterminate environments. Ye [32] and Jiang and Ye [33] introduced NN nonlinear and linear programming models and their preliminary solution methods. However, existing methods for solving complex NN optimization problems imply some difficulty and calculational complexity in their solution process. Inspired by the previous solution methods, this paper first selects the models of practical applications in production process, such as inventory models and production planning models. Then, NN nonlinear and linear mathematical models and their solution methods (Matlab built-in function “fmincon()” and operations of NNs) are built with indeterminacy I as our preliminary application study. Finally, real examples of NN linear programming (NN-LP) and NN nonlinear programming (NN-NP) problems illustrate the feasibility of the proposed methods. The advantage of the proposed methods is that the optimization calculations are simple and effective in practical applications.
The remainder of this paper is organized as follows. Section 2 depicts some concepts and their operations of NNs. Section 3 first introduces NN-NP problems with an inventory mathematical model and model formation and then uses two methods (Matlab built-in function “fmincon()” and operations of NNs) to solve the NN-NP problems in indeterminate setting. Section 4 presents NN-LP problems with the production planning mathematical model and model formation and then applies two methods regarding the Matlab built-in function “fmincon()” and operations of NNs to solve the solutions in the NN-NP problems and to show the simplicity and effectiveness of the proposed NN-LP methods. Conclusions and future research are provided in Section 5.
2. Mathematical Preliminaries2.1. Some Concepts and Their Operations of Neutrosophic Numbers (NNs)
The concept of NN was first proposed by Smarandache [34, 35], which consists of two parts (a determinate part and an indeterminate part). He defined the mathematical expression form z = x + yI for x, y ∈ R, where R represents all real numbers and I is indeterminacy. So, it is conveniently used in indeterminate environments.
For example, consider that a NN is z = 13 + 5I. Then, its determinate part value is 13 and its indeterminate part value is 5I. When I ∈ [0, 0.5], it is equivalent to z ∈ [13, 15.5] for sure z ≥ 13.
Let z1 = x1 + y1I and z2 = x2 + y2I be two NNs. Then, Smarandache [34, 35] gave their operations of NNs in the following:
z1 + z2 = x1 + x2 + (y1 + y2)I.
z1– z2 = x1– x2 + (y1– y2)I.
z1 × z2 = x1x2 + (x1y2 + x2y1 + y1y2)I, in particular, when z1 = 0 and z2 = I, we get the equation with 0 × I = 0.
z12 = (x1 + y1I)2 = x12 + (2x1y1 + y12)I, in particular, when z1 = I, we get the equation with I^2 = I.
z1/z2=x1+y1I/x2+y2I=x1/x2+x2y1−x1y2/x2x2+y2I for x2 ≠ 0 and x2 ≠ –y2.
3. Neutrosophic Number Nonlinear Programming (NN-NP)3.1. NN-NP Mathematical Model
The usual mathematical model of NN-NP is represented in the following form [36, 37]:(1)Min Fx,Is.t.gix,I≤ci,i=1,2,...,p,hjx,I=0,j=1,2,...,q,x∈Zn,where g1x,I,g2x,I,...,gpx,I,h1x,I,h2x,I: Zn⟶Z(Z is the set of all NNs), and I ∈ [IL, IU](the interval range of I).
3.2. Inventory Mathematical Model [38]3.2.1. Notations
The following notations are used in the inventory model.
Three decision variables:
D: demand/unit/time
Qp: production quantity/batch
Cs: setup cost/unit/time
Except the above cost variable CS, three other cost variables are
Cta: total average cost/unit/time
Ctp: total production cost/cycle
Ch: time depending on holding cost/unit/item
Other time and space variables:
T: every cycle of length
Q(t): inventory level at time t (t ≥ 0)
S: total storage space area
s0: space area/unit/quantity.
3.2.2. Assumptions
The inventory model is developed by considering the following assumptions:
Only one item is involved in the inventory system.
The replenishment occurs with the near instantaneous response.
The startup time can be ignored.
The demand rate at any time is constant.
The total production cost Ctp is related to the setup cost CS and production quantity QP.
Holding cost is the time depended function.
3.3. Model Formation
As shown in Figure 1, in every time period T, the value of the production quantity Q(t) decreases from Qp to zero. The slope of the line is constant negative D and denoted by dQt/dt=−D0≤t≤T.
Crisp inventory model.
The total average cost of the cycle T (denoted by Cta) consists of three sections: setup cost (denoted by C1), holding cost (denoted by C2), and production cost (denoted by C3).(2)Cta=C1+C2+C3.
Because we have the equation Q(t) = Qp−Dt, we obtain the cycle T, T=Qp/D.(3)C1=CsT=CsDQp,C2=∫0TChQtdtT=∫0TetQtdtT=eQp3/6D2TC3=CtpT=fCs−xQp−yT=fDCsxQp1+yeQp26D.
Based on equations (2) and (3), we obtain the following equation:(4)Cta=CsDQp+eQp26D+fDCsxQp1+y.
So, the inventory model is constructed as follows:(5)Min CtaD,Cs,Qp=CsDQp+eQp26D+fDCsxQp1+y,s.t. s0Qp≤S,D,Cs,Qp>0.
3.4. Solution Corresponding to Matlab Built-In Function “fmincon()”
In order to conveniently calculate the solutions, we simplify some parameters and set some constants with history records, where e = 18, f = 5, x = 1, y = 3, s0 = 200, and S= 1100. When we assume D = x1, Cs = x2, and Qp = x3, we can obtain the following mathematical model:(6)Min Ctax1,x2,x3=x1x2x3+18x326x1+5x1x2x34,s.t. 200x3≤1100with maximum allowable tolerance 400x1,x2,x3≥0.
Assume a1 = x1−40.496I, a2 = x2 + 0.058I, and a3 = x3−2I; then, equation (6) can be expressed in the following form:(7)s.t.200x3−2I≥1100with maximum allowable tolerance 400,x1−40.496I>0,x2+0.058I>0,x3−2I>0.
According to the de-neutrosophication technique proposed by Ye [39] and considering I = 0 or 0.5 or 1 as the minimum or moderate or maximum indeterminacy, we can obtain three optimal solutions as follows:
x1∗ = 80.615, x2∗ = 0.097, x3∗ = 7.500, and fx∗,I = 4.187 for I = 0.
x1∗ = 60.367, x2∗ = 0.126, x3∗ = 6.5, and fx∗,I = 4.343 for I = 0.5.
x1∗ = 40.119, x2∗ = 0.155, x3∗ = 5.5, and fx∗,I = 4.525 for I = 1.
Clearly, using the indeterminacy I ∈ [0, 1], different optimal results are revealed. The optimal solutions of the optimization problem are x1∗ = [40.119, 80.615], x2∗ = [0.097, 0.155], and x3∗ = [5.5, 7.5] for fx∗,I = [4.187, 4.525], which show the interval optimal ranges.
3.5. Solution Corresponding to Operations of NNs
According to the front optimal solutions, we assume a1 = x1 + y1I = 80.615−40.496I, a2 = x2 + y2I = 0.097 + 0.058I, and a3 = x3−y3I = 7.500−2I, and then we give the results by equation (9):(8)Min Ctaa1,a2,a3=a1a2a3+18a326a1+5a1a2a34=x1+y1Ix2+y2Ix3+y3I+18x3+y3I26x1+y1I+5x1+y1Ix2+y2Ix3+y3I4=x1x2+x1y2+x2y1+y1y2Ix3+y3I+3x32+2x3y3+y32Ix1+y1I+5x1+y1Ix2+y2Ix32+2x3y3+y32I2=x1x2+x1y2+x2y1+y1y2Ix3+y3I+3x32+2x3y3+y32Ix1+y1I+5x1+y1Ix2+y2Ix34+22x3y3+y32x32+2x3y3+y322I5x1x2x34+5x2x34y1−x1x222x3y3+y32x32+2x3y3+y322+y2x34+y222x3y3+y32x32+2x3y3+y322x2x34x2+y2x34+y32x3+y32x32+2x3y3+y32I.
Because a1 = 80.615−40.496I, a2 = 0.097 + 0.058 I, and a3 = 7.500−2I, we can get x1 = 80.615, y1 = −40.496, x2 = 0.097, y2 = 0.058, x3 = 7.5, and y3 = −2. Then, we calculate the three costs, respectively, as follows:
Setup cost:(9)C1=x1+y1Ix2+y2Ix3+y3I=1.407+0.087I.
Holding cost:(10)C2=18x3+y3I26x1+y1I=2.093+0.169I
Production cost:(11)C3=5x1+y1Ix2+y2Ix3+y3I4=1.047+0.082I
Then, we add the three costs and obtain the total cost Cta with equation (2) as follows:(12)Cta=C1+C2+C3=x1+y1Ix2+y2Ix3+y3I+18x3+y3I26x1+y1I+5x1+y1Ix2+y2Ix3+y3I4=1.047+0.087I+2.093+0.169I+1.047+0.082I=4.187+0.338I.
So, the calculational results validate that the same solution is obtained by using the two methods of both the Matlab built-in function “fmincon” and the operations of NNs, which are x1∗ = [40.119, 80.615], x2∗ = [0.097, 0.155], and x3∗ = [5.5, 7.5] for fx∗,I = [4.187, 4.525]. We also obtain every cost C1 = [1.047, 1.134], C2 = [ 2.093, 2.262], and C3 = [1.047, 1.129], which are the interval optimal ranges.
4. Production Planning Mathematical Model4.1. NN-LP Mathematical Model
The usual mathematical model of NN-LP is similar to mathematical model (1), so we omit it.
4.2. Production Planning Mathematical Model4.2.1. Notations
The following notations are used in the production planning model.
Nine decision variables:
a1 to a6: product quantities of six plans of type I
a7 to a8: product quantities of two plans of type II
a9: product quantities of two plans of type III
Objective function:
b: maximum profit
4.2.2. Assumptions
The production planning model is developed by considering the following assumptions:
Every product must pass two working procedures: A and B.
The startup time of two working procedures can be ignored.
Product quantities are only affected by validity time of machines.
The demand rate at any time is constant.
4.3. Model Formation
As shown in Table 1, we consider an application in production planning studied by Hu [40]. A company manufactures three types of products: Types I, II, and III. All types must pass two working procedures: A and B. We consider that procedure A can be operated on machine A1 or A2 and procedure B can be operated on the machines B1, B2, and B3. Type I can be operated on all machines of procedure A and procedure B; Type II can be operated on all machines of procedure A and only machine B1 of procedure B; Type III can be operated on only machine A2 of procedure A and machine B2 of procedure B. Our aim is to schedule the optimal production planning, which can pursue for the maximum profits. All used data are listed in Table 1, including required procedure time of every working procedure, processing fees, material cost, and selling price per unit. So, Type I has six plans to produce products, along with (A1, B1) or (A1, B2) or (A1, B3) or (A2, B1) or (A2, B2) or (A2, B3), respectively. Similarly, we consider the product quantities of the six plans a1, a2, a3, a4, a5, and a6, respectively. Type II has two plans to produce products, along with (A1, B1) or (A2, B1), and Type III has one plan to produce products, along with (A2, B2). We consider the product quantities of the remaining three plans a7, a8, and a9, respectively. So, we can get the following objective function.
So, we get the followed production planning mathematical model:(14)max bs.t. 4.5+1.71×x1+x2+x3+8+I×x7≤5600+700I,6.7+1.8I×x4+x5+x6+8.6+1.4I×x8+11−I×x9≤0.02+0.01I,5.6−0.1I×x1+x4+7.8+1.2I×x7+x8≤0.05+0.02I,3.5+2.5I×x2+x5+10+2I×x9≤0.10+0.02I,6.7+1.3I×x3+x6≤0.04+0.02I.
4.4. Solution regarding Matlab Built-In Function “fmincon()”
According to the de-neutrosophication technique proposed by Ye [37] and considering I = 0 or 0.5 or 1 as the minimum or moderate or maximum indeterminacy, we can obtain three optimal solutions as follows:
Clearly, using the indeterminacy I ∈ [0, 1], different optimal results are revealed. The optimal solutions of the optimization problem are x1∗ = [0, 0], x2∗ = [0, 778.508], x3∗ = [465.936, 625], x4∗ = [0, 0], x5∗ = [0, 677.953], x6∗ = [0, 56.452], x7∗ = [0, 146.667], x8∗ = [386.667, 474.359], andx9∗ = [0, 583.333] for fx∗,I = [762.717, 1297.389], which shows the interval optimal ranges.
4.5. Solution regarding Operations of NNs
According the front optimal solutions, we next calculate the nine relation formulas of the indeterminacy I and variables a1, a2, a3, a4, a5, a6 , a7, a8, and a9. For example, let us calculate a3 = 465.936 + 159.064I. Firstly, according to three points (0, 465.936), (0.5, 578.231), and (1, 625), we obtain the linear equation (a3 = 159.06I + 476.86). Next we amend the intercept of trend curve on the vertical coordinate. The other linear equations are obtained in the same way. So, a1 = x1 + y1I = 0 + 0I = 0, a2 = x2 + y2I = 778.508−778.508I, a3 = x3 + y3I = 465.936 + 159.064I, a4 = x4 + y4I = 0 + 0I = 0, a5 = x5 + y5I = 677.953−677.953I, a6 = x6 + y6I = 56.452−56.452I, a7 = x7−y7I = 0 + 146.667I, a8 = x8 + y8I = 474.359−87.692I, and a9 = x9 + y9I = 0 + 583.333I; then, we calculate the results of equation (13) as follows:
So, these calculational results validate that the same solution is obtained by using the two methods of both the Matlab built-in function “fmincon()” and the operations of NNs, which are x1∗ = [0, 0], x2∗ = [0, 778.508], x3∗ = [465.936, 625], x4∗ = [0, 0], x5∗ = [0, 677.953], x6∗ = [0, 56.452], x7∗ = [0, 146.667], x8∗ = [386.667, 474.359], andx9∗ = [0, 583.333] for fx∗,I = [762.717, 1297.389] and show the interval optimal ranges.
5. Conclusion
This paper first introduced some concepts and their operations of NNs with indeterminacy I. Next, we built a mathematical model with constrained conditions and then constructed the corresponding inventory model and production planning model. Finally, we obtained the optimal solutions by using the two methods of the Matlab built-in function “fmincon()” and the operations of NNs to solve the NN-NP and NN-LP problems with constrained conditions as preliminary application study in indeterminate setting. The final results show that the two methods obtained the same effective solutions, but the former needs the Matlab built-in function along with the simple calculational process, while the latter needs a lot of operations of NNs along with the complex calculational process. Some contributions in this study are that (1) different methods can obtain the same optimal results, (2) the NN-NP and NN-LP methods provided the new application ways for engineering management, (3) the NN-NP and NN-LP methods are more suitable than other ones under uncertain environments as the generalization of traditional programming methods, and (4) the two approaches can obtain the interval solutions for avoiding determinate solutions of traditional programming methods.
Obviously, the proposed NN-LP and NN-NP methods can handle indeterminate and/or determinate mathematical programming problems, which are the generalization of existing uncertain or certain linear and nonlinear programming methods. As the preliminary application study in this paper, however, there exist a lot of mathematical solution methods and proof problems along with some complexity/difficulty in the nonlinear programming problems which need to be studied further. Hence, as our future works, one is to further analyze the two presented methods of this paper from the mathematical problems, such as the convexity problem in the nonlinear programming, the stability and solution range problems regarding the changeability of NNs, and the sensitivities of NNs on the solution results, and then NN-LP and NN-NP approaches will be extended to other fields, such as engineering design and management science.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
KumarS.KumarN.LiuS.An inventory model for deteriorating items under inflation and permissible delay in payments by genetic algorithm20163110.1080/23311975.2016.12396052-s2.0-85044826904OliveiraA.OliveiraT.SeijasmacíasA.The uniform distribution product: an approach to the (Q, r) inventory model using R201845228429710.1080/02664763.2016.12755312-s2.0-85008397837MahataP.MahataG. C.Kumar DeS.Optimal replenishment and credit policy in supply chain inventory model under two levels of trade credit with time- and credit-sensitive demand involving default risk2018141314210.1007/s40092-017-0208-82-s2.0-85020067059WangsaI. D.HuiM. W.An integrated vendor–buyer inventory model with transportation cost and stochastic demand20185429530910.1080/23302674.2017.12966012-s2.0-85057133147HedenstiernaC. P. T.DisneyS. M.Avoiding the capacity cost trap: three means of smoothing under cyclical production planning201820114916210.1016/j.ijpe.2018.04.0082-s2.0-85046700965ChoH.-M.JeongI.-J.A two-level method of production planning and scheduling for bi-objective reentrant hybrid flow shops201710617418110.1016/j.cie.2017.02.0102-s2.0-85013175279JiangC.ZhangZ. G.ZhangQ. F.HanX.XieH. C.LiuJ.A new nonlinear interval programming method for uncertain problems with dependent interval variables2014238124525310.1016/j.ejor.2014.03.0292-s2.0-84901231013PergherI.De AlmeidaA. T.A multi-attribute decision model for setting production planning parameters20174222423210.1016/j.jmsy.2016.12.0122-s2.0-85009227328VeresnikovG. S.PankovaL. A.ProninaV. A.Uncertain programming in preliminary design of technical systems with uncertain parameters2017103364310.1016/j.procs.2017.01.0072-s2.0-85016001039JiangC.LongX. Y.HanX.TaoY. R.LiuJ.Probability-interval hybrid reliability analysis for cracked structures existing epistemic uncertainty2013112-11314816410.1016/j.engfracmech.2013.10.0092-s2.0-84888048727ZhangB.PengJ.Uncertain programming model for uncertain optimal assignment problem20133796458646810.1016/j.apm.2013.01.0322-s2.0-84875373434LiuB. D.ChenX. W.Uncertain multiobjective programming and uncertain goal programming201531010.1186/s40467-015-0036-6ChenL.PengJ.ZhangB.Uncertain goal programming models for bicriteria solid transportation problem201751495910.1016/j.asoc.2016.11.0272-s2.0-85007307513WangL.ChenZ.YangG.SunQ.GeJ.An interval uncertain optimization method using back-propagation neural network differentiation202036611306510.1016/j.cma.2020.113065WangL.ChenZ.YangG.An interval uncertainty analysis method for structural response bounds using feedforward neural network differentiation20208244946810.1016/j.apm.2020.01.059SmarandacheF.1998Rehoboth, MA, USAAmerican Research PressSmarandacheF.2014Columbus, OH, USASitech & Education PublishingSmarandacheF.2013Columbus, OH, USASitech & Education PublisherKongL.WuY.YeJ.Misfire fault diagnosis method of gasoline engines using the cosine similarity measure of neutrosophic numbers20158434610.5281/zenodo.22567YeJ.Fault diagnoses of steam turbine using the exponential similarity measure of neutrosophic numbers20163041927193410.3233/ifs-1519032-s2.0-84962507429YeJ.Multiple-attribute group decision-making method under a neutrosophic number environment201625337738610.1515/jisys-2014-01492-s2.0-84978099928YeJ.Bidirectional projection method for multiple attribute group decision making with neutrosophic numbers20172851021102910.1007/s00521-015-2123-52-s2.0-84947794084SmarandacheF.2015Brussels, BelgiumEuropaNovaYeJ.YongR.LiangQ. F.HuangM.DuS. G.Neutrosophic functions of the joint roughness coefficient (JRC) and the shear strength: a case study from the pyroclastic rock mass in Shaoxing City, China20162016482570910.1155/2016/48257092-s2.0-84973334352YeJ.ChenJ.YongR.DuS.Expression and analysis of joint roughness coefficient using neutrosophic number functions2017826910.3390/info80200692-s2.0-85031326589ChenJ.YeJ.DuS.YongR.Scale effect and anisotropy analyzed for neutrosophic numbers of rock joint roughness coefficient based on neutrosophic statistics201791020810.3390/sym91002082-s2.0-85031324513ChenJ.YeJ.DuS.YongR.Expressions of rock joint roughness coefficient using neutrosophic interval statistical numbers20179712310.3390/sym90701232-s2.0-85027375328YeJ.Neutrosophic linear equations and application in traffic flow problems201710413310.3390/a100401332-s2.0-85038636251YeJ.Aggregation operators of neutrosophic linguistic numbers for multiple attribute group decision making201651Springer Plus169110.1186/s40064-016-3247-52-s2.0-84990851663FangZ.YeJ.Multiple attribute group decision-making method based on linguistic neutrosophic numbers20179711110.3390/sym90701112-s2.0-85022193737YeJ.Multiple attribute decision-making methods based on the expected value and the similarity measure of hesitant neutrosophic linguistic numbers201810245446310.1007/s12559-017-9535-82-s2.0-85039845189YeJ.Neutrosophic number linear programming method and its application under neutrosophic number environments201822144639464610.1007/s00500-017-2646-z2-s2.0-85019672338JiangW.YeJ.Optimal design of truss structures using a neutrosophic number optimization model under an indeterminate environment2016149397SmarandacheF.2013Craiova, RomaniaSitech & Education PublishingSmrandacheF.2014Craiova, RomaniaSitech & Education PublishingYeJ.Neutrosophic number linear programming method and its application under neutrosophic number environments201822144639464610.1007/s00500-017-2646-z2-s2.0-85019672338YeJ.CuiW.LuZ.Neutrosophic number nonlinear programming problems and their general solution methods under neutrosophic number environments2018713KarC.MondalB.RoyT. K.An inventory model under space constraint in neutrosophic environment: a neutrosophic geometric programming approach20182193109YeJ.An improved neutrosophic number optimization method for optimal design of truss structures201814329530510.1142/s17930057185001872-s2.0-85052193866HuY.GuoY.2012FourthBeijing, ChinaTsinghua University Press