Issue 
Mechanics & Industry
Volume 18, Number 6, 2017



Article Number  606  
Number of page(s)  12  
DOI  https://doi.org/10.1051/meca/2017004  
Published online  21 March 2018 
Research Article
Multiobjective robust design optimization of the mechanism in a sewing machine
^{1}
LGM, ENIM, University of Monastir,
Monastir, Tunisia
^{2}
Department of Mechanical Engineering, American University of Sharjah,
Sharjah, UAE
^{*} email: najlawibilelali@gmail.com
Received:
13
June
2016
Accepted:
4
January
2017
This work deals with the multiobjective robust design optimization of a needlebarandthreadtakeuplever (NBTTL) mechanism used in sewing machines. A multiobjective imperialist competitive algorithm combined with the Monte Carlo method is developed and used for the robust multiobjective optimization of the NBTTL mechanism. This robust optimization minimizes simultaneously the needle jerk, the maximum angular velocity ratio and their standard deviations when the design parameters are known with a given uncertainty. By comparison with Juki 8700, the obtained results showed that the developed algorithm can give a robust design of a sewing machine with the same performance as the those obtained by deterministic methods but less sensitive to the design parameters uncertainties.
Key words: Robust design / sewing machine / uncertainty / imperialist competitive algorithm / Monte Carlo simulation
© AFM, EDP Sciences 2017
1 Introduction
Many researchers focused on optimizing the designs of sewing machines in order to improve their mechanical performances. Some works, treated the effect of the design parameters (DPs) of the needlebarandthreadtakeuplever (NBTTL) mechanism on their performances. Payvandy and Ebrahimi [1] used a monoobjective function with penalty weighting factors to minimize the coupler tracking error and the Jerk of the thread takeup lever. Ebrahimi et al. [2] presented an optimization problem based on the imperialistic competitive algorithm (ICA) for optimizing the needle velocity and the variation of needle acceleration in a link drive mechanism of a sewing machine. Komarec et al. [3] develop a monoobjective optimization of the needle transfer mechanism by means of the software Pro/Engineer Wildfire 4 with the aim to obtain a reduction of needle wear.
It is a common practice in sewing machine design optimization to consider the nominal value of DPs as input variables. However, the sewing machine DP values are usually known with a certain uncertainty. This uncertainty is caused by several factors, i.e., manufacturing errors, geometric errors and material properties. To estimate the effect of DP uncertainty on the performance of a mechanical system, several methods have been used. In particular, the Monte Carlo simulation (MCS) is a popular and simple tool in estimating the effect of uncertainties on the output of a process [4]. Chandrashekhar and Ganguli [5] developed an optimization method based on fuzzy logic for damage detection in the steel beam having material uncertainty. To evaluate the frequency variations of a beam, they performed a probabilistic analysis using the MCS. Cloupet et al. [6] used the MCS method to estimate the wear variability of rolling contact in industrial machinery caused by uncertain friction coefficient. Motevalli et al. [7] applied the MCS to study the impact of water inflow uncertainty on the performance of both single and multireservoir systems.
Other works have studied robust design product where uncertainties of the DPs are considered. Cheng and Li [8] have developed a hybrid differential evolution and sequential quadratic programming method for the robust mechanical structures under uncertain DP. Kalantari et al. [9] have developed a hybrid robust evolutionary algorithm by combining the NSGAII process with a local search method. This strategy is used to optimize composite structures under an uncertain fibre angle and a lamina thickness. Sun et al. [10] developed a robust optimization algorithm based on successive Taguchi approach for the crashworthiness design of hybrid steelaluminum tailor welded thinwalled structures. Ghanmi et al. [11] studied the robust design optimization of complex mechanical structures using the genetic algorithm (GA). The robustness, defined by the ratio of the mean value to the standard deviation, is treated as an objective function.
This paper deals with the multiobjective robust design optimization of the NBTTL mechanism used in sewing machines. A combined algorithm based on the multiobjective imperialist competitive algorithm (MOICA) and the MCS is proposed and used. The obtained results are analyzed and compared to the commercialized sewing machine Juki 8700.
The remainder of this paper is organized as follows. In Section 2, the functioning of the NBTTL mechanism is described. Section 3 presents the objective functions to be considered in the NBTTL mechanism design. Section 4 deals with the deterministic multiobjective design optimization. In Section 5, a combined algorithm is developed and used in the multiobjective robust design optimization. Finally some concluding remarks are presented.
2 The NBTTL mechanism
Figure 1 presents the NBTTL mechanism of a sewing machine. The thread takeup lever mechanism is the fourbar linkage OABC. OEF represents the slidercrank mechanism in which point F denotes the needle. The function of the thread takeup lever, in the stitch formation process, is to ensure appropriate thread feeding. The function of the needle, which is fixed to the needle bar, is to penetrate the fabric. The rotation of the input link (OA) is transformed to the translation of the needle bar through the slidercrank mechanism. The displacement of the needle bar is represented by the distance “S”.
Fig. 1 The NBTTL mechanism. 
3 Kinematics of the NBTTL mechanism
To optimize the NBTTL mechanism, we need to find:

The input–output relation between the crank OA angular rotation and the coupler point, D, displacement;

The needle jerk in point F.
In this section, the modeling of the NBTTL mechanism is presented. The angle θ_{4}, the velocity and the acceleration describing the input of the machine are given by the driving motor.
The needle displacement, with respect to (OX, OY) frame, is given by (Fig. 1): (1)
By taking the time derivative of equation (1), we can find the needle jerk: (2) (3)
In textile engineering, sewing machines have been basically designed for work at high speed. This can lead to high vibration. To decrease the vibration of the sewing machine, it is recommended to reduce the jerk of its needle given by: (4)
NJ defines the first objective function to be minimized.
Generally, the maximum angular velocity ratio (MAVR) of the mechanism is treated geometrically [12]. In this work we will give a closed form solution of the MAVR. Referring to Figure 2, H is the intersection of the two lines (CB) and (OA). It represents the instantaneous rotation centre of link r_{3} with respect to link r_{1}. The closed loop equations of the virtual loop ABH can be given by: (5)
Equation (5) can be projected into their X and Y components. (6)
By using equation (6), we can express λ as follows: (7)
From equations (7) and (8), the MAVR is given by: (9) Minimizing the MAVR yields the minimization of the input to output torque ratio. Minimizing this ratio allows the minimization of accelerations and hence minimizing the inertial forces of the mechanism. Hence, The MAVR represents the second objective function.
Fig. 2 Position of a virtual loop of a fourbar linkage. 
4 Deterministic multiobjective optimization design of the NBTTL mechanism
The multiobjective optimization problem can be defined as minimizing simultaneously the NJ and the MAVR under some constraints. In order to ensure that the input crank of the thread takeup lever can make a complete rotation, the mechanism should satisfy the Grashof law. Moreover, the set of input link angles θ_{4} has to vary in an increasing order (θ_{4}(i) < θ_{4}(i + 1)). The displacement of the needle, determined experimentally, is limited (30 mm ≤ S ≤ 62 mm) and the design variables have to be limited to some bounds due to the practical limitations in the construction of the mechanism. Hence, the multiobjective optimization problem can be formulated as: (10) (11)
The search domain D(DP) of each DP is listed in Table 1.
To solve this optimization problem, the MOICA method [13] will be used.
Range of the design variables.
4.1 The optimization method
The MOICA is a stochastic optimization method developed in [13]. This method is based on the ICA algorithm [13–15]. Figure 3 shows the diagram of the MOICA.
The MOICA begins with a generation of a random initial population. Each element of the population is a country. The country cost is fixed based on the evaluation of the objective functions. The “imperialist” is the most powerful country and the “colonies” are the other countries. The initial empires, composed by several colonies and one imperialist, are formed by the initial population. Each colony starts moving towards the more powerful imperialist (Fig. 4). In this movement, θ and X are numbers generated uniformly X ∼ U(0, β × d) and θ ∼ U(− γ, γ). d is the distance between the imperialist and the colony and β must be greater than 1. γ is a parameter representing the direction deviation [13–15].
In each empire, the random mutation and the crossover operators can improve the colonies quality by adopting more powerful new ones. Therefore, the powerful colony can permute positions with the imperialist. Based on the fast nondominated sorting approach [16], the nondominated imperialists are stored in the Pareto front.
The sum of the imperialist power with a mean power proportion of its colonies, gives the empire power. During the imperialists competition all empires, based on their power, try the possession of the other empires colonies (Fig. 5). The weak empires lose their colonies and collapse. The MOICA algorithm stops when all the weak empires collapse except the powerful one. The MOICA parameters are given in Table 2.
Fig. 3 The MOICA flowchart. 
Fig. 4 Moving colonies toward their relevant imperialist. 
Fig. 5 Imperialistic competition. 
The MOICA parameters.
4.2 Results and discussions
Using the MOICA algorithm, the Pareto front (Fig. 6) can be determined to find the optimal nondominated solutions of the multiobjective optimization problem defined by equation (10) under the constraint presented by equation (11).
In Figure 6, each nondominated solution represents an optimal design vector of the NBTTL mechanism. One can note that the NBTTL having the best MAVR, has the worst NJ and vice versa. Table 3 presents the design vectors of the extreme deterministic optimal solutions SD1 and SD2 taken from the Pareto front. For comparison reason, we give in the same table the design vector of the known sewing machine Juki 8700 [17].
Figure 7 shows the variation of the angular velocity ratio (AVR), for the extreme determinist optimal solutions SD1 and SD2, as a function of the input link angle θ_{4}. For the sake of comparison, we present in the same figure the variation of the AVR obtained by the sewing machine Juki 8700.
From Figure 7 one can note that, compared to commercialized sewing machine Juki 8700, we can reduce the MAVR by about 80% by using the optimal SD1 solution. However, if we consider SD2 that minimizes the NJ we have a reduction of the MAVR by 9.5%.
Figure 8 shows the evolution of the needle Jerk for SD1 and SD2 solutions and Juki 8700. One can note that by considering the SD2 design vector, we decrease the needle jerk down by 70% compared to the one of the Juki 8700.
The obtained results show that the MAVR and the Needle Jerk can be reduced by 80% and 70%, respectively, in the case of the optimized mechanism, compared to Juki 8700 machine. However, in practice the DPs are not known precisely and can present uncertainties around their nominal values. These uncertainties can generate fluctuations of the NBTTL performances evaluated by NJ and MAVR. Therefore, it's interesting to find a robust optimal NBTTL mechanism where the performances have low sensitivity to the DP uncertainties.
Fig. 6 The deterministic optimal solutions. 
Design vector of SD1, SD2 and Juki 8700.
4.3 The NBTTL performance sensitivity
In this section, we investigate the effect of the DP uncertainties on the NBTTL performance. For this goal, each DP is modeled by a normal distribution given by a mean value and a standard deviation. The MCS is performed for every DP where each evaluation consists of a specified number of runs (Fig. 9). For each MCS, all the deterministic design variables are fixed at their nominal values. The uncertainties are selected randomly from their statistical distributions supposed to vary within ±10% of the DPs nominal values [18,19].
With the MCS, we perform n = 10^{4} simulations to determine the mean values of the mechanism performances and their respective standard deviations σ_{NJ}, σ_{MAVR}.
The performances of the obtained deterministic solutions are compared to those of the Juki 8700 machine. The determinist and Juki performances vulnerability to the DPs uncertainties are determined by introducing their corresponding DP under the MCS.
Figures 10 and 11 show the sensitivity of the NBTTL performances (MAVR and NJ) to the design parameters uncertainties for the extreme deterministic optimal solutions (SD1 and SD2). The variations of the same NBTTL performances, obtained by the sewing machine Juki 8700, are depicted in Figure 12.
Table 4 presents the mean and the sensitivity percentage of the NBTTL performances to the DP uncertainties for SD1, SD2 and Juki 8700. The sensitivity percentage of the NBTTL performance u (u = NJ, MAVR) is evaluated as follows: (12)
One can note that the performances of the deterministic optimal solutions (SD1 and SD2) are very sensitive to the DP uncertainties. In fact, a variation of 37% and 35% in the values of NJ and MAVR, respectively, are noted for uncertainties of 10% of the design parameters nominal values. Moreover, the Juki 8700 has a much lower NBTTL performances sensitivity compared with SD1 and SD2 solutions. The performances sensitivity of SD1 and SD2 solutions can reach the double of the Juki 8700 performances sensitivity. Therefore, a robust design of the NBTTL mechanism, that considers the mean value and the standard deviation of the machine performances, is necessary.
In what follows, we try to improve the deterministic optimization in order to obtain robust NBTTL performances to the DPs uncertainties. A multiobjective robust design, where the sensitivity of the NBTTL performances is considered as objective functions, will be conducted.
In order to have a robust design of the NBTTL mechanism, a combined algorithm that uses the MOICA and the MCS (MOICA–MCS) is developed and used.
Fig. 9 The flowchart of MCS. 
Fig. 10 Performance vulnerability of the solution SD1. 
Fig. 11 Performance vulnerability of the solution SD2. 
Fig. 12 Performance vulnerability of the Juki machine. 
Comparison of the performances results of SD1, SD2 and Juki 8700.
5 Robust design of the NBTTL mechanism
A robust design of the NBTTL mechanism should have the minimum values of the objective functions NJ, MAVR and also these performances should have the minimum sensitivity generated by the DP uncertainty. This robustness can be quantified by the standard deviation of the objective function, i.e., high robustness of the solution means low standard deviation, for a given uncertainty of the DP. It is worth mentioning that all the solutions have to satisfy all the constraints imposed by the designer. Consequently the robust optimal design problem can be expressed as follows: (13) (14)
5.1 The MOICAMCS algorithm
In this section, the MOICA–MCS Algorithm (Fig. 13) that combines the MOICA and the MCS is developed and used to obtain robust optimal design of the NBTTL mechanisms. In fact, after the evolution stage, the MOICA sends a generation of countries to the MCS. The MCS performs n simulations for every country and sends back to the MOICA the mean values, , and their standard deviations, σ_{NJ}, σ_{MAVR}. Then, these countries undergo the remaining steps of the MOICA algorithm. The MOICAMCS parameters are the same for MOICA (Tab. 2). The number of MCS n is fixed to 10^{4}.
Fig. 13 The MOICA–MCS algorithm. 
5.2 Results and discussions
The MOICA–MCS algorithm yields the nondominated solutions presented in the Pareto front (Fig. 14). Each point represents a NBTTL mechanism design vector with the corresponding values of the four objective functions. The solutions are more robust if the corresponding performance standard deviation is reduced. To represent graphically the four objective functions, different colors were used to represent the values of the fourth objective function, σ_{NJ}.
The first remark is that and tend to vary in the opposite direction of σ_{NJ} and σ_{MAVR}, respectively. This means that reducing the average of an objective function leads inevitably to increasing its variability [20,21].
Table 5 presents the DPs of four extreme solutions from the Pareto front, as shown in Figure 14. By comparing these solutions, one can notice that the first solution presents the advantage of having the minimum of , the second one has the minimum of , the third one has the minimum of σ_{MAVR} and the fourth one has the minimum of σ_{NJ}.
To analyze the advantage of the robust optimal solutions compared with the deterministic ones, different sections of the robust Pareto front (given by Fig. 14), corresponding to different values of σ_{MAVR}, are presented in Figure 15. The deterministic Pareto front is also presented in the same figure.
For the comparison, the performance vulnerability of the robust solution (SR1) and the determinist one (SD1) are given in the same figure (Fig. 16). One can note that the SR1 solution has a comparable value of to the one obtained by the deterministic solution SD1. However, the performance of SR1 is much less sensitive to the DPs uncertainties (Tab. 6). Moreover, the robust solution SR1 has a slightly higher value for (3% percentage difference), compared to the deterministic solution SD1, but with a much lower sensitivity to the DPs uncertainties.
Compared to the Juki 8700, one can note that SR1 has the lowest NJ with a sensitivity divided by about six (Tab. 6). Moreover, the value of and its sensitivity for SR1 are significantly improved.
From Figure 15, one can note also the coincidence of a robust and a deterministic solution (SR2 + SD2). The performance vulnerability of these two solutions are given in Figure 17 and Table 7. One can remark that for the same values for and of the deterministic solution SD2, the robust solution SR2 presents a much lower variability to DPs uncertainties.
Compared to the Juki 8700, one can note that SR2 solution is much more robust to the DPs uncertainties (Tab. 7).
One can conclude that finding a deterministic optimal solution can lead to solutions highly sensitive to the uncertainty of the design parameters. Therefore, the robust optimization strategy is recommended when seeking optimal solutions where the DPs are known with a certain uncertainty. This procedure was applied to the case of the NBTTL mechanism and we showed that the obtained solutions are highly robust to the DP uncertainties.
Fig. 14 The robust optimal solutions. 
Design vector of the extremely robust solutions.
Fig. 15 Some configurations of robust solutions. 
Fig. 16 Performance vulnerability of the solutions SD1 and SR1. 
Performances comparison of SR1, SD1 and Juki one.
Fig. 17 Performance vulnerability of the solutions SD2 and SR2. 
Performances comparison of SR2, SD2 and Juki one.
6 Conclusion
This paper deals with the multiobjective robust design optimization of the NBTTL mechanism. First, the mechanism is presented and analyzed to calculate the two objective functions, i.e., the jerk and the angular velocity ratio. Then, a deterministic optimization was conducted. It was shown that the optimal deterministic solutions are highly sensitive to the design parameter uncertainties. An algorithm, that combines the multiobjective imperialist competitive algorithm and the Monte Carlo Simulation, is developed and used for the multiobjective robust design optimization of the NBTTL mechanism. In the robust design, the NBTTL performances and their corresponding standard deviations are minimized, simultaneously.
The obtained results showed that the developed algorithm can yield robust designs of the sewing machine with similar performances as the those of the obtained by the deterministic optimization. However, with the robust design approach the performances of the obtained solutions are highly insensitive to the design parameter uncertainties.
References
 P. Payvandy, S. Ebrahimi, Optimization of the thread takeup lever mechanism in lockstitch sewing machine using the imperialistic competitive algorithm, J. Text. Polym. 3 (2015) 12–18 [Google Scholar]
 S. Ebrahimi, P. Payvandy, Optimization of the link drive mechanism in a sewing machine using imperialist competitive algorithm, Int. J. Cloth. Sci. Technol. 26 (2014) 247–260 [CrossRef] [Google Scholar]
 J. Kamerec, J. Beran, M. Lima, M. Machado, J.M. Machado, J.P. Silva, Finding the optimal setting of the sewing needle transfer mechanism using simulation software, in: TRS 2012 41st Textile Research Symposium, Portugal, 2012 [Google Scholar]
 B. Peherstorfer, T. Cui, Y. Marzouk, K. Willcox, Multifidelity importance sampling, Comput. Methods Appl. Mech. Eng. 300 (2016) 490–509 [CrossRef] [Google Scholar]
 M. Chandrashekhar, R. Ganguli, Uncertainty handling in structural damage detection using fuzzy Logic and probabilistic simulation, Mech. Syst. Signal Process. 23 (2009) 384–404 [CrossRef] [Google Scholar]
 S. Cloupet, L. Chevalier, H. Chollet, Probabilistic simulation of wear in loaded rolling contact, Mech. Ind. 8 (6) (2007) 609–630 [Google Scholar]
 M. Motevalli, A. Zadbar, E. Elyasi, M. Jalaal, Using MonteCarlo approach for analysis of quantitative and qualitative operation of reservoirs system with regard to the inflow uncertainty, J. Afr. Earth Sci. 105 (2015) 1–16 [CrossRef] [Google Scholar]
 S. Cheng, M. Li, Robust optimization using hybrid differential evolution and sequential quadratic programming, Eng. Optim. 47 (1) (2015) 87–106 [CrossRef] [Google Scholar]
 M. Kalantari, C. Dong, Ian J. Davies, Multiobjective robust optimization of unidirectional carbon/glass fibre reinforced hybrid composites under flexural loading, Compos. Struct. 138 (2016) 264–275 [CrossRef] [Google Scholar]
 G. Sun, J. Fang, X. Tian, G. Li, Q. Li, Discrete robust optimization algorithm based on Taguchi method for structural crashworthiness design, Expert Syst. Appl. 42 (2015) 4482–4492 [CrossRef] [Google Scholar]
 S. Ghanmi, M. Guedri, M. L. Bouazizi, N. Bouhaddi, Robust multilevel and multiobjective optimization of complex mechanical structures, Mech. Ind. 11 (2010) 393–400 [Google Scholar]
 M. Khorshidi, M. Soheilypour, M. Peyro, A. Atai, M. Shariat Panahi, Optimal design of fourbar mechanisms using a hybrid multiobjective GA with adaptive local search, Mech. Mach. Theor. 46 (2011) 1453–1465 [CrossRef] [Google Scholar]
 B. Najlawi, M. Nejlaoui, Z. Affi, L. Romdhane, An improved imperialist competitive algorithm for multiobjective optimization, Eng. Optim. 48 (11) (2016) 1823–1844 [CrossRef] [Google Scholar]
 S. Hosseini, A. Al Khaled, A survey on the Imperialist competitive algorithm metaheuristic: Implementation in engineering domain and directions for future research, Appl. Soft Comput. 24 (2014) 1078–1094 [CrossRef] [Google Scholar]
 M.R. Ghaderi, M. Aghakhani, A. Eslampanah, K. Ghaderi, The application of imperialist competitive algorithm for optimization of deposition rate in submerged arc welding process using TiO_{2} nanoparticle, J. Mech. Sci. Technol. (2015) 357–364 [CrossRef] [Google Scholar]
 K. Deb, A. Pratap, S. Agarwal, T. Meyarivan, A fast and elitist multiobjective genetic algorithm: NSGAII, IEEE Trans. Evol. Comput. 6 (2) (2002) 182–197 [Google Scholar]
 B. Najlawi, M. Nejlaoui, Z. Affi, L. Romdhane, Multiobjective design optimization of the NBTTL Mechanism, in: International conference on acoustics and vibration, Tunisia, 2016 [Google Scholar]
 M. Vallerio, D. Telen, L. Cabianca, F. Manenti, J.V. Impe, F. Logist, Robust multiobjective dynamic optimization of chemical processes using the Sigma Point method, Chem. Eng. Sci. 140 (2016) 201–216 [CrossRef] [Google Scholar]
 M. Martinez, S. Kotov, D. De Vleeschouwer, D. Pas, H. Pälike, Testing the impact of stratigraphic uncertainty on spectral analyses of sedimentary series, Clim. Past 12 (2016) 1765–1783 [CrossRef] [Google Scholar]
 G.C. Marano, S. Sgobba, R. Greco, M. Mezzina, Robust optimum design of tuned mass dampers devices in random vibrations mitigation, J. Sound Vib. 313 (2008) 472–492 [CrossRef] [Google Scholar]
 P. Nimmegeers, D. Telen, F. Logist, J. Van Impe, Dynamic optimization of biological networks under parametric uncertainty, BMC Syst. Biol. 10 (2016) 86 [CrossRef] [Google Scholar]
Cite this article as: N. Mohamed, N. Bilel, Z. Affi, L. Romdhane, Multiobjective robust design optimization of the mechanism in a sewing machine, Mechanics & Industry 18, 606 (2017)
All Tables
All Figures
Fig. 1 The NBTTL mechanism. 

In the text 
Fig. 2 Position of a virtual loop of a fourbar linkage. 

In the text 
Fig. 3 The MOICA flowchart. 

In the text 
Fig. 4 Moving colonies toward their relevant imperialist. 

In the text 
Fig. 5 Imperialistic competition. 

In the text 
Fig. 6 The deterministic optimal solutions. 

In the text 
Fig. 7 Variation of the Angular Velocity Ratio for the three cases presented in Table 3. 

In the text 
Fig. 8 Variation of the needle jerk for the three cases presented in Table 3. 

In the text 
Fig. 9 The flowchart of MCS. 

In the text 
Fig. 10 Performance vulnerability of the solution SD1. 

In the text 
Fig. 11 Performance vulnerability of the solution SD2. 

In the text 
Fig. 12 Performance vulnerability of the Juki machine. 

In the text 
Fig. 13 The MOICA–MCS algorithm. 

In the text 
Fig. 14 The robust optimal solutions. 

In the text 
Fig. 15 Some configurations of robust solutions. 

In the text 
Fig. 16 Performance vulnerability of the solutions SD1 and SR1. 

In the text 
Fig. 17 Performance vulnerability of the solutions SD2 and SR2. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.