## Services on Demand

## Article

## Indicators

## Related links

- Cited by Google
- Similars in Google

## Share

## South African Journal of Industrial Engineering

##
*On-line version* ISSN 2224-7890

*Print version* ISSN 1012-277X

### S. Afr. J. Ind. Eng. vol.26 n.3 Pretoria Nov. 2015

#### http://dx.doi.org/10.7166/26-3-802

**GENERAL ARTICLES**

**Model and algorithms of the fuzzy three-dimensional axial assignment problem with an additional constraint**

**C.-J. Lin ^{I}; K.-T. Ma^{II}**

^{I}Department of Industrial Engineering and Management Ta Hwa University of Science and Technology, Taiwan lcj@tust.edu.tw

^{II}Department of Industrial Engineering and Engineering Management National Tsing Hua University, Taiwan d9734802@oz.nthu.edu.tw

**ABSTRACT**

This study constructs a practical fuzzy three-dimensional axial assignment model, and proposes two efficient algorithms to solve the model. In our case, the model is applied to team performance management in a company to promote the performance of all members in a team. Two algorithms, namely the index-based branch and bound (B&B) algorithm and the *f-g *trade-off algorithm, which is a hybrid of the trade-off and B&B concepts, are proposed. A numerical example is presented to illustrate these two algorithms. The computational results show that the proposed algorithms are sufficiently efficient and accurate. Two special cases are also discussed.

**OPSOMMING**

'n Praktiese, wasige, driedimensionele aksialetoekenningsmodel word voorgehou en twee doeltreffende algoritmes om die model op te los word voorgestel. Die model word toegepas op die spanprestasiebestuur in 'n maatskappy met die doel om die spanlede se vertoning te verbeter. Twee algoritmes, naamlik die indeksgebaseerde vertak-en-begrens algoritme en die *f-g *kompromie algoritme, word voorgestel. 'n Numeriese voorbeeld word ter illustrasie van die twee algoritmes voorgehou. Die resultate toon dat die voorgestelde algoritmes doeltreffend en akkuraat is. Twee spesiale gevalle word ook bespreek.

**1 INTRODUCTION**

The three-dimensional (3D) axial assignment problem is an extension of the linear assignment problem - i.e., the 2D assignment problem. The 3D axial assignment problem is actually NP-hard [1]. It is usually stated as follows: assign *n *workers to *n *jobs on *n *machines to minimise the total cost [2]. There should be one-to-one correspondences between workers and machines, workers and jobs, and machines and jobs. These correspondences make up a cubic assignment. The model can be formulated as the following 0-1 integer programming problem [1-5]:

where *c _{ijk} *is the cost of assigning worker

*i*to job

*j*on machine

*k,*and all

*c*

_{ijk}*s construct a cubic matrix. If*

*x*

*= 1, then worker*

_{ijk}*i*is assigned to do job

*j*on machine k. If

*x*= 0, then worker

_{ijk}*i*is either not assigned to job

*j*or not assigned to machine

*k.*In a deterministic environment, the 3D axial assignment problem is important in operations research. Applications of the 3D axial assignment problem arise in many areas, such as scheduling capital investments, addressing a rolling mill, military troop assignment, and satellite coverage optimisation [1,6]. Algorithms for the 3D axial assignment have been developed. However, the coefficient matrix of Equation 1.1 is not totally unimodular to ensure an integer solution [7]; thus it cannot be treated as efficiently at continuous problem-solving for its optimum as in a 2D assignment problem. Hence, algorithms to search for the exact solution of a 3D axial are mostly implicit enumeration methods [8], such as the branch and bound (B&B) method [3,9]. Some papers have developed good bounds for the B&B method, especially the lower bound for minimising the problem [6,8,10]. Hung and Lim [11] proposed a local search genetic algorithm-based method to solve the 3D axial assignment problem, but they could not guarantee that an optimal solution could be obtained.

Although all *c*_{ijk}s are deterministic numbers in Equation 1.1, the cost of many real-world applications may not be deterministic numbers. Many researchers have begun to study the 2D assignment problem in a fuzzy environment. These situations include various fuzzy assignment problem (FAP) models and algorithms [12-15], sensitivity analysis of FAP [16], quadratic FAPs [17], multi-criteria FAPs [18-20], traffic FAPs [21], two-objective FAPs [20,22], group decision-making [23], and multi-job FAPs [24]. Solution procedures of these 2D FAPs all take advantage of the totally unimodular property for efficiency. However, no paper discusses the modelling and methods for a fuzzy 3D axial assignment problem. Therefore, this paper presents a fuzzy 3D axial assignment model, and proposes two B&B-based algorithms to solve the model. The next section describes the related fuzzy theory and construction of the fuzzy 3D axial assignment model. Section 3 presents two algorithms for solving this model and a numerical example. Section 4 describes computational results, and Section 5 offers a conclusion to the paper.

**2 MODEL CONSTRUCTION**

**2.1 Basic concept of related fuzzy theory**

Fuzzy theory is useful for accessing many real situations of human life. In this section, we first recall several related definitions of fuzzy theory, which are used in this paper.

**Definition 2.1.1. **A crisp set *A *is a proper subset of a universe of discourse *X *such that ** A **⊆

**; the membership function of**

*X**A*can be defined as [25]:

**Definition 2.1.2. **A fuzzy set *Ã *is a subset of a universe of discourse X, which is characterised by a membership function *u~ _{A}*(x),

*x*

*e*

*X*, which maps

*x*to a value within the range of [0,1] [26].

According to **Definitions 2.1.1. **and **2.1.2., ***x *completely belongs to ** A **if the value of the membership function

*u*(x) = 1;

_{Ã}*x*does not belong to

**if**

*A**u*x) = 0.

_{Ã}(**Definition 2.1.3. **The λ-cut of *Ã *of the universe of discourse X is defined as .

Therefore, the λ-cut is a subset of the universe of discourse X, where *u _{Ã}*(x)

__>__

*λ*

*.*

**Definition 2.1.4. **The largest value of * λ*of the λ-cut of

**, which is non-empty, is called the height of**

*Ã**Ã*, denoted as

*h(*

**) [26].**

*A***Definition 2.1.5. **A fuzzy set *Ã* of the universe of discourse X is normal if h( ** A **)=1; otherwise, a fuzzy set

**of the universe of discourse X is subnormal [15,26,27].**

*A***Definition 2.1.6. **A fuzzy set *Ã*of a universe of discourse *X *is convex if and only if .

**Definition 2.1.7. **A fuzzy set ** Ã**of a universe of discourse

*X*is a fuzzy number or fuzzy interval if

**is both normal and convex [26].**

*A***Definition 2.1.8. **A fuzzy set ** Ã**of a universe of discourse

*X*is a subnormal fuzzy number or subnormal fuzzy interval if

**is both subnormal and convex [27].**

*Ã***Definition 2.1.9. **A subnormal fuzzy interval [15,27] is a subnormal fuzzy subset of the universe of discourse X, which is both subnormal and convex.

Accordingly, a fuzzy number and fuzzy interval are a special case of subnormal fuzzy number and subnormal fuzzy interval, such that the maximal height of the membership function *u _{Ã}*(x) is 1, respectively.

**2.2 Model formulation**

Suppose a project team consists of *n *workers and a manager. The workers and the manager are responsible for performing *n *jobs and for minimising the total costs, respectively. Each worker is assigned to one and only one job to perform it. In many cases, the job quality achieved depends on the assigned worker's skill, the machine used, and the input-cost of the job. Generally, a higher job quality indicates greater job cost. The job quality achieved is defined as the performance of the worker. Hence, workers desire high-input job costs and a suitable machine to perform the assigned job well. Thus the cost of worker *i *performing job *j *on machine k, which is denoted as *c _{m}*, is positively related to the performance of worker .

Let *q**. _{Ijk}*(0

__<__7<q

_{ijk}

__<__1) be the highest possible quality of job

*j*performed by worker

*i*on machine k. Generally,

*q..*should be between 0.5 and 1.0. The membership function of is then defined as the linear monotone increasing function in Equation 2.1 and Figure 1. In the function, α

_{k}*is the minimal cost for worker*

_{ijk}*i*to perform job

*j*on machine

*k,*and β

*is the minimal cost required to reach the maximal job quality*

_{ijk}*q*When the job input cost is between and β

_{ijk}*, the job cost is linearly and positively related to the job quality achieved. However, any expense exceeding β*

_{ijk}*., is wasteful because quality can no longer be enhanced. Without loss of generality, it is assumed that 0 < α*

_{ijk}*«β*

_{ijk}*. Condition*

_{ijk}*x*= 1 is added to Equation 2.1 because there is no real expense if

*x*= 0 in any feasible solution.

_{ijk}

Notation «α_{ijk} β_{ijk} > is used to denote . The cubic matrix is shown as

Additionally, all α_{ijk}*S* form cubic matrix *[**α** _{ijk}, *all β

_{ijk},ε form cubic matrix [β

_{ijk}], and all q

_{ijk}s form cubic matrix [q

_{ijk}]. Define the cubic matrix [γ

_{ijk}] as

where *γ** _{ijk}, *is the slope of the membership function.

The manager's duty is to keep the total cost (c_{T}) as low as possible. Numbers *a *and *b *are assumed to be constants chosen subjectively by the company. When the total cost is lower than α, the manager's performance is 1, when the total cost is greater than b, the manager's performance is 0, and when the total cost is between *a *and *b, *the total cost is linearly and negatively related to the manager's performance. Thus, the membership function of total cost is defined as the linear monotonically decreasing function in Equation 2.3), which is shown in Figure 2. Notation <a, *b> *denotes the fuzzy interval It is suggested that *a *should be a number less than or equal to the minimum of Equation 1.1, with *a _{ijkk}*as its parameter, and

*b*should be a number larger than or equal to the maximum of Equation 1.1, with B

*as its parameter.*

_{ijk}

Because the performance of worker *i *(π_{i}) is the job quality achieved, then

The performance of the manager (*π*_{m}) is

Usually, the company takes the lowest performance of the team members, composed of all workers and the manager, as the team performance for promoting overall performance. Thus the performance of every team member must be equally emphasised. To maximise the team performance is then to maximise the minimal performance of all team members. Therefore, the study uses Bellman-Zadeh's criterion [28], which maximises the minimum of all membership functions corresponding to that solution. The objective function is then defined as

or

where *x _{ijk}*is an element of feasible solution x of Equation 1.1. The fuzzy 3D axial assignment problem model is represented as follows:

Using the membership functions from Equation 2.1 and Equation 2.3, Equation 2.7 can be represented as the following equivalent model:

where * *is the λ-cut of and is the corresponding total cost . In

Equation 2.8, the and λ are all decision variables. Hence, Equation 2.8 is a mixed nonlinear programming model.

In our case, Equation 2.8 is applied to team performance in human resource management that equally emphasises the performance of the manager and workers. However, Equation 2.8 can also be used for situations that equally emphasise the quality and quantity, time and cost, humanity and technology, and soon.

**2.3 Model transformation**

Suppose ω = {( /, *j, *k) | *x _{ijk}.. *= 1} is the corresponding cubic assignment of the feasible solution

**x**of (1.1) and the discussion is confined based on ω. Then, Equation 2.8 can be simplified as

**Theorem 2.1. **Let λ_{x} be the optimal objective function value of (2.9). Then,

** Proof. **Rewriting Equation 2.9 into a linear programming model results in

The dual problem of Equation 2.10 is then

Obviously, Equation 2.10 produces > 0 for all (i, j, *k)e *ω. Using the complementary slackness theorem yields the following:

Based on Equation 2.12 and Equation 2.2, Equation 2.11 can be simplified as follows:

The dual problem of Equation 2.13 is

Rewriting Equation 2.14 as

completes the proof.

Because λ_{χ} is the maximal value of Equation 2.9 by giving the feasible solution **x **of Equation 1.1, the maximum of λ_{χ} must be the optimal solution of Equation 2.8. Equation 2.8 is rewritten as

Obviously, Equation 2.16 is a special multi-fractional max-min 0-1 programming problem.

**3 SOLUTION ALGORITHM**

This section presents two efficient algorithms for solving the exact solution of Equation 2.16. The constraint of Equation 2.16 is the same as that of Equation 1.1 - i.e., the lack of a totally unimodular property. The methods in the literature for solving FAP are difficult to implement for Equation 2.16. Hence, implicit enumeration methods are employed. One of the proposed algorithms is modified from the index-based B&B algorithm [8-9], which quickly solves Equation 1.1 [9]. The other is a hybrid of the efficient trade-off concept [29] and index-based B&B algorithm. A numerical example is shown to clarify these algorithms.

**3.1 The index-based B&B algorithm**

This algorithm performs implicit enumeration of all feasible solutions for solving Equation 2.16. Figure 3 shows part of the tree graph of a 4x4x4 problem used to illustrate the algorithm. The proposed algorithm proceeds from tree nodes from top to bottom and from left to right until all nodes are checked or eliminated by a fathoming test. For each node (*.*, *j, *k) residing at level /, a *path *leads to it from the root node. The variables in the path are fixed at 1, and the other variables at level *t *(0< *t< /) *equal 0. Thus, each node is associated with a set of feasible solutions. For example, node (2, 3, 2) at level 2 in Figure 3 corresponds to the feasible solutions *x *= *x *= 1; other variables at level 1 and 2 are 0, while other variables are undetermined.

Consider the following notations:

Assume that the algorithm has currently reached node (p, *j *, *k *) (p>0) and that nodes (1, *j , *k_{j}), (2, *j _{2}, *k

_{2}),...,(p,

*j*,

*k*) compose the path - Path (p) - from the root node to the current node. In addition, let

**Theorem 3.1. **Function UPV provides an upper bound of the optimal objective value of Equation 2.16 corresponding to Path (p).

**Proof. **Using Equations 3.3, 3.5, 3.6, and 3.8 yields

Suppose nodes (p+1, *j *, *k *), (p+2, *j *, *k *+_{2}),..., and (n, *j, k _{n}) *compose the optimal solution of Equation 2.16 corresponding to Path (p) and call it Path (n). Then,

Using Equation 3.1 produces . Thus,

Similarly, it can be verified that

Using Equation 3.12 and Equation 3.13 produces

and completes the proof.

The following theorem is derived from Equation 2.16.

**Theorem 3.2 **If (i, j, *k) *is a node of Path (p), then *q _{ijk}. *gives an upper bound of the optimal objective function value of Equation 2.16 corresponding to Path (p).

Let z* represent the largest known objective function value, which is associated with a feasible solution. Using Theorems 3.1 and 3.2, the solution procedure of the index-based B&B algorithm is described as follows:

**Step 1: **Find ** **and ,- for

*i*=1, 2, n.

**Step 2: **Compute SLA(i) and *SLG(i) *for *i *=1, 2, n+1. Let TA(0)= SLA(1) and TG(0)= SLG(1).

**Step 3: **Use *x _{111}*=

*x*= .=

_{222}*x*= 1 as the initial feasible solution and its objective function value as z*.

_{nnn}**Step 4: **Let p= 1 and proceed downward to node (1, 1, 1) from the root node.

**Step 5: **Suppose node (p, *j *, *k *) is reached.

(1) If *q*_{p}j_{pkp}*< *z* or UPV (p)< z*, then node (p, *j *, *k *) is fathomed and moves toward the immediate right node of (p, j_{p}, k_{p}) at level p. If node (p, j_{p}, k_{p}) is the last node of level p, then move upward to the first node right of (p-1, *j, *k) at level p-1 and update *p <- *p-1.

(2) If *q*_{p}j_{pkp}*> *z*, UPV (p)> z*, and *n, *then move downward from (p, j_{p}, k_{p}) to the first node at level p+1 and update *p - *p+1.

(3) If *q*_{p}j_{pkp}*> *z*, *UPV *(p)> z*, and p= n, then replace z* with UPV (p). Move upward

to the first node right of (p-1, *j *, *k *) at level p-1 and update *p <- p-1. *

**Step 6: **If p= 0, then stop; otherwise, proceed to Step 5.

The main concept of the index-based B&B algorithm is as follows. Whenever the algorithm reaches a new node *(p, j _{p}, k_{p}), *the two upper bounds (UPV

*(p)*and

*q*) are calculated for that path before proceeding forward to any branches leaving from node

_{pjpkp}*(p, j*If both UPV (p) and

_{p}, k_{p}).*q*

*are larger than z*, then a better feasible solution may be found by proceeding from node (p,*

_{p}j_{pkp}*j*,

*k*). However, if UPV (p) or

*q*are less than or equal to z*, no better feasible solution can be obtained by proceeding forward from node

_{pjpkp}*(p, j*That is,

_{p}, k_{p}).*(p, j*is fathomed. The algorithm locates the exact optimal solution.

_{p}, k_{p})**3.2 The f-g trade-off algorithm**

Because the performances of team members are equally emphasised to maximise the team performance, the *f-g *trade-off algorithm is proposed. The *f-g *trade-off algorithm balances the maximum of the manager's performance and minimum of the workers' performances to maximise team performance. In the algorithm, some procedures are repeated to solve updated problems with different parameters. In each solution, the maximum of the manager's performance decreases (or remains unchanged), and the minimum of the workers' performances simultaneously increases until the maximal team performance is reached.

To obtain the maximal manager's performance, Equation 3.14 is solved.

where *t *is an index of iteration and *t *= 0 indicates the original problem parameter. Equation 3.14 is a variant of the 3D axial fractional assignment problem. To solve Equation 3.14 efficiently when *t *= 0, the concept of the fraction cost penalty method (FCPM) [30] is used to obtain a good initial feasible assignment. Cubic matrix defined as the cost coefficients of Equation 1.1, but its objective function is changed to the maximal case. An approach of the cubic FCPM is applied to [ψ ]. In [ψ^], the penalty for each row, column, and height was calculated, defined as the difference between the two largest elements of the row, column, and height. In the row, column, or height with the largest penalty, the largest ψ _{ijk}* *must be selected and its row, column, and height crossedout. The procedure is repeated to obtain a feasible assignment and its corresponding objective function value. The obtained objective function value and assignment are z* and ω*, respectively. Other than in the cubic FCPM, z* and ω* represent the largest known objective function value and associated assignment, respectively.

To optimise Equation 3.14, the revised branch and bound (R-B&B) algorithm is defined. It is used repeatedly in the f-g trade-off algorithm. The R-B&B algorithm is reduced from the index-based B&B algorithm.

__The R-B&B algonthm:__

All steps of the R-B&B algorithm are the same as those in the index-based B&B algorithm, except that the condition *q** _{p}j_{pkp}*< z* is deleted from (1) of Step 5 and

*q*

*> z* is deleted from (2) and (3) of Step 5.*

_{p}j_{pkp}Whenever the R-B&B algorithm is performed, it produces the optimal objective function value f^{(t)} of Equation 3.14 and its associated assignment a)^{(t)}. Equation 3.15 is then used to compute g^{(t)}:

where *g ^{(t)}*is the minimum of all workers' performance at iteration

*t.*The largest known performance of the team z* can then be updated by using

*f*and

^{(t)}*g*

^{(t)}.To promote z*, the cubic matrix [α* _{ijk}*] in (3.14) is updated; the updated rule of [α ] is

(1) if *q _{ijk}*

__>__z*, α

_{ijk}remains unchanged.

(2) If *q**.j _{k}*

__<__z*,

*α*

*is replaced with M (M-->oo) - i.e., fathom out that is useless.*

The rule means that f^{(t)} decreases (or remains unchanged) and g^{(t)} increases at each iteration for z* to reach the maximum. Therefore, the flowchart in Figure 4 is used to depict the proposed *f-g *tradeoff algorithm.

**3.3 Numerical example**

The following numerical 3x3x3 fuzzy 3D axial assignment problem illustrates the two proposed algorithms. Table 1 shows the job costs and quality.

According to Equation 2.16, the numerical example can be formulated as

Using Equation 2.2, α* _{ijk}* and β

*are transformed into γ*

_{ijk}*as shown in Table 2,*

_{ijk}*a*= 309, and

*b*= 612. The algorithms are then used to solve Equation 3.16.

(1) The index-based B&B algorithm __Iteration ____0:__

**Step 1: **Find *α** _{1}*= 101, α

_{2}= 104, α

_{3}= 104, Yj = 28.17, γ

_{2}=24, and γ

_{3}=26.03.

**Step 2: **Determine SM(1) = 309, SM(2) = 208, SM(3) = 104, SM(4) = 0, SLG(1) = 78.2,* SLG(2) *= 50.03, SLG(3) = 26.03, SLG (4) = 0, *TA(0) *= 309, and TG(0)= 78.2.

**Step 3: **Use *x _{111}*=

*x*=

_{222}*x*= 1 as the initial feasible solution and obtain z*= 0.60.

_{333}**Step 4: **Let *p *= 1 and proceed downward to node (1, 1, 1) from the root node.

**Step 5: **Identify TA(1) = 309 and TG(1) = 95.34. Because q = 0.64 > 0.6 = z* and UPV(1) = 0.76 > 0.6 = z*, proceed downward to node (2, 2, 2) and update *p *= 2.

**Step 6: **Let *p *= 2 and proceed to Step 5.

__Iteration ____1:__

**Step 5: **Because q_{222}= 0.6 __<__ z*, node (2, 2, 2) is fathomed. Proceed toward node (2, 2, 3).

**Step 6: **Let *p *= 2 and proceed to Step 5.

__Iteration ____2:__

**Step 5: **Determine q_{223} = 0.69 > 0.6 = z*, *TA(2) *= 309, and TG(2) = 113.37. Because UPV(2)= 0.728 > 0.6= z*, proceed downward to node (3, 3, 2) and update *p *= 3.

**Step 6: **Let *p *= 3 and proceed to Step 5.

__Iteration ____3:__

**Step 5: **Determine q_{332}= 0.69 > 0.6= z*, TA(3) = 358, and TG(3) = 118.04. Because UPV(3)= 0.604> 0.6= z* and p= 3= n, replace z*= 0.604 and move upward to node (2, 3, 2) and to the first node right of (2, 2, 3) of level 2. Update *p *= 2.

**Step 6: **Let *p *= 2 and proceed to Step 5.

The remaining procedures repeat these iterations and obtain the optimal solution *x _{111}*=

*x*=

_{223}*x*= 1 and optimal value z*= 0.604.

_{332}(II) The f-g trade-off algorithm

__When t = 0:__

Equation 3.14 produces Equation 3.17.

After using the cubic FCPM and R-B&B program to solve Equation 3.17, ω^{(0)} = (x , x_{222}, x_{333}), f^{(0)}=0.662, and g^{(0)} = 0.6. Let ω*= ω^{(0)} and z*= Min (0.662, 0.6)= 0.6 = g^{(0)}.

__When t = ____1:__

Update by replacing α_{131}, α_{132}, α_{222}, α_{231}, α_{232}, and α_{313} with M and obtain

Solving Equation 3.17 again by using [ α] ] produces ω^{(1)} = (x_{111}, x_{223}, x_{332}), f^{(1)}=0.604, and g^{(1)}= 0.64. Because f^{(1)} = 0.604 > 0.6= z*, update ω*= ω^{(1)} and z*= Min (f^{(1)}, g^{(1)})= 0.604.

__When t = ____2:__

Update ] by replacing α_{131} , α_{132} , α_{222} , α_{231} , α_{313} , and α with M and obtain

Solving Equation 3.17 again by using produces ω^{(2)}= (x_{111}, x_{223}, x_{332}), f^{(2)}= 0.604, and g^{(2)} = 0.64. Because f^{(2)}= 0.604 = z*, the algorithm stops. The current ω*= (x_{111}, x_{223}, x_{332}) and z*= 0.604 are the optimal solution and optimal objective function value, respectively.

**4 COMPUTATIONAL RESULTS AND DISCUSSION**

The computational test verifies the efficiency of the proposed algorithms. Both proposed algorithms were coded in Visual Basic 6.0 and compared with the LINGO 11 package by using the global optimal function. They were all run on an Intel Pentium III-500 CPU.

Problem sizes from 3x3x3 to 10x10x10 were tested. For each problem of NxNxN size, uniform distribution numbers between 10 and 10+20N were randomly generated as α_{ijk} s, and uniform distribution numbers between 6N to 10N were randomly generated as the difference between β_{ijk} s and α_{ijk} s. Uniform distribution numbers between 0.6 and 1.0 were* *generated as *q _{ijk}s. *For each problem size, 30 problems were generated. Table 3 shows the average computational times (in seconds) of the two proposed algorithms and the LINGO package. The results show that both proposed algorithms have shorter computational times than the LINGO package.

In general cases, the manager or any worker may be the person with the lowest performance and may determine the team's performance. Here, two special cases are presented. In the first case, the company is assumed to support the team with sufficient funds, e.g., 10 times of b. The computational results of the three approaches are shown in Table 4. The results show that the index-based B&B algorithm becomes more efficient. The time taken by the LINGO package also reduces. The *f-g *trade-off algorithm, on the contrary, becomes inefficient. This is because the optimal values f^{(t)} of Equation 3.14 approach 1, such that the fathom speed of a* _{ijk}* decreases. To discuss this further theoretically, b→∞ moves the membership function in Figure 2 close to the horizontal line μ

*=1 and*

_{T}(c_{T})

i.e., the manager's performance approaches 1. The manager's performance is greater than the maximal of q_{ijk}. Thus the manager never determines the team performance; team performance is decided by a worker. The workers also obtain sufficient expense and reach the highest quality (performance). Therefore, Equation 2.16 is simplified to

Equation 4.1 is a variant of the 3D axial bottleneck assignment problem. Hence, a proper algorithm for the 3D axial bottleneck assignment problem will improve the efficiency if b-cc is known in advance.

In the second case, *b *is small; that is, the team is underfinanced. For example, define the manager's performance as 0 when the total cost is greater than (a+b)/2. The computational results of the three approaches are shown in Table 5. The results show that the *f-g *tradeoff algorithm is the most efficient. Theoretically, when

where Z_{[}_{α]*} is the optimal value of the crisp 3D axial assignment problem (min) using *[a _{ijk}*] as the cost coefficient, the

*b*is too small for the manager's performance to reach the minimum of q

_{ijk}. Team performance is then determined mainly by the manager, despite the max-min objective function equalising the performance of the

*n*workers and the manager. Thus Equation 2.16 is simplified to Equation 3.14 with t=0. Hence, a proper specific algorithm for the 3D axial fractional assignment problem will improve the efficiency.

**5 CONCLUSIONS**

This paper first formulates a fuzzy 3D axial assignment problem model to overcome the actual uncertain environment, and then focuses on developing algorithms to obtain the exact solutions of the proposed model efficiently. The initially-constructed fuzzy 3D axial assignment problem is a mixed nonlinear programming problem. For simplicity, it is transformed into a special multi-fractional max-min 0-1 programming problem. An index-based B&B algorithm and the *f-g *trade-off algorithm are developed. Computational results show that the proposed methods are both more efficient than LINGO in every test. In addition, two special cases are discussed. One is when the company amply supports the team with funds, while the other is when the team is underfunded. With respect to computational time, the index-based B&B algorithm is better for the former case. However, the *f-g *trade-off algorithm is better for the latter case. These two special cases are further considered from a theoretical view. For the case with unbounded funds, the model is reduced to a 3D axial bottleneck assignment problem; for the case with rare funds, the model is reduced to a 3D axial fractional assignment problem. Hence, the developments of algorithms for these two problems would be interesting and practical directions for future research. Moreover, a sensitivity analysis procedure for identifying the largest sensitivity range - one that maintains the current optimal assignment - is worth studying.

**REFERENCES**

[1] **Burkard, R., Dell'Amico, M. & Martello, S. **2012. *Assignment problems. *1^{st} edition, SIAM. [ Links ]

[2] **Pentico D.W. **2007. Assignment problems: A golden anniversary survey. *European Journal of Operational Research, *176(2), pp. 774-793. [ Links ]

[3] **Balas, E. & Saltzman, M.J. **1991. An algorithm for the three-index assignment problem. *Operations Research, *39(1), pp. 150-161. [ Links ]

[4] **Poore, A.B. & Robertson, A.J. **1997. A new Lagrangian relaxation based algorithm for a class of multidimensional assignment problems. *Computational Optimization and Applications, *8(2), pp. 129-150. [ Links ]

[5] **Pierskalla, W.P. **1968. The multidimensional assignment problem. *Operations Research, *16(2), pp. 422-431. [ Links ]

[6] **Kim, B.J., Hightower, W.L., Hahn, P.M., Zhu, Y.R. & Sun, L. **2010. Lower bounds of the axial three-index assignment problem, *European Journal of Operational Research, *202(3), pp. 654-668. [ Links ]

[7] **Adlakha V. & Arsham, H. **1998. Managing cost uncertainties in transportation and assignment problems. *Journal of Applied Mathematics & Decision Sciences, *2(1), pp. 65-104. [ Links ]

[8] **Burkard, R.E. & Rudolf, R. **1993. Computational investigations on 3-dimensional axial assignment problems. *Belgian Journal of Operations Research, Statistics & Computer Science, *32(1-2), pp. 85-98. [ Links ]

[9] **Pasiliao, E.L., Pardalos, P.M. & Pitsoulis, L.S. **2005. Branch and bound algorithms for the multidimensional assignment problem. *Optimization Methods & Software, *20(1), pp. 127-143. [ Links ]

[10] **Sergeev, S.I. **2006. The three-dimensional assignment and partition problems - New lower bounds. *Automation & Remote Control, *67(2), pp. 242-250. [ Links ]

[11] **Huang, G. & Lim, A. **2006. A hybrid genetic algorithm for the three-index assignment problem, *European Journal of Operational Research, *172(1), pp. 249-257. [ Links ]

[12] **Lia, F., Xu, L., Jina, D.C. & Wang, H. **2012. Study on solution models and methods for the fuzzy assignment problems. *Expert Systems with Applications, *39(12), pp. 11276-11283. [ Links ]

[13] **Lin, C.J. & Wen, U.P. **2004. A labeling algorithm for the fuzzy assignment problem. *Fuzzy Sets and Systems, *142(3), pp. 373-391. [ Links ]

[14] **Mukherjee, S. & Basum, K. **2012. Solution of a class of Intuitionistic fuzzy assignment problem by using similarity measures. *Knowledge-Based Systems, *27, pp. 170-179. [ Links ]

[15] **Lin, C.J. **2013. Assignment problem for team performance promotion under fuzzy environment. *Mathematical Problems in Engineering, *2013, pp. 1-10. [ Links ]

[16] **Lin, C.J., Wen, U.P. & Lin, P.Y. **2011. Advanced sensitivity analysis of the fuzzy assignment problem. *Applied Soft Computing, *11(8), pp. 5341-5349. [ Links ]

[17] **Liu, L.Z. & Li, Y.Z. **2006. The fuzzy quadratic assignment problem with penalty: New models and genetic algorithm. *Applied Mathematics and Computation, *174(2), pp. 1229-1244. [ Links ]

[18] **Belacela, N. & Boulasselb, M.R. **2001. Multicriteria fuzzy assignment method: A useful tool to assist medical diagnosis. *Artificial Intelligence Medicine, *21(1-3), pp. 201-207. [ Links ]

[19] **Dubois, D. & Fortemps, P. **1999. Computing improved optimal solutions to max-min flexible constraint satisfaction problems. *European Journal of Operational Research, ***118**(1), pp. 95-126. [ Links ]

[20] **Huang, D.K., Chiu, H.N., Yeh, R.H. & Chang, J.H. **2009. A fuzzy multi-criteria decision making approach for solving a bi-objective personnel assignment problem. *Computers & Industrial Engineering, *56(1), pp. 1-10. [ Links ]

[21] **Ridwan, M. **2004. Fuzzy preference based traffic assignment problem. *Transportation Research Part C: Emerging Technologies, *12(3-4), pp. 209-233. [ Links ]

[22] **Feng, Y. & Yang, L. **2006. A two-objective fuzzy k-cardinality assignment problem. *Journal of Computational & Applied Mathematics, *197(1), pp. 233-244. [ Links ]

[23] **Bashiri, M. & Badri, H. **2011. A group decision making procedure for fuzzy interactive linear assignment programming. *Expert Systems with Applications, *38(5), pp. 5561-5568. [ Links ]

[24] **Liu, L. & Gao, X. **2009. Fuzzy weighted equilibrium multi-job assignment problem and genetic algorithm. *Applied Mathematical Modelling, *33(10), pp. 3926-3935. [ Links ]

[25] **Mutingi, M. & Mbohwa, C. **2014, A fuzzy-based particle swarm optimisation approach for task assignment in home healthcare. *The South Africa Journal of Industrial Engineering, *25(3), pp. 84-95. [ Links ]

[26] **Klir, G.J., Clair, U.S. & Yuan, B. **1997, *Fuzzy set theory: Foundations and applications. *1^{st }edition, Prentice Hall PTR. [ Links ]

[27] **Baruah, H.K. **2011, Theory of fuzzy sets: The case of subnormality. *International Journal of Energy, Information and Communications, ***2**(3), pp. 1-8. [ Links ]

[28] Bellman, R.E. & Zadeh, L.A. 1970. Decision-making in a fuzzy environment. *Management Science, ***17B**(4), pp. 141-164. [ Links ]

[29] **Geetha, S. & Vartak, M.N. **1989. Time-cost tradeoff in a three dimensional assignment problem. *European Journal of Operational Research, *38(2), pp. 255-258. [ Links ]

[30] **Joshi, V.D. & Gupta, N. **2009. A heuristic for obtaining a better initial solution for the linear fractional transportation problem. *Bulletin of Australian Society of Operations Research, *28(4), pp. 45-49. [ Links ]

* Corresponding author