Issue 
Int. J. Simul. Multidisci. Des. Optim.
Volume 11, 2020



Article Number  1  
Number of page(s)  10  
DOI  https://doi.org/10.1051/smdo/2019020  
Published online  15 January 2020 
Research Article
Posture prediction and optimization for a manual assembly operation involving lifting of weights
^{1}
Department of Production Engineering, VSSUT, Burla 768018, Odisha, India
^{2}
Department of Mechanical Engineering, CET, Bhubaneswar 751003, Odisha, India
^{*} email: biswa.pp@gmail.com
Received:
26
January
2019
Accepted:
2
December
2019
The present work combines ergonomics with the posture prediction in the assembly process to avoid musculoskeletal issues of human operator. For improved productivity the operator should be in a better work environment and in sound health. The purpose of this paper is to provide a different perspective to avoid ergonomic risk factors in manual assembly. Here, a human is modeled as 20DOF as modeled in robotic analysis and simulated in a virtual environment. In the present study, two objective cost functions i.e. joint discomfort function and energy expenditure function have been employed for evaluating the optimized posture. For posture prediction, a combined multiobjective optimization (MOO) method is used and the objective cost functions are minimized i.e. less joint discomfort and less energy in MOO method required to do the manual assembly operation and consequently, the results are compared and finally the movements are tested using REBA technique.
Key words: Ergonomics / posture prediction / joint discomfort function / energy function / MOO / REBA
© B. Rout et al., published by EDP Sciences, 2020
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Precise simulation of assembly posture is necessary of the posture prediction for ergonomic assessment. However, manual manipulation of human models is required by the ergonomics simulation tools, which results in errors and less efficiency. The accuracy of the posture prediction method in nonautomated and semi automated workplace should satisfy the ergonomic condition. Most importantly, this method helps in providing a better understanding in predicting realistic postures for varied assembly conditions for the ergonomic analysis of the human performance.
This work proposes development of a generalized posture prediction method in assembly job simulation for the manual assembly tasks. At first, a human model is proposed, where the human body represents a system of kinematic chain with a number of revolute joints connecting the series of links. The model is same as in robotics; the open chain mechanism is used for analysis and application which are connected by joints serially to the links [1]. For finding a suitable posture prediction procedure, a human's anthropometric characteristic, for example, masses of body segments, body dimensions, are taken into consideration. We have also taken some assumptions to simulate the human body: they are flexion and extensions of the neck; rotational movement of trunk is not considered [2]. The discomfort and energy expenditure in joints are to be less in order to avoid musculoskeletal problems in human operators. In present work we have considered the human for posture prediction in different path in lifting weight in a packaging industry.
For the model to achieve the optimized posture prediction the use of multi objective optimization (MOO) technique used to find the effectiveness of the algorithm using multiple objectives in assembly task is described. In this work we have taken two objective functions that are energy expenditure function and joint discomfort function. The former deals with human postures to optimize effort of the job and the later deals with the survival capacity to predict fatigue of the body.
The joint limitations in human body parts, the work space limitations and reference point for the human are the constraints considered for optimization problem. The basic constrain taken for optimization of posture prediction problems is a hand that will be in contact with the targeted object. Also, we have considered different paths for lifting the weight.
Weight of the object, working distance, joint limitations, posture analysis and clearance analysis visibility, assembly force and clearance, are the common ergonomic requirements of a manual industry. In vision analysis parts, or components or operations must be able to see when assembled. In posture analysis weight of the object, working distance and assembly force are calculated. If weight limits are exceeded, a weight reduction of the part must be made or an adequate lifting device be developed.
Various researchers have employed methods which include experimental, statistical, and analytical approach in posture prediction of human with various degrees of freedom [3–13] for modeling. Inverse dynamics [3] and image processing tools [4] are also used by various researchers to predict posture and optimization using GA [9–11], ANN [15], MOO [16–23] to attain the posture prediction of human.
Rapid Entire Body Assessment (REBA) is ergonomic tool for posture analysis. REBA is developed using postural analysis and musculoskeletal risk factor in variety task with a quantitative indicator, calculating the average of each involved muscle activities in joints. In this work the optimized posture is validated using REBA [31].
From the literature reviews it is observed that, yet almost no study concerning 20DOF has been reported in posture prediction of human that would ensure scope for researchers. Novelty aspect, the present study focuses on a 20DOF human model for posture prediction. The objective cost functions taken are joint discomfort function and energy expenditure rate. MOO is used as optimization method to predict posture in assembly operation where different weights are to be lifted.
2 Human modeling
Human body consists of foot, leg, spine, torso, head, hand and wrist joined by various joints which have relative motion when humans do some work similar to the links of the robotic manipulator system. Depending upon DOF of joints and weight of the individual body parts each can be modeled as the link and joint of robotic arm [14].
In this work we have taken as 20 degrees of freedom human model for posture prediction as in our previous work. In the musculoskeletal system of human, different joints are considered as revolute joints joining the links or bones which are considered as rigid bodies. Different joints considered are ankle joint, knee joint, hip joint, trunk joint, shoulder joint, elbow joint and wrist joints. Out of which shoulder joint is modeled as 3 DOF, wrist is modeled as 2 DOF and other joints are 1 DOF each. In Figure 1 human model and the most neutral position is shown. In this work the movement of trunk and spine is not considered.
In this case the lengths of various joints for two subjects are: L_{1} is taken 460 mm and 400 mm, L_{2} is 400 mm and 380 mm, L_{3} is 140 mm and 120 mm, L_{4} is 680 mm and 580 mm, L_{5} is 360 mm and 320 mm, L_{6} is 190 mm and 150 mm, L_{7} is 340 mm and 300 mm, L_{8} is 310 mm and 260 mm, L_{9} is 190 mm each, L_{10} is 150 mm and 140 mm, respectively.
DH method is a proven method for mathematical modeling in the field of robotics. This mathematical modeling (DH method) is used for transformation of adjacent joint position with respect to the link. In this method various factors used are two link factors i.e. θ_{i} (link twist angle) and d_{i} (link distance) and two joint factors a_{i} (joint distance) and α_{i} (joint angle) to explain the position and orientation. In human biomechanics this DH method is used as in robotics assuming the body parts as link to be connected by joints. The standard DH algorithm is available in literature [24].
Neutral position is the most comfortable position of human and each joint is having a limiting value i.e. lower limit and upper limit. DH parameters and joint limitations are for the human model as shown in Table 1.
With the help of above parameters describe in Table 1 the position of different joints and relative angle between each body parts under consideration can be calculated which gives the posture prediction keeping the foot at fixed point and tip of the hand as the target point. The position of joint are depends upon the objective functions considered. The model is developed keeping the foot position unchanged, and changing the other joints to lift the job from its required position. The change in joint parameter is such that discomfort at the joint is less and with less energy expenditure.
Fig. 1 Showing the human model with neutral position. 
DH parameter and limiting position for 20 DOF model.
3 Posture prediction
In this work a 20 DOF model is used to predict the posture of human operator for a lifting of weight in a manual assembly work. For this we have taken two objective cost functions, they are joint discomfort function and energy expenditure function.
3.1 Joint discomfort function
Discomfort function generally measures the difference of CG to the desired position with respect to the neutral position. In Table 1 each joint has a lower limitation and upper limitation. The discomfort value is small when it is at neutral position, if the joints are moving from the neutral position to the limiting value, the discomfort level increases. The basic aim is to reduce this discomfort level.
It is found that the joint discomfort function is minimum in its neutral position where no movement is there and is higher when going towards the upper limit of the joint and lower limit of the joint. QU and QL are especially designed provisions for the upper limit and the lower limit of each joint. The neutral position of the digital human is given in Figure 1. The values at neutral positions and at the upper/lower limits of each joint are specified in Table 1. QU and QL are defined mathematically and discomfort function are derived [25,26].(1) (2) (3) (4)
3.2 Energy function
The total rate of energy expenditure (E), which is expressed in Watts per kilogram, is defined as follows:(5)
where E_{W} = mechanical power of muscle, E_{B} = basal metabolic rate (BMR).
Normally, in the joint space only the force is considered for energy expenditure in human model. Following Kim et al. [28], the definition of mechanical power is given as dot product of joint torque (τ_{i}) and joint velocity (q_{i}). The total mechanical power E_{M} is as follows:(6)
In Table 2, joint torque limits are given.
The BMR i.e. basal metabolic rate is the metabolic rate of a human being at rest [27]. It means the minimum quantity of energy necessary for functioning of human, without doing any kind of external work. Hase et al. [29] has presented the following BMR model:(7)where BW is the body weight (kilogram).
For optimization of posture prediction, the method is used, which determines different position of joints with respect to reference point such that the objective cost function is minimum. This can be expressed mathematically as(8)where k, Total number of objective functions, m, Total number of inequality constraints, e, Total number of equality constrains, q ϵ E ^{DOF} = vector of design variables.
The optimal posture can be obtained by solving the optimization condition. For this case the problem can be defined as:(9)
In this work, cost functions refer to joint discomfort function and energy function. qi ^{U} and qi ^{L} define the upper and lower limits for qi joint in that order. Limits are given for not to assume an unrealistic posture [24].
Multi Objective Optimization (MOO) method is a multi criteria decisionmaking process. MOO is used where more than one objective function are there. It is a welldeveloped method for posture prediction. Here we have taken two objective functions: joint discomfort function and energy expenditure function. For MOO method for posture prediction is redefined as [30]:(10)where, f _{1} N is the joint discomfort function and f _{2} N is the energy function. ε is the positive infinitesimal value which is equal to zero. The P (hip) can be found out, joint movement and anthropometric data. Weight factors w _{1} and w _{2} are given to discomfort function and energy expenditure function, respectively.
For optimization, MATLAB ® optimization toolbox is used. The function “fmincon” is used to find the minimum discomfort function and energy expenditure function. It finds the constrained minimum of a scalar function of different variables from initial point to target point. It uses SQP algorithm.
Joint torque limits.
4 Procedure
Posture prediction of manual assembly is highly required in human models. The predicted posture required for a particular task has to satisfy a large number of constraints due to the task condition. In manual assembly visual demand has a great impact. Visual demand gives clear picture of location of the task to secure the assembly efficiency and the quality. Here in this work same object is to be lifted by two subjects having weight 75 kg and 58 kg, respectively, taken for assembly operations. Here we have taken a load of 3 kg to be lifted by two subjects to do a particular task.
The starting point and target points are same for both the subjects as the height and length of the body parts are different for both the subjects. The objective cost functions described in previous section gives rise to different postures which is under consideration in this work.
Figure 2 gives the midsagittal plane along with angulations. As there is no distance between hand and object, P (hand) is considered to be the object point. The distance between P (eye) and object is known as visual demand (d). When the point P (hip) is found out the virtual demand is obtained by inclination q_{4}. In this case P (hip) controls q_{1}, q_{2}, q_{3} and q_{4}.
Fig. 2 Model with angulations and visual demand. 
5 Results and discussion
To test and analyze the procedures described in previous section the different path followed by the operator considered in Figure 3 shows initial position as point 1 and the final point as point 7. The work piece is to be lifted from point 1 and be kept to point 7. For this purpose, we have taken three paths. The first path is specified by points 1234567, which is the longest path, where 1234 is the vertical path and points 4567 represents horizontal path. The second path is specified by 1121314157, which is the shortest path. Path three is the in between path of both path1 and path2. Path three is specified by 18910117.
At first the posture prediction is used with optimized discomfort function in all three paths. Discomfort function in all the three paths is shown in Figure 4 for both the subjects. In this figure the discomfort value obtained from the objective function are depicted. In Figure 4a, it is described for path1 for both the subjects. Similarly, in Figure 4b and c it is described for path2 and path3, respectively. From the results which are presented in Figure 4a–c for subject1 discomfort is higher in path1 and path2 and it is opposite in path3 i.e. discomfort is higher in subject2. It is established that the discomfort is higher in initial conditions and it is found to be gradually decreased.
Then the posture prediction is used with optimized energy expenditure function in all three paths. In Figure 5 energy functions for three paths are shown for both the subjects. In this figure the energy expenditure value obtained from the objective function are depicted. In Figure 5a it is described for path1 for both the subjects. Similarly, in Figure 5b and c it is described for path2 and path3, respectively. From the results in Figure 5a–c, it is evident that higher level of energy expenditure is in subject1 as compared with subject2. As far as energy expenditure is concerned it is dependent on hip angle. If hip angle i.e. Ө_{4} increases, energy expenditure decreases and if, Ө_{4} decreases energy expenditure decreases, so energy expenditure is directly proportional to hip angle.
In above figures discomfort function and energy expenditure for both subjects in different paths are represented. For subject1 discomfort functions, energy function and the multi objective is shown for all the three paths in Figure 6.
Figure 6a–c gives the discomfort function, energy expenditure function and multi objective function for subject1 for three different paths i.e. path1, path2 and path3. Figure 7a–c gives the discomfort function, energy expenditure function and multi objective function for subject2 for three different paths i.e. path1, path2 and path3. From above figures normalized graphs of objective cost functions, it is found that when discomfort function is considered the curve lies in the top of each figure, while for energy expenditure function, the curve lies on the bottom of graph. When multi objective is considered the curve lies in between discomfort and energy function. So, in optimization of discomfort, we neglected the total energy required which implies the energy consumption is high. Similarly, in energy expenditure optimization curve, energy expenditure is optimized but not the discomfort, and discomfort function is high. In MOO, both the objective functions are considered and the curve lies in between. So, MOO should be used for prediction of posture. And also, it is marked that when discomfort function increases energy expenditure decreases and if discomfort decreases energy expenditure increases. From Figure 8a–c it is evident that for this particular task subject2 is better as compared to subject1 that implies a short heighted person is better for this work. These postures are verified by REBA technique. Pareto front is given for both the subjects in Figures 9 and 10.
REBA procedure involves dividing body segments to coding with respect to different planes. Scoring system for muscle activities is determined for different postures. Coupling activity of load carrying capacity is taken into consideration. REBA scoring system is given by Hignett and McAtamney (2000) [31].
In Table 3, REBA score is given for three different paths taken by the operator for the same work. Initial point score is same for all the paths and the final point score is also same for all three paths. Scores of intermediate points are different for all. From the above table the REBA scores are in between medium risk for all path points.
Fig. 3 Paths of operation. 
Fig. 4 Discomfort function in 3 paths. 
Fig. 5 Energy expenditure in 3 paths. 
Fig. 6 Multi objective for subject1 in path1 (a), path2 (b), path3 (c). 
Fig. 7 Multi objective for subject2 in path1 (a), path2 (b), path3 (c). 
Fig. 8 Multi objective for path1 (a), path2 (b) and path3 (c). 
Fig. 9 Pareto front for Subject1. 
Fig. 10 Pareto front for Subject2. 
REBA score in different path. 1 = Negligible risk; 2–3 = Low risk, change may be needed; 4–6 = Medium risk, further investigation, change soon; 7–10 = High risk, investigate and implement change; 11+ = Very high risk, change immediately.
6 Conclusion
In this work a model with 20 DOF of human model is considered for lifting the load in a workspace within reachable limit of the standing operator. The model is developed keeping the ankle joint as the reference point and assumed to be fixed, and changing the other joints to lift the job from its required position to desired position. The work uses to find out the posture of the test examples with these objective measures separately and also in combination as MOO case. The two subjects (of different weight and height) are considered as test examples. Three procedures are used to find out the optimized posture in both the test cases. It is found that the discomfort function is only the measure of the joints of the forearm part as it is found that the hip position almost remains of the same height. In minimum energy expenditure rate hip is at a lower position and using MOO the hip remains in between the both previous parts. The change in joint parameter is such that discomfort at the joint is less and with less energy expenditure which also depicts that the movement is under medium risk in the REBA scale.
Based on the extensive experimental work the following outcomes can be drawn. A short heighted person is better for this operation. The multi objective is a better procedure for posture prediction. When discomfort function increases energy expenditure decreases and if discomfort decreases energy expenditure increases. This paper can be helpful in managing the operators in manual assembly operations in various workplaces. In the extended work a greater number of cost functions can be taken for analysis in different work places.
References
 B. Rout, P.P. Tripathy, R.R. Dash, D. Dhupal, Optimization of posture prediction using MOO in brick stacking operation, in Computational Intelligence in Data Mining (Springer Science and Business Media LLC, Cham, 2020) [Google Scholar]
 B. Rout et al., Optimal posture prediction in brick stacking operation for reducing ergonomic risk factor, Int. J. Recent Technol. Eng. 8, 1835–1841 (2019) [Google Scholar]
 V. Lippi, Prediction in the context of a humaninspired posture control model, Robot. Auton. Syst. 107, 107–63 (2018) [CrossRef] [Google Scholar]
 S.D. Farahani et al., Human arm posture prediction in response to isometric endpoint forces, J. Biomech. 48, 48–4178 (2015) [Google Scholar]
 J.H. Kim et al., Lifting posture analysis in material handling using virtual human, in Proceedings of International Mechanical Engineering Congress and Exposition IMECE2005, November 5–11, 2005, Orlando, Florida USA [Google Scholar]
 B. Rout et al., Optimization of posture analysis in manual assembly, Int. J. Mech. Prod. Eng. Res. Dev. 8, 751–764 (2018) [Google Scholar]
 A.M. Karim, Y. Jingzhou, M. Timothy, B. Steven, M. Anith, Z. Xianlian, P. Amos, A. Jasbir. Towards a new generation of virtual humans, Int. J. Human Factors Model. Simul. 1, 1–2 (2006) [Google Scholar]
 J. Yang, T. Sinokrot, M.K. Abdel, A general analytic approach for santosa upper extremity workspace, Comput. Ind. Eng. 54, 54–242 (2008) [Google Scholar]
 B. Rout et al., Effective work procedure design using discomfort and effort factor in brick stacking operation – a case study, IOP Conf. Ser. Mater. Sci. Eng. 310 , 012020 (2018) [CrossRef] [Google Scholar]
 A.K. John, K. Krishnakumar, Performing multiobjective optimization on perforated plate matrix heat exchanger surfaces using genetic algorithm, Int. J. Simul. Multisci. Des. Optim. 8 , A3 (2017) [CrossRef] [Google Scholar]
 M. Chandramouli, G.R. Bertoline, A genetic algorithm based optimization framework to visualize, evaluate, and modify 3D space configurations in Desktop VR, Int. J. Simul. Multisci. Des. Optim. 5 , A01 (2014) [CrossRef] [Google Scholar]
 S. Pheasant, C.M. Haslegrave, Bodyspace: Anthropometry, Ergonomics, and the Design of Work (Taylor and Francis, Boca Raton, FL, 2006) [Google Scholar]
 K. Hase, N. Yamazaki, Development of threedimensional wholebody musculoskeletal model for various motion analyses, JSME Int. J. Ser. C 40, 40–25 (1997) [Google Scholar]
 D.B. Chaffin, Improving digital human modelling for proactive ergonomics in design, Ergonomics 48, 48–478 (2005) [CrossRef] [Google Scholar]
 A. Gholipour, N. Arjmand, Artificial neural networks to predict 3D spinal posture in reaching and lifting activities: applications in biomechanical models, J. Biomech. 49, 49–2946 (2016) [CrossRef] [Google Scholar]
 R.T. Marler, J.S. Arora, J. Yang, H.J. Kim, K. AbdelMalek, Use of multiobjective optimization for digital human posture prediction, Eng. Optim. 41, 41–925 (2009) [CrossRef] [Google Scholar]
 L. Ma, W. Zhang, D. Chablat, F. Bennis, F. Guillaume, Multiobjective optimisation method for posture prediction and analysis with consideration of fatigue effect and its application case, Comput. Ind. Eng. 57 , 1235–1246 (2009) [CrossRef] [Google Scholar]
 R.T. Marler, J.S. Arora, Survey of multiobjective optimization methods for engineering, Struct. Multidiscip. Optim. 26, 26–369 (2004) [CrossRef] [MathSciNet] [Google Scholar]
 Y. Xiang et al., Human lifting simulation using a multiobjective optimization approach, Multibody Syst. Dyn. 23, 23–431 (2010) [CrossRef] [Google Scholar]
 Z. Mi, J. (James) Yang, K. AbdelMalek, Optimizationbased posture prediction for human upper body, Robotica 27, 27–607 (2009) [Google Scholar]
 E.N. Horn, Optimizationbased dynamic human motion prediction, Master Thesis, University of Iowa, 2005 [Google Scholar]
 F.C. Anderson, M.G. Pandy, Static and dynamic optimization solutions for gait are practically equivalent, J. Biomech. 34, 34–153 (2001) [CrossRef] [PubMed] [Google Scholar]
 L. Ma, W. Zhang, D. Chablat, F. Bennis, F. Guillaume, Multi objective optimisation method for posture prediction and analysis with consideration of fatigue effect and its application case, Comput. Ind. Eng. 57, 57–1235(2009) [Google Scholar]
 J. Denavit, R.S. Hartenberg, A kinematic notation for lowerpair mechanisms based on matrices, ASME J. Appl. Mech. 23, 23–215 (1955) [Google Scholar]
 R.T. Marler, S. Rahmatalla, M. Shanahan, K. AbdelMalek, A new discomfort function for optimizationbased posture prediction, in Proceedings of the SAE Human Modeling for Design and Engineering Conference, 2005, p. 2680 [Google Scholar]
 Q.L. Zou, Q.H. Zhang, J.Z. Yang, J. Gragg, An inverse optimization approach for determining weights of joint displacement objective function for upper body kinematic posture prediction, Robotica 30, 30–389 (2012) [Google Scholar]
 B.R. Umberger, K.G. Gerritsen, P.E. Martin, A model of human muscle energy expenditure, Comput. Methods Biomech. Biomed. Eng. 6, 6–99 (2003) [CrossRef] [Google Scholar]
 J.H. Kim, K. AbdelMalek, J. Yang, T. Marler, K. Nebel, Lifting posture analysis in material handling using virtual humans, Proc. ASME Int. Conf. Manuf. Sci. Eng. 16, 16–1445 (2005) [Google Scholar]
 K. Hase, N. Yamazaki, Development of threedimensional wholebody musculoskeletal model for various motion analyses, JSME Int. J. Ser. C Dyn. Control Robot. Des. Manuf. 40, 40–25 (1997) [Google Scholar]
 L. Ma, W. Zhang, D. Chablat, F. Bennis, F. Guillaume, Multiobjective optimisation method for posture prediction and analysis with consideration of fatigue effect and its application case, Comput. Ind. Eng. 57, 57–1235 (2009) [Google Scholar]
 S. Hignett, L. McAtamney, Rapid Entire Body Assessment (REBA), Appl. Ergon. 31, 31–201 (2000) [CrossRef] [Google Scholar]
Cite this article as: Biswaranjan Rout, Rati Ranjan Dash, Debabrata Dhupal, Posture prediction and optimization for a manual assembly operation involving lifting of weights, Int. J. Simul. Multidisci. Des. Optim. 11, 1 (2020)
All Tables
REBA score in different path. 1 = Negligible risk; 2–3 = Low risk, change may be needed; 4–6 = Medium risk, further investigation, change soon; 7–10 = High risk, investigate and implement change; 11+ = Very high risk, change immediately.
All Figures
Fig. 1 Showing the human model with neutral position. 

In the text 
Fig. 2 Model with angulations and visual demand. 

In the text 
Fig. 3 Paths of operation. 

In the text 
Fig. 4 Discomfort function in 3 paths. 

In the text 
Fig. 5 Energy expenditure in 3 paths. 

In the text 
Fig. 6 Multi objective for subject1 in path1 (a), path2 (b), path3 (c). 

In the text 
Fig. 7 Multi objective for subject2 in path1 (a), path2 (b), path3 (c). 

In the text 
Fig. 8 Multi objective for path1 (a), path2 (b) and path3 (c). 

In the text 
Fig. 9 Pareto front for Subject1. 

In the text 
Fig. 10 Pareto front for Subject2. 

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.