/ / Article
Special Issue

## Complex System Modelling in Engineering Under Industry 4.0

View this Special Issue

Research Article | Open Access

Volume 2021 |Article ID 2071282 | https://doi.org/10.1155/2021/2071282

Yingda Li, "Simplified Beam Element Model of Badminton Batting Process Based on Motion Differential Equation", Complexity, vol. 2021, Article ID 2071282, 10 pages, 2021. https://doi.org/10.1155/2021/2071282

# Simplified Beam Element Model of Badminton Batting Process Based on Motion Differential Equation

Revised12 May 2021
Accepted17 Jun 2021
Published25 Jun 2021

#### 1. Introduction

Finite element analysis is a process of solving the essence of a simple problem after it is complicated. It can simulate the load conditions of a real physical system and so on. It is a practical application of mathematical methods. Applying it to the finite element analysis of badminton racket and badminton, it can accurately simulate the working condition, mechanics, displacement changes of the moment when badminton ball slaps badminton racket, especially for the analysis of the complex process of hitting badminton, which has the characteristics of high accuracy and wide applicability. As the process of badminton hitting mainly involves displacement and force, the mathematical expression of the relationship between force and displacement can be accurately carried out by using the equation of motion [3]. By combining the finite element analysis method with the differential equation of motion, the mechanical analysis of badminton batting process can clearly show the relationship between the load and displacement in the process of badminton batting [4].

The main contributions of this paper are as follows: firstly, the finite element model of string and beam elements of badminton racket is constructed, which can be used to analyze the micromechanics and displacement of badminton batting process. Secondly, based on the established finite element model, the static numerical calculation of badminton racket net surface is carried out, and the net tension of badminton beam element is simulated.

The innovation of this research mainly has three points: the first is to construct the dynamic differential equation of badminton and carry out the mechanical analysis of the badminton batting process using badminton racket; the second is to construct the badminton net model and verify the correctness and feasibility of the net model through the simulation calculation of beam element and rod element; the third is to analyze the different badminton characteristics and the influence of the tension on the return speed of badminton. This research can let people from a more professional point of view analyze the new badminton batting moment; also it provides a certain theoretical basis for the production of badminton racket. At the same time, it is feasible to apply the finite element analysis method to the process analysis of badminton racket hitting, which also provides a reference for the process analysis of similar sports. From the perspective of theoretical analysis, this paper provides ideas for the finite element analysis of sports equipment. It can be seen that the finite element analysis method can achieve good results in the design and optimization of sports equipment.

In modern industrial applications, the finite element analysis method is also widely used. Based on the finite element analysis software, Bhelsaikar et al. designed a new type of continuous manipulator. Through the finite element analysis software, it is helpful for the prototype manufacturing and testing of robot gripper [8]. Badia et al. proposed an enhanced finite element space technology based on cell focusing technology to solve the ill-posed problems in the traditional three kinds of nonadaptive finite element technology and obtained the optimal finite element solution by dividing different elements [9]. Hiptmair et al. proposed a finite element method that can be applied to incompressible fluid dynamics equations. The second-order semi-implicit time stepping method is used to establish the energy law. Through numerical verification, the correctness of theoretical prediction and the effectiveness of preprocessor of the method are also illustrated [10]. Wang et al. carried out finite element analysis on the manufacturing process of dielectric materials in electromagnetic engineering. Due to the comprehensive influence of processing technology, environmental temperature, and personnel operation, there are many spatial uncertainties in the manufacturing of dielectric materials in electromagnetic engineering. After introducing the interval field model, the uncertainty of this material can be described, and the upper and lower bounds of electromagnetic response can be calculated by using the interval finite element method [11]. Spherical material is difficult in engineering design. The design of circular material plate by finite element analysis is helpful to test and measure with finite element analysis software, so as to understand the stiffness of the designed material [12]. Karthikeyan and Chandru analyzed the wave propagation of annular elliptic thin plate with different voltage materials. By compiling the computer program of finite element method, the numerical calculation of different shell elements was carried out. The model is simple and effective and is suitable for the material design problems that need to be solved accurately [13]. Andena et al. constructed a three-dimensional finite element model of track and field runway and carried out quasistatic compression test on laboratory samples through the model, so as to continuously adjust the parameters of the model, in order to continuously adjust and improve the vibration reduction performance of track and field runway [14]. Yang and others have carried out risk assessment on the training process of badminton players, designed and constructed a platform of Internet of things to analyze the badminton process, and improved the badminton training skills according to the results of the study [15].

From the literature research related to badminton, racket design, and finite element analysis method, we can see that the current research on badminton is more from the analysis of movement and the physiological process of athletes; even the research on badminton is from the material design of badminton; the finite element analysis design is in the engineering design and manufacturing. It can achieve very good results. To sum up, this paper proposes to use the finite element analysis method to build the finite element model of badminton and analyze the mutual mechanics and displacement relationship between badminton and badminton racket based on the motion equation, so as to determine the relevant material parameters of badminton and carry out the relevant mechanical verification.

#### 3. Construction of Basic Finite Element Model of Badminton Batting Process Based on Motion Differential Equation

##### 3.1. Simplified Finite Element Model of Badminton Based on Motion Differential Equation

The batting moment of badminton is a very short process. The contact time between racket and ball is defined as batting time. The speed before the contact between racket and badminton is the initial speed. The speed after the contact is the return speed. The ratio of return speed and initial speed is the speed ratio of batting. In the construction of the finite element model, the stiffness of badminton, racket cable tension, and axial stiffness need to be considered. The vibration equation of the badminton can be simplified as one of the following equations:

In the badminton racket model, the stiffness perpendicular to the net is assumed to be nonlinear, so the displacement along the direction perpendicular to the net is also nonlinear. Let be three undetermined coefficients, the vertical displacement of the central point of the mesh is , and the vertical force on the central point of the mesh is ; then, the formula of vertical displacement in the mesh direction can be obtained:

In order to simplify the model, the impact point of the net surface is regarded as the centroid, and the net surface of the racket is a rigid body. The simplified mechanical model of batting can be obtained from Figure 1, as shown in Figure 2. The racket can be regarded as a combined system of linear spring and mass. In the process of badminton and spherical impact, it is a process of storing the hitting energy in the racket net, which will make the ball and badminton racket vibrate.

Using Newtonʼs second law, we can transform the badminton motion differential equation:

During the process of badminton catching, the impact of badminton on the racket is the resilience of the tennis surface of the racket; thus, we can get

Suppose that the velocity of the badminton player is , the displacement of the ball is , and the displacement of the contact point between the badminton and the racket is :

##### 3.2. Parameter Determination of Finite Element Model of Badminton and Racket

The mechanical properties of four common badminton rackets are shown in Table 1.

 Sample number Linear density (g/km) Specific strength (Cn  g/km−1) Initial modulus (Cn  g/km−1) Maximum strength (N) Elongation after fracture (%) Sample 1 441 51.8 173.5 228.4 43.1 Sample 2 472 51.8 112.1 244.5 45.5 Sample 3 527 36.1 115.6 190.0 45.6 Sample 4 531 24.0 114.9 127.3 46.6

In Table 1, the initial modulus of fiber refers to the linear stress-strain ratio of the initial part of the load elongation curve of the fiber used in badminton racket. In order to facilitate the calculation, four kinds of badminton racket samples are converted according to the unified diameter of 0.7 mm, from which the calculation parameters of the model can be obtained (see Table 2 for details).

 Sample number Mass density Youngʼs modulus Poissonʼs ratio Sample 1 1380 1.6e9 0.374 Sample 2 1370 1.6e9 0.376 Sample 3 1227 1.4e9 0.376 Sample 4 1146 1.9e9 0.378

Thirdly, the tension of badminton racket net is selected to represent the tightness of racket line. Different people have different preferences and applicability for the ball line weight. The specific corresponding relationship and use characteristics are shown in Table 3. Among them, professional players are more suitable for medium- and high-pound tension rackets. This type of racket has better ball control effect and can clearly feel the control of their own strength on the ball.

 Pound rating Pounds Tension of stay wire Use characteristics Low weight ≤20 ≤90 The racket line is very loose, and the ball has a significant residence time on the racket. It is suitable for pulling and hanging the ball. Medium and low pounds 21 94.5 Good elasticity, a little sense of retention, attack power will be weakened, suitable for amateurs 22 99.0 23 103.5 Medium pound 24 108.0 Amateurs use the most, can play a better performance, and has a significant high elasticity 25 112.5 Medium high pound 26 117.0 It is hard to pull the back court, and it is weak to hit the ball. It is very comfortable to control the ball. When the strength is large, it can have a very fast speed and can accurately feel the size of its own strength. 27 121.5 Big pound ≥28 ≥126 It is easy for professional athletes to control the ball and ensure the match between strength and speed

#### 4. Construction of Net Beam Element Model of Badminton Racket

##### 4.1. Construction of Differential Equations of Motion for Beam Element and Bar Element

In order to build the basic model of badminton racket beam element mesh, based on the previous relevant parameters of badminton, comparative experimental analysis of bar element and beam element is carried out. According to the definition of mechanics of materials, the relationship between force and displacement of a single string in a badminton racket is shown in the following equation:

In equation (6), is the load value of the middle point of the chord, and the value is 1n; is the cross-sectional area of the chord, and the value is 3.85  10−7 m2; is 1.4 GPa; is 1/2 chord length, and the value is 0.1 M; and is the displacement of the middle point of the badminton racket in the load displacement. The formula regards the stress of a single badminton racket as static load. On this basis, the beam element and rod element of the badminton racket can be selected for comparative analysis. The stress analysis is shown in Figure 3.

##### 4.2. The Finite Element Model Construction of the Net Beam Element and the Rod Element of Badminton Racket

The displacement of the nodes at both ends of the chord line of badminton racket is limited. Under the setting of each parameter in equation (6), 40 beam element finite element models are constructed and numerically calculated. The results of numerical calculation are compared with the analytical results calculated by the formula. The results are shown in Table 4.

 Load value (n) Analytical solution of displacement (mm) Numerical solution of displacement (mm) Error value (%) 0.0000 0.00 0.00 0.00 0.0156 3.07 2.81 8.47 0.0312 3.86 3.61 6.48 0.0469 4.42 4.20 4.98 0.0703 5.06 4.81 4.94 0.1055 5.81 5.52 4.99 0.1582 6.62 6.42 3.02 0.2373 7.61 7.41 2.63 0.3560 8.72 8.51 2.41 0.5339 9.96 9.78 1.81 0.8009 11.44 11.27 1.49 1.0000 12.31 12.18 1.06

Further analysis from Figure 4 shows that when the load value is 1n, the error between the analytical solution and the numerical solution is the smallest, which is 1.06%. The beam element does not change completely in the process of force, which shows that it is feasible to analyze the beam element when the chord line of badminton racket is thin.

The displacement of nodes at both ends of badminton chord is limited, and 1n load is applied at the midpoint of chord. Referring to the previous parameters, let be 1n, be 1.4 GPa, be 3.85 × 10−7 m−2, and be 0.1 m. When using the finite element method for analysis, considering the lack of a load balance force in the -direction, a constraint with minimal stiffness can be selected, and the error comparison of springs with different stiffness can be carried out according to the spring stiffness values of 0.1, 1, 2, 3, and 10. The results are shown in Figure 5.

It can be seen from Figure 5 that when the stiffness value is 2, the midpoint displacement value is close to the value of the analytical solution. It can be seen that the construction of the model can be applied to the numerical calculation of badminton net. Therefore, the finite element software is used to carry out the comparative analysis between the load displacement of the badminton club element and the analytical solution. The calculation results are shown in Table 5.

 Load value (n) Analytical solution of displacement (mm) Numerical solution of displacement (mm) Error value (%) 0.0 0.00 0.00 0.00 0.1 5.71 5.48 4.03 0.2 7.14 7.02 1.68 0.3 8.23 8.11 1.46 0.4 9.06 8.89 1.88 0.5 9.76 9.65 1.13 0.6 10.37 10.29 0.77 0.7 10.92 10.88 0.37 0.8 11.42 11.38 0.35 0.9 11.88 11.85 0.25 1.0 12.31 12.29 0.16

It can be seen from Figure 6 that the displacement numerical solution and analytical numerical solution of the badminton racket bar element are similar to the simulation results of the beam element. When the load value is 1n, the error value is the smallest.

From the simulation analysis of beam element and bar element, it can be seen that the error values of single-chord simulation results of badminton are in a small level, so it is necessary to compare and analyze the numerical calculation of badminton racket net surface of beam element and bar element.

#### 5. Model Construction and Numerical Results of Beam Element and Bar Element

##### 5.1. Static Numerical Calculation of Badminton Racket Net

The racket line of badminton is a state of overlapping. In the moment of hitting, some net lines on badminton racket will be separated, which makes the net line no longer interact effectively. In order to make the process analysis more close to the real situation, it is necessary to consider releasing some node constraints when building the finite element model. An elliptic finite element model with long axis and short axis of 28 cm and 22 cm is constructed. 430 intersections are set up. A position near the middle of racket is selected as the hitting point, then 19 nodes will be constrained to release, and the rest 411 intersections will have the force effect. In order to compare the energy value of beam element and rod element, some nodes are deleted to compare the displacement of badminton mesh. The maximum 50 N load is applied to the mesh surface. The displacement comparison of badminton net surface is shown in Table 6:

 Beam element Rod element Load value (n) Node not deleted (mm) Deleted node (mm) Load value (n) Node not deleted (mm) Deleted node (mm) 0.000 0.00 0.00 0.000 0.00 0.00 0.196 3.84 3.86 0.101 2.36 2.36 0.392 5.07 5.11 0.338 4.11 4.11 0.585 5.93 5.95 0.514 4.92 4.92 0.874 6.87 6.93 0.768 5.86 5.86 1.319 7.98 8.05 1.152 6.99 6.99 1.974 9.25 9.29 1.712 8.22 8.22 2.967 10.69 10.78 2.547 9.67 9.67 4.451 12.51 12.55 3.848 11.34 11.34 6.684 14.19 14.23 5.772 13.15 13.15 10.024 16.28 16.37 12.978 17.62 17.62 15.472 18.64 18.76 19.457 20.14 20.14 22.547 21.44 21.49 29.214 23.45 23.45 33.724 24.57 24.63 43.175 26.93 26.93 50.000 28.11 28.15 50.000 28.17 28.17

The displacement curve of beam element and rod element with some nodes deleted and nodes not deleted are drawn respectively, as shown in Figure 7.

In Figure 7, A and B, respectively, represent the displacement curves of some nodes not deleted or deleted in the beam element; C and D, respectively, represent the displacement curves of some nodes not deleted or deleted in the rod element. Through the comparative analysis of Figure 7, it can be seen that after deleting some nodes, the net stiffness of badminton racket is reduced, and the displacement change of chord midpoint is the largest under the same load. Considering that in the moment of badminton hitting, some nodes will not be forced, so in this simulation analysis, the deletion of nodes is more in line with the actual situation. Further comparison of the displacement curve difference after deleting the node shows that the difference between the midpoint displacement of the beam element and the bar element is 0.28%, and the difference between the midpoint displacement after deleting the node is 0.43%.

##### 5.2. Simulation Results of Net Tension of Badminton Beam Element

The string tension is added to the net surface of badminton racket, so we can get the formula of analytical solution as follows:

In the finite element model analysis, the applied load at the middle point is 1n, the chord length is 0.2 m, and the value of in equation (7) is 539n. After fixing the two ends of the chord, observe the displacement curve of the beam element and the bar element model under different tension of the stay wire. The results are shown in Table 7.

 Tension of stay wire (N) Analytical solution of displacement (mm) Displacement of beam element (mm) Error value (%) Displacement of rod element (mm) Error value (%) 0.108 12.14 12.06 0.66 12.01 1.07 0.524 11.74 11.68 0.51 11.63 0.94 1.054 11.18 11.14 0.36 11.04 1.25 2.134 10.15 10.12 0.30 10.02 1.28 5.347 7.42 7.35 0.94 7.29 1.75 10.248 4.48 4.47 0.22 4.31 3.79

The curve of the simulation calculation results of the pull wire tension of badminton beam element is shown in Figure 8. A is the analytical solution of displacement, B and C are the displacement value and error value of beam element, respectively, and D and E are, respectively, the displacement value and error value of the rod element. According to the displacement curve of the pull tension, it can be seen that with the increasing tension of the cable applied, the displacement value of the midpoint of the chord line on the badminton racquet net surface decreases gradually, which indicates that the displacement of the action point of the mesh surface will be reduced gradually under the same position and the same load. That is, when a constant force is applied to the chord, the completeness near the midpoint position is the lowest.

The error of beam element is the smallest under different tension, so it is more reasonable to use beam element for numerical simulation.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.

#### References

1. R. Abedanzadeh, K. Becker, and S. Mousavi, “Both a holistic and external focus of attention enhance the learning of a badminton short serve,” Psychological Research, vol. 17, pp. 1–9, 2021. View at: Publisher Site | Google Scholar
2. M. Naoto and K. Koji, “Comparison of modular control during smash shot between advanced and beginner badminton players,” Applied Bionics & Biomechanics, vol. 2018, Article ID 6592357, 6 pages, 2018. View at: Publisher Site | Google Scholar
3. R. W. L. Lee, C. Y. Andy, and T. W. L. Wong, “Application of analogy in learning badminton among older adults: Implications for rehabilitation,” Motor Control, vol. 23, no. 3, pp. 384–397, 2019. View at: Google Scholar
4. M. Á. Gómez, F. Rivas, A. S. Leicht, and J. M. Buldú, “Using network science to unveil badminton performance patterns,” Chaos, Solitons & Fractals, vol. 135, Article ID 109834, 2020. View at: Publisher Site | Google Scholar
5. L. K. Huang, L. N. Huang, Y. M. Gao et al., “Electrical impedance myography applied to monitoring of muscle fatigue during dynamic contractions,” IEEE Access, vol. 99, p. 1, 2020. View at: Google Scholar
6. M. F. Che Hassan, M. U. Rosli, and M. A. Mohd Redzuan, “Material selection in a sustainable manufacturing practice of a badminton racket frame using Elimination and Choice Expressing Reality (ELECTRE) Method,” Journal of Physics: Conference Series, vol. 1020, Article ID 012012, 2018. View at: Publisher Site | Google Scholar
7. S. Mori, K. Tanaka, S. Nishikawa, R. Niiyama, and Y. Kuniyoshi, “High-speed and lightweight humanoid robot arm for a skillful badminton robot,” IEEE Robotics & Automation Letters, vol. 3, no. 3, pp. 1727–1734, 2018. View at: Publisher Site | Google Scholar
8. A. H. Bhelsaikar, D. Roy, and V. Atpadkar, “Estimation of dynamic characteristics of a novel non-parallel detachable-jaw robotic gripper using finite element method,” Journal of Physics: Conference Series, vol. 1831, no. 1, Article ID 012010, 2021. View at: Publisher Site | Google Scholar
9. S. Badia, F. Verdugo, and A. F. Martín, “The aggregated unfitted finite element method for elliptic problems,” Computer Methods in Applied Mechanics and Engineering, vol. 336, pp. 533–553, 2018. View at: Publisher Site | Google Scholar
10. R. Hiptmair, L. Li, S. Mao et al., “A fully divergence-free finite element method for magnetohydrodynamic equations,” Mathematical Models and Methods in Applied Sciences, vol. 28, no. 4, pp. 659–695, 2018. View at: Google Scholar
11. Z. Wang, C. Jiang, B. Ni, C. Wang, J. Zhong, and T. Fang, “An interval finite element method for electromagnetic problems with spatially uncertain parameters,” Science China Technological Sciences, vol. 63, no. 1, pp. 25–43, 2020. View at: Publisher Site | Google Scholar
12. C. Itu, A. Öchsner, S. Vlase, and M. I. Marin, “Improved rigidity of composite circular plates through radial ribs,” Proceedings of the Institution of Mechanical Engineers, Part L: Journal of Materials: Design and Applications, vol. 233, no. 8, pp. 1585–1593, 2019. View at: Publisher Site | Google Scholar
13. S. Karthikeyan and S. Chandru, “Linear modal analysis of ring shaped-elliptical plate using finite element method Linear modal analysis of ring shaped-elliptical plate using finite element method,” Journal of Physics Conference Series, vol. 1597, no. 1, Article ID 12041, 2020. View at: Google Scholar
14. L. Andena, S. Aleo, F. Caimmi et al., “Modelling the cushioning properties of athletic tracks,” Sports Engineering, vol. 21, no. 4, pp. 453–463, 2018. View at: Publisher Site | Google Scholar
15. J. Yang, H. Xia, Y. Wang, and H. Tian, “Simulation of badminton sports injury prediction based on the internet of things and wireless sensors,” Microprocessors and Microsystems, vol. 81, no. 1, Article ID 103676, 2021. View at: Publisher Site | Google Scholar

Copyright © 2021 Yingda Li. 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.