Research Article  Open Access
Process Design of a Ball Joint, Considering Caulking and PullOut Strength
Abstract
A ball joint for an automobile steering system is a pivot component which is connected to knuckle and lower control arm. The manufacturing process for its caulking comprises spinning and deforming. In this study, the process was simulated by flexible multibody dynamics. The caulking was evaluated qualitatively through numerical analysis and inspecting a plastically deformed shape. The structural responses of a ball joint, namely, pullout strength and stiffness, are commonly investigated in the development process. Thus, following the caulking analysis, the structural responses were considered. In addition, three design variables related to the manufacturing process were defined, and the effects of design variables with respect to pullout strength, caulking depth, and maximum stress were obtained by introducing the DOE using an L_{9} orthogonal array. Finally, the optimum design maximizing the pullout strength was suggested. For the final design, the caulking quality and the pullout strength were investigated by making six samples and their tests.
1. Introduction
An automobile ball joint is a pivoting element used to allow rotational motion between the parts of the steering and suspension systems. Depending on the car development, in general, most automobile ball joints are made of a socket, bearing, plug, and ball stud. The manufacturing process assembling these parts is called the caulking process.
Existing studies [1–4] have tried to simulate the caulking process, because it can affect the performance of the ball joint. Based on the simulation results, it is possible to evaluate the caulking process, qualitatively, by inspecting the plastically deformed shape. The numerical analysis of a ball joint is a set of coupled problems of rigid body and flexible body. In this research, 3D flexible multibody dynamic analysis for a caulking process using cold forging was simulated by DAFUL [5], which adopts an implicit integration method. In general, it is known that results calculated from an implicit solver have less noise and are more stable [5–7].
The structural responses of stiffness and pullout strength were then calculated, to check if the design satisfies the related requirements. The analysis was sequentially performed, following the caulking process. In this process, the deformation and stress results obtained from the analysis were saved. Sequential analysis has a strong advantage, in that it can be analyzed by considering the deformed shape and residual stress. The axial and lateral forces were applied to obtain its stiffness. This axial stiffness and this lateral stiffness can be measured by the axial displacement and lateral displacement, respectively, which are generated in the ball stud when the specified loads are applied. The pullout strength means the required force to pull the ball stud out from the ball joint assembly. A low pullout strength can deteriorate the structural stability and safety performances [2, 3].
Three process variables in caulking were set up as design variables. These were the revolution speed, downward distance, and downward time of the roller. The process variables in the existing products were determined through empirical case study, not by the systematic design methodology proposed in the research. That is, the process variables in the existing products were set up by changing them 5 to 10 times, performing dynamic analyses, and selecting one adequate case. The effects of design variables with respect to pullout strength, deformed depth of socket, and maximum stress generated in the socket were obtained by introducing the DOE. The DOE is a statistical technique used to study the effects of multiple variables simultaneously [8–10]. In this research, the extent of the DOE is confined to conducting orthogonal array experiments and the ANOVA and selecting optimum levels. Furthermore, conducting experiment in this study means performing numerical simulation. In this research, orthogonal array was utilized, in which it is assumed that there is no interaction between design variables.
Then, nine flexible multibody dynamic analyses using DAFUL were performed. Based on the simulation results, the design sensitivity over the design range was determined. In addition, the optimum process variables were suggested. Because of such design characteristics and the calculation time of flexible multibody dynamics, no numerical optimization algorithm could be applied. The process variables in this study were suggested by applying DOE scheme, which could replace 27 dynamic analyses with 9 analyses thus reducing the number of analyses by onethird while achieving the similar result. This research focuses on the improvement of structural responses, considering the caulking process and guaranteeing the caulking quality.
2. Initial Design and Analysis of the Ball Joint
2.1. Initial Design of the Ball Joint
The ball joint of the present study is the part mounted on a pickup truck being produced at company. The ball joint investigated in this research is a part that is connected to a knuckle and a lower control arm. The ball joint serves as a flexible pivot element for the steering system. The initial design of the ball joint is shown in Figure 1. The ball joint is made of a socket, bearing, plug, and ball stud as shown in Figure 1.
A plug prevents the components from being separated during manufacture and operation of the ball joint. The bearing, which has relatively much less stiffness among the components of the ball joint, acts as lubrication and buffering. The socket serves as the body of the temporarily assembled ball joint and plays a role in covering the interior parts with the plug through plastic deformation. The ball stud, which induces rotation in all directions, is made by assembling the upper ball and bearing. Each part and the assembly shape of the ball joint are represented in Figure 1 [1–3].
2.2. Equations for Dynamic Analysis [2, 6, 7]
The virtual work, , done by a generalized force is represented as where is the virtual displacement and is the generalized force. The Lagrange multiplier method yields the governing equation for dynamic analysis, called the equations of motion, as follows: where is the mass matrix, is the stiffness matrix, is the constraint equation, and is a Lagrange multiplier. The positive level constraint equations are represented as Then, the following equations for each step are defined by applying the tangent space method to (2) and (3): To solve (4), NewtonRaphson method is applied. Thus, The solution is updated as
2.3. Flexible Multibody Modeling of the Ball Joint
The parts of the ball joint were modeled with hexahedral elements as shown in Figure 2, while two rollers were modeled as rigid bodies. The number of elements of the socket, bearing, plug, and ball stud is 9,790, 3,598, 1,620, and 17,669, respectively. The materials of the socket, bearing, plug, and ball stud are SM45C, nylon, SPC1, and SCM435 [1–3], respectively. The stressstrain curve of each material is represented in Figure 3.
For the contact of the ball joint, there are the contact surfaces between 4 parts, composed of the bearingball stud, socket bearing, socket plug, and bearing plug. Each contact surface is defined as a 3dimensional side. These contacts are defined as “Flex to Flex” to define the contact condition between flexible bodies in DAFUL. In contrast, the contact between roller and socket is defined as “Flex to Rigid” to consider the contact surface between the roller, modeled as a rigid body, and the socket modeled as a flexible body.
2.4. Simulation for the Caulking Process and Prediction of the PullOut Strength and Stiffness
The manufacturing process of a ball joint is called the caulking process. A caulking machine, as shown in Figure 4, is used to assemble the parts of a ball joint. The parts of the socket, bearing, plug, and ball stud are sequentially positioned in the caulking machine. The lower part of the socket is fixed to a jig. Then, two rollers push down in a vertical direction, to compress the assembled ball joint and plastically deform the upper part of the socket. The roller is set to 300 rpm. The caulking process was simulated by using DAFUL. The movement of each roller is constrained to a CJ (cylindrical joint). The boundary conditions for caulking analysis are represented in Figure 5(a).
(a) Caulking machine
(b) Assembled parts
(a) Boundary conditions
(b) Stress contour
(c) Time versus stress curve at the socket
The caulking analysis of the initial design requires 38 hours of runtime on a 3 GHz PC. The maximum stress in the ball joint was 871 MPa, which was generated at the surface of the socket. The value was determined from the plastic deformation region. In contrast, the maximum stress at the plug was 130 MPa. Both stresses are shown in Figure 5(b). The stress variation due to time at the socket is represented in Figure 5(c). The ball stud should be securely and stably held in the socket after the caulking process. The performance can be evaluated qualitatively through inspection. Based on the deformation results, it was seen that the caulking of the initial design was successful.
The force to pull out the ball stud of the assembly must be greater than a specified value in order to maintain car performance. The specified value is determined according to the design requirement set by the manufacturer. If the force required to pull out the ball stud is below the specified value, it is considered as a risk that the ball joint could fail during driving. That would make the car lose suspension performance of a lower control arm. The boundary condition for the pullout strength analysis, initial shape, and pulled out shape is shown as in Figures 6(a), 6(b), and 6(c), respectively.
(a) Boundary condition
(b) Initial shape
(c) Pulled out shape
(d) Time versus contact force curve
The specified displacement was given to the center of the ball joint as the loading condition. In this analysis, the contact force between the ball stud and plug was defined as the pullout force when the ball stud started to be pulled out. The maximum contact force was calculated as 39 kN. The variation of contact force according to time is represented in Figure 6(d), which satisfies the requirement.
Axial and lateral forces are applied to the ball stud, to determine the stiffness of the ball joint. The two force conditions are represented in Figures 7(a) and 7(b). The displacement variations are shown in Figures 7(c) and 7(d). The maximum displacement of each case is much less than its criterion, which is very marginal with respect to the allowable value, therefore satisfying the design requirements for both types of stiffness. Thus, the stiffness performance was excluded from the design constraint, when the DOE was applied.
(a) Axial stiffness condition
(b) Lateral stiffness condition
(c) Time versus axial displacement curve
(d) Time versus lateral displacement curve
3. Application of Design of Experiments to CAE Based Design
3.1. Definition of Process Variables and Determination of Responses
The roller speed, , the downward distance of the roller, , and the downward time, , were set up as the process variables to find an optimal condition. The downward distance of the roller means the length between the lowest point of the roller and the highest point of the socket. Thus, the process variables, and , determine the descending velocity of the roller.
The most interesting response was the pullout strength, which was measured by force unit. In addition, the caulking depth and the maximum stress were considered as characteristics. The caulking depth is the bent length of the socket. The longer it is, the better the quality of the caulking becomes. In contrast, the maximum stress is the stress generated in the plastically deformed socket.
3.2. DOE Using Orthogonal Array
The overall design process to find an optimum condition is as follows. First, the smallest orthogonal array was selected as the orthogonal array with minimum experiments that can assign all the design variables to their columns. For each process variable, the number of levels was set to three. The second level was set up as the initial variable. The first and third levels were fixed by the lower and upper ones, respectively, around the initial value [8].
The levels of design variables for an orthogonal array were determined as shown in Table 1. Then, an appropriate orthogonal array was selected. For a problem with three variables and three levels, the orthogonal array is recommended, in which two variables are assigned for the first two columns and the remaining variable for the fourth column. The orthogonal array can replace 3^{3} fullcombinational experiments.

Second, 9 dynamic analyses were performed as indicated in Table 2. An experiment in this study means one dynamic analysis using DAFUL [6]. It is assumed that the interaction caused by variables can be ignored. That is because just one caulking analysis requires computation time of between 35 and 45 hours, depending on the contact conditions. Almost 15 days were required to finish 9 analyses using DAFUL, including the caulking and pullout analyses. Each response was obtained as in Table 2.

Third, based on the analyses results, the relative importance of each variable to each response could be obtained through an ANOVA. The ANOVA tables for pullout strength, caulking depth, and maximum stress are shown as Tables 3, 4, 5, 6, 7, and 8, in which significance was established at the 95% confidence level. For the pullout strength, it is seen from Tables 3 and 4 that the most sensitive variable is , and the variable is insignificant. Therefore, for the caulking depth performance, the variable is considered as the major influence. Considering the maximum stress, only the variable was determined as the significant one. However, the variation of the stress values over the interested range is not large, in comparison with the other responses. Thus, the response of the maximum stress was not considered, when determining an optimum process. The average characteristics for each level called the factor effects are shown in Figure 8. For example, the average pullout force for the 1st level of design variable was calculated as (22.5 + 19.3 + 8.1)/3. It means that the design variable with large difference between levels has a large influence on the experiments.






(a) Effect for pullout force
(b) Effect for deformed depth
(c) Effect for maximum stress in socket
The fourth step is to determine the optimum level of each design variable. For the pullout force, the optimum level is determined as the level with maximum mean from Table 9. Considering only the pullout strength and the significant variables, an optimum setting is , since the level with the largest pullout force can be an optimum level. However, an optimum setting considering only the caulking depth is . The process variable has the same optimum level for the pullout strength and the caulking depth. Then, the process variable has an inconsistent optimum level. However, it is seen that the main effect of for the pullout force is much larger than that for the caulking depth. Thus, an optimum setting of significant variables was selected as . Also, the variable was selected as . Thus, the optimum levels were determined as . Then, the predicted response at the optimum setting was determined as [8–10] where is a predicted value of characteristic, is the overall mean, and , , and are mean of the main effects of , , and at , and th level, respectively. The predicted response determined from (7) and the true responses determined from DAFUL are summarized in Table 10. Table 10 shows that the predicted responses at the optimum calculated through DOE are comparable with the true values determined from DAFUL. It can be seen that the difference between the estimated value and the true value is relatively small. Both of them satisfy the requirement specified by the manufacturer.


3.3. Test Results for PullOut Strength
Based on the suggested optimum design, six samples were made to validate the simulation results of the caulking quality and the pullout strength. From the deformation shapes of the sockets, it was seen that the caulking processes in all the samples were successful. Then, the pullout strength test was carried out for six specimens. The loading and boundary conditions are the same as the ones of the simulation. The pullout forces for six specimens are summarized in Table 11. The error between the average value obtained from six tests and the simulation result is about 17%. The specimens after pullout strength test are shown as in Figure 9.

4. Conclusions
In this study, the caulking process of a ball joint mounted on a pickup truck was simulated by applying the flexible multibody dynamic software called DAFUL. Then, the structural responses of pullout strength and stiffness were calculated, following the caulking analysis. Sequential analysis has the strong advantage that it can be analyzed by considering the deformed shape and residual stress. Furthermore, three variables in the caulking process were defined as design variables to find their optimum settings. The pullout strength and the caulking depth were considered in this design stage. Because the flexible multibody dynamic analysis for the caulking process requires a much computation time, an orthogonal arraybased design was adopted to economize the computation time. By performing experimental design, the sensitivity of each process variable to each characteristic was obtained. Then, by the rule of the DOE, the optimum setting of three variables was determined as A_{2}B_{3}C_{1}. Based on the suggested optimum design, the caulking quality and the pullout strength were investigated by making six samples and their tests. The method used in this research can achieve the equivalent result using the existing method requiring twothird less number of analyses.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research was financially supported by the Ministry of Education Science and Technology (MEST) and the National Research Foundation of Korea (NRF) through the Human Resource Training Project for Regional Innovation (2012H1B8A2026078).
References
 K. H. Lee and S. C. Hwang, “Structural dynamic analysis of a ball joint,” AIP Conference Proceedings, vol. 1499, pp. 394–398, 2012. View at: Google Scholar
 B. H. Jang and K. H. Lee, “Analysis and design of a ball joint, considering manufacturing process,” Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, vol. 228, no. 1, pp. 146–151, 2014. View at: Google Scholar
 B. H. Jang, B. S. Sin, and K. H. Lee, “The process design of a ball joint using the DOE, considering caulking and pullout strength,” in Proceedings of the International Symposium on Green Manufacturing and Applications, Honolulu, Hawaii, USA, June 2013. View at: Google Scholar
 J. Park, S. Choi, K. Na, and Y. Kim, “Application of FE analysis for optimal design of caulking process,” Journal of Materials Processing Technology, vol. 198, no. 1–3, pp. 471–477, 2008. View at: Publisher Site  Google Scholar
 Virtual Motion, DAFUL 4.1 User's Manual, 2012.
 C. H. Lee and D. S. Bae, “A parametric generalized coordinate formulation for flexible multibody system dynamics,” in Proceedings of the ECCOMAS Thematic Conference on Multibody Dynamics, vol. 29 June2, Warsaw, Poland, June 2009. View at: Google Scholar
 J. K. Lee, C. H. Lee, and D. S. Bae, “A parametric generalized coordinate formulation for mechanical joint elements in dynamic systems,” Proceedings of the Institution of Mechanical Engineers C: Journal of Mechanical Engineering Science, vol. 228, no. 6, pp. 1063–1076, 2014. View at: Publisher Site  Google Scholar
 K. H. Lee, J. W. Lee, J. S. Park, and G. J. Park, “An optimization algorithm using orthogonal arrays in discrete design space for structures,” Finite Elements in Analysis and Design, vol. 40, no. 1, pp. 121–135, 2003. View at: Publisher Site  Google Scholar
 M. Teruo, Taguchi Methods Benefits, Impacts, Mathematics, Statistics, and Applications, ASME, New York, NY, USA, 2011.
 R. Mead, Statistical Principles for the Design of Experiments: Applications to Real Experiments, Cambridge University Press, New York, NY, USA, 2012.
Copyright
Copyright © 2014 BongSu Sin and KwonHee Lee. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.