Research Article  Open Access
Jingzhao Yang, Jiuzhi Dong, Xiuming Jiang, Yunjun Chen, "Reducing the Collision Damage Done to the Tips of Steel Needles during Integrated Piercing by Using Shape Optimization with Feature Selection", Mathematical Problems in Engineering, vol. 2020, Article ID 5078590, 19 pages, 2020. https://doi.org/10.1155/2020/5078590
Reducing the Collision Damage Done to the Tips of Steel Needles during Integrated Piercing by Using Shape Optimization with Feature Selection
Abstract
Recent research has shown that, during integrated piercing, the yarn tension can displace the needles from the centers of the holes in the piercing template. To reduce the damage done to the needle tips when the needles and the hole walls collide when the needle displacement is too large, this paper proposes a method for optimizing the needle shape that satisfies the strength constraint while targeting minimum needle displacement. First, the optimal objective function and strength constraint conditions for the tip displacement of the integrated puncture needle are established, which are affected by many factors. Then, the feature selection method of machine learning is used to reduce the dimensionality of the objective function after data reduction, and the feasible region of key features is reduced to avoid falling into the local best point in the optimization iteration. Finally, machine vision is used to measure experimentally the tip displacement of a needle array. The results show that the average tip displacement is reduced by 18.16–31.42% and the collision risk is reduced from 43.14% to 5.54%. It shows that the optimization method of needle shape based on feature selection is effective for reducing needle tip impact damage during integrated piercing.
1. Introduction
Carbon/carbon (C/C) composite preforms have received much attention in recent years because of their excellent properties of high specific strength, high specific modulus, and good ablation resistance. Integrated piercing (IP) on 3D fabrics has good overall structure and fiber volume fraction; they are excellent substrates for highperformance and heatinsulating C/C composite materials and are used widely in aircraft brake pads and rocket nozzles after curing C/C composites [1–5]. The distance between the needle tip of the steel needle array and the hole center of the template is determined by the structure and geometric parameters of the steel needle array [6–11]. The distance directly affects the uniformity of fiber volume fraction and the distribution of the damage degree between the tip and the hole wall.
Truss structure optimization with application diversity is one of the most popular problems in civil engineering [12–19]. Truss optimization includes three types: section shape, geometry, and topology. The integrated piercing process of woven carbon cloth is affected by many factors, such as fabric structure, needle geometry parameters, array layout, and template parameters. As a special application form of truss optimization in the field of textile machinery, the shape optimization of steel needle array aiming at reducing the impact damage of needle tip is a highly nonlinear optimization problem.
Optimization methods can be generally divided into two categories: (1) gradient descent method represented by gradient derivation and (2) metaheuristic algorithm combining random method and local search. Gradient descent method is the most popular machine learning optimization algorithm in recent years. However, choosing the appropriate learning rate when using the gradient descent method to avoid falling into the local best point is a problem. Improvements made for the application of gradient descent include minibatch gradient descent, annealing gradient descent, and accelerated stochastic gradient descent [20–22]. On the contrary, research studies based on natural heuristic technologies, such as simulated annealing, genetic algorithm, tabu search, ant colony optimization, particle swarm optimization, cuckoo search, leapfrog, flower pollination algorithm, fruit fly optimization algorithm, and ideal gas molecular movement, have made important progress. Researchers have applied and improved these algorithms in the field of truss structure optimization [12, 14, 19, 23]. However, there is currently no specific metaheuristic algorithm that can be applied to the optimal solution of all structures and geometric dimensions. Proper selection of algorithms and parameters for a specific optimization problem is a timeconsuming process, and engineering applications need to combine experience and trial and error. Naderi et al. [23] proposed a multistage metaheuristic method called topup technology, which can solve the problem that metaheuristic algorithm falls into local optimum. In this paper, a feature selection method is proposed to reduce the feasible region of optimization problems under constraints. It can avoid the metaheuristic algorithm falling into the local optimal solution range of nonkey factors in the iterative solution process.
Neither the gradient descent method nor the element heuristic algorithm can be used to solve the problem well. In this paper, a singleneedle displacement model satisfying the strength requirement is proposed, which takes into account the fabric structure parameters, needle shape parameters, and punching template parameters. The paper then simplifies the model and retains the parameters that are highly correlated with the target feature parameters, and it optimizes the calculation examples given within the allowable range of process parameters. Finally, the results of this optimization are applied to an experiment, and data samples of the needle tip displacement of a steel needle array are obtained using machine vision. The mean and variance of the needle displacement in the needle array from the corresponding hole center of the piercing template are analyzed to test the accuracy and applicability of the theoretical model.
2. Displacement Modeling Method of Needle for Integrated Piercing
During the IP of woven carbon cloth, the displacement of the steel needles away from the hole centers of the piercing template should be reduced as much as possible under the condition of meeting the constraint of the needle material strength. First, in modeling the yarn tension, three factors are considered, namely, (i) the horizontal needle winding, (ii) the Zdirection slippage, and (iii) the compensation of yarn tension caused by the Zdirection bending in different fabric structures. Second, the constraint condition between the parameters of the steel needles and the tension is established on the premise that the influence of the number of carbon layers on the steel needles is ignored. Finally, the deflection model of a combined variablesection cantilever beam under concentrated load is used to simplify the calculation of the displacement of a single IP steel needle under the tension of the yarn.
2.1. Tension Modeling of Representative UnitCell Yarn for Integrated Piercing of Woven Carbon Cloth
During IP, the needle array disturbs the yarn pattern in the carbon cloth and squeezes into the fiber gaps, causing the fiber to stretch horizontally and vertically, as described by Zhu [4]. The fiber tension acts on the steel needles to cause horizontal displacement, and if that displacement is too large, then the needles will collide with the hole walls and damage the needle tips. To avoid that problem, it is very important to study the representative unitcell (RPC) yarn tension model of the interaction between the equidistant and closely arranged steel needles and the yarn.
2.1.1. Analysis of Integrated Piercing of Woven Carbon Cloth
The IP process can be summarized as (i) arrangement of needle array, (ii) integral puncture, (iii) pressurization compaction, and (iv) needle replacement [24, 25]. The detailed process of IP preform is shown in Figure 1.
(a)
(b)
(c)
(d)
(e)
(f)
Each needle in the array corresponds to a pinhole in the piercing template. Then, when the piercing template moves downward, the equidistant dense steel needles touch the carbon cloth, as shown in Figure 2. The IP equipment has two sets of orthogonal direction carbon cloth conveying devices, which work alternately in sequence [26].
As shown in Figure 3, the IP process has four stages. Figure 3(a) shows the woven carbon cloth is placed on the top of the needle array; at this time, the carbon cloth experiences zero net force and remains stationary. Figure 3(b) shows the fiber straightening caused by the carbon cloth being pressed onto the needle tips when the piercing template is moved slowly downward to touch the static carbon cloth. Figure 3(c) shows that, as the piercing plate continues to move downward, the needle tips penetrate the carbon cloth and occupy the spaces therein. In short, the area of the carbon cloth in which the steel needles are inserted changes gradually from the smallest (the needle tips) to the largest (the steel needle rods) and then remains constant. The two processes described above are shown in Figures 3(c) and 3(d).
(a)
(b)
(c)
(d)
There are two reasons for the complex and irregular bending and elongation of most fibers when they enter the gaps of the needle array. First, in the ideal state, the fibers in the warp and weft directions are evenly distributed in the gaps of the steel needle array, but it is difficult to achieve in practice. Second, when the manipulator holds the woven carbon cloth, the consequent wrinkling and deformation prevent the warp and weft yarns from remaining absolutely orthogonal, which is determined by the loose structure of the woven carbon cloth.
The fiber displacements occur during IP include horizontal displacement (Figure 4(a)) and vertical displacement (Figure 4(b)). The horizontal displacement is caused by the needle tips squeezing the fiber, the vertical displacement is caused by the frictional force of the fiber, and Figure 4(c) shows a plan view of the simultaneous horizontal and Zdirection displacements. Although in some ways similar to [5], the present study refines the model and considers the fiber slip in the vertical direction.
(a)
(b)
(c)
2.1.2. Bending Model of Fiber around Needle considering ZDirection Slip
(1) Maximum Elongation of Bending Fibers around Adjacent Steel Needles. In the second case of fiber winding around the needle, as shown in Figure 4(a), the fiber has the largest elongation. As shown in Figure 4, the yarn paths between the three steel needles i_{21}, i_{22}, and i_{23} have symmetry; therefore, only the length of the fiber B1B2 between i_{21} and i_{22} is analyzed in Figure 5, and the dotted line between C_{2}B_{2}′ is the projection of C_{2}B_{2} on the plane.
(a)
(b)
The fiber length B_{1}B_{2} after bending and stretching between two adjacent needles can be obtained from Figure 5 as follows:
The foursegment lengths that makeup B_{1}B_{2} are described by equations (2)–(5), and these formulas are derived from the geometric relationship in Figure 5:where is the tip radius of the contact point between the yarn and tip 1, R_{2} is the same for tip 2, and d_{kj} is the distance between the two holes adjacent to the center of the piercing template.where d_{k} is the hole diameter.where h is the distance from the yarn to the piercing template at any time.
(2) Compensation for Bending and Stretching of Adjacent Steel Needle Fibers around the Needle. The close contact of the needle array with the carbon cloth makes it impossible to obtain compensation from the fiber segment at a far distance when the fiber is bent and stretched. Therefore, fiber compensation is considered to occur only in an RPC of the fabric structure in this study. In theory, plainweave, twillweave, or satinweave carbon cloth can be used as an IP woven carbon cloth. However, considering the high modulus and brittleness of the carbon fiber, to reduce the fiber abrasion during weaving of the woven carbon cloth, it is preferable to use a satinstructured carbon cloth, as suggested in [4]. In the present paper, 8/3 satin fabrics are taken as an example, and their organizational structure is shown in Figure 6. The curve of the fabric structure phase is expressed by a sinusoidal curve with some adjustments in [27].
In Figure 6, the length s of the bent fiber can be obtained by integrating the fiber path curve at (0, L_{xj}):where Z is the fiber path curve, C is the crimp height of the fabric structure phase, L_{xj} is the spacing of adjacent warp or weft yarns in the fabric, and τ is the correction coefficient of the yarn spacing.
2.1.3. Fiber Tension between Adjacent Steel Needles
From fiber material science, we havewhere E_{c} is the carbon cloth modulus of elasticity and A_{z} is the carbon cloth fiberbundle crosssectional area. Satinweave carbon cloth is soft and deforms easily, and the crimp height during the orthogonal stacking process is affected by engineering factors, so the adjustment factor τ should be considered. During the IP, the needle tip only directly contacts the adjacent fibers, so the adjustment coefficient η related to the yarn k number, the fabric structure, and the needlepoint diameter should be introduced. Therefore, from equations (1)–(9), the tension model during the IP process can be obtained as follows:where d_{ned} is the needle diameter.
From the knowledge about fiber materials, it is known that the yarn elongation determines the yarn tension received by the needles. However, equations (1)–(9) show that a number of factors determine the tension of the yarn during IP. To ensure needle stability, the tension is calculated according to the maximum possible elongation. To simplify the calculation, R_{1}, R_{2}, and d_{ned} are made equal herein, and thus, the value of h in equation (5) is actually a multivariable numerical value that is very difficult to calculate accurately because of the interaction between the piercing template and the needles in a certain range. Based on the engineering experience in [4], h = 0.5 mm can be obtained.
By substituting the parameters of Table 1 into Equation (10), the yarn tension F_{max} was found to be 14.69 N, from which the yarn strain ε was calculated easily as 1.67 × 10^{−2} based on the values of the tension F_{max}, the elastic modulus E_{c}, and the crosssectional area A_{z}. It is known from engineering experience that yarn with a strain of less than 2 × 10^{−2} is stretched safely.

2.2. Displacement Modeling of Steel Needle under Yarn Tension
To study how the yarn tension influences the needle displacement in IP, the deflection model of a composite variable crosssectional cantilever beam under variable concentrated load was used to study how the needle displacement changes under the influence of the yarn resultant force. The displacement modeling of a single steel needle under the action of yarn tension is divided into three steps, namely, (i) force analysis of the needle, (ii) strength constraints, and (iii) displacement modeling.
2.2.1. Force Analysis of Steel Needle under Yarn Tension
For the convenience of analysis, the IP needle force model is rotated 90° to the right, as shown in Figure 7. The coordinate system is established with the needle tip O as the origin, horizontally left as the X direction, and vertically upward as the Y direction. Then, the crosssectional radius y of the needle tip section changes with x, and the curve of the needle tip is given an exponential form [4], which considers the processing feasibility.
(a)
(b)
(c)
Figure 7(a) shows that the needle tip curve passes through point B′ with coordinates (L_{1}, d_{ned}), from which we can derivewhere L_{1} is the length of the needle tip, d_{ned} is the radius of the needle rod, m is the exponent in the power function representing the tip shape curve, and x is the distance from any point to the apex of the tip.
For Figure 7(b), when the yarn acts at any point C′ of the tip, the direction of F is the normal direction of the needle tip curve of the contact point, and F is decomposed into F_{1} in the horizontal direction and F_{2} in the vertical direction. As shown in Figure 7(c), the F_{2} action point is translated to the intersection point C, which is the intersection point of the neutral layer of the steel needle and F_{2}, and then, F_{1} can be translated to the tip vertex F_{1}′. Next, a couple of M1 must be added, which is determined by the translation theorem of the force, and the couple moment of component F_{2} to point A is M_{2}.
The position of the contact between the fiber and the IP needle changes with the IP process, and the direction of the force F is always perpendicular to the tangent line of the contact point C′. As shown in Figure 7(b), if the angle between F and F_{2} is θ, then
From Figure 7(b), the tension components F_{1} and F_{2} are calculated as
Also, from Figure 7(b), the relationship between F1 and F2 is
2.2.2. Strength Constraint Condition When Bottom of Steel Needle Is Fixed
According to the force analysis of Figure 7(c), the cross section of A′A″ is a dangerous section. To analyze the position of the dangerous point in the section A′A″, the stress distributions of (a) the force component F_{2}, (b) the additional force moment, and (c) the axial force component F_{1} are plotted in Figure 8.
(a)
(b)
(c)
Figure 8 shows that point A′ and point A″ are the dangerous points of the needle, and the maximum stress condition at A′ meeting the strength requirement iswhere A is the crosssectional area of the contact point between the yarn and the needle tip, σ is the allowable stress of the steel needle, and n is a safety factor. Combining equations (11) and (16) gives
Similarly, the maximum stress condition at A″ satisfying the strength requirements is
Combining equations (11) and (18) gives
2.2.3. Modeling of Needle Displacement during Integrated Piercing
The needle displacement from the center of the corresponding hole in the piercing plate is affected by many factors, such as the fabric structural parameters, the needle shape parameters, and the orifice parameters. Herein, a model for calculating the IP needle displacement under carbon fiber yarn tension is converted into one for calculating the deflection of the combined variable crosssectional cantilever beam under concentrated load. The authors in [28–30] consider how to calculate the deflection of various variablesection cantilever beams under different loads. Unlike [8], the present paper considers the deflection caused by the cumulative cornering of segmental steel.
Figure 9 shows the calculation and analysis of the deflection of the cantilever beam of the piercing steel needle. Figure 9(a) shows the needle force as the fiber bends around the needle, in which AB is a constant crosssectional needle bar, BO is a variable crosssectional needle tip, and F_{2} acts at any point C. Figure 9(b) shows the additional force moment required for F_{2} to translate to point B, which is determined by the principle of force translation, namely,
(a)
(b)
(c)
(d)
(e)
(f)
Calculating the forced needle displacement is equivalent to calculating the deflection of the combined variablesection cantilever beam. For small deformation, it is analyzed in a stepwise manner, as shown in Figures 9(c)–9(f). The deflection of a composite variable crosssectional cantilever beam has four parts, namely,where is the deflection under the concentrated load F_{2} after the tempered section AB (Figure 9(c)), is the deflection of the concentrated load F_{2} after the tempered segment BO is translated to point B and F_{2}′ (Figure 9(d)), is deflection due to additional moment M of force F_{2} after translation to point B (Figure 9(e)), and is the deflection of the O point caused by the rotation at the connection point B (Figure 9(f)).
Beam deflection can be calculated by the integral method, the superposition method, the graph multiplication method, or the Mohr integral method, among others [31, 32]. The double integral method and the singularpoint function method are mainly used to calculate the complete deflection curve of the beam. However, only the displacement of the maximum deflection point of the cantilever beam is calculated herein, and the maximum displacement of the curved cantilever beam can be calculated by the method of moment area [31, 33]. Figure 10 shows the piecewise deflection analysis.
(a)
(b)
(c)
(1) Displacement v_{1}and Rotation Angle θ_{1}. The action of the fiber on the tip of the variable cross section is equivalent to F_{2} acting on any point C′ of the variable crosssectional cantilever beam. In this case, the method of using the moment area is complicated, and the midsection moment of inertia of the straightline cantilever beam after the second integral can be replaced by I_{x} for analysis, expressed as
The maximum angle at point O of the cantilever beam in Figure 9(c) is
The maximum displacement at point O of the cantilever beam in Figure 9(c) is
(2) Displacement v_{2}and Rotation Angle θ_{2}. For Figure 9(d) under small deformation, the moment area method giveswhere E is the modulus of elasticity of the steel needle, and
(3) Displacement v_{3}and Rotation Angle θ_{3}. For Figure 9(e) under small deformation, this study selects the free end B of the cantilever beam as the origin and substitutes M_{x}, I, and AB to obtain
(4) Calculation of Cantilever Beam Deflection v. The displacement of L_{1} in Figure 9(f) is due to the rotation angles of Figures 9(d) and 9(e). Therefore,
The calculated maximum needle deflection obtained by substituting equations (24), (26), (28), and (29) into equation (25) is
Then, substituting equation (14) into (30) gives
3. Optimization of Needle Shape Parameters Based on Selection of Key Features
The main difficulty in optimizing the needle shape is that the model is multiobjective and multivariable. In machine learning, variable parameters are known as features, and a very important issue is determining how each feature influences the optimization objective. Herein, the machine learning method of feature selection is used to evaluate and sort the correlation between features and the optimization objective, and these features with large correlation are selected for example optimization. Then, the optimal needle shape parameters that meet the strength and displacement conditions are determined by deflection calculation and crossvalidation.
Al Nuaimi et al. [34] review the traditional feature selection algorithms and then scrutinize the current algorithms that use streaming feature selection to determine their strengths and weaknesses. The authors in [35] provide the basic concepts necessary to build an ensemble for feature selection, as well as reviewing the uptodate advances and commenting on future trends. The authors in [31] provide an overview of feature selection techniques and the instability of the feature selection algorithm.
3.1. Selection of Key Characteristic Variables for Optimization Problems
Feature selection usually considers (i) whether a feature is divergent and (ii) the relationship between a feature and the target, and those features with high relevance to the target should be selected. Feature selection comes in three forms, namely, filter, wrapper, or embedded. The common feature selection algorithms are compared in Table 2, showing that the Pearson coefficient method is suitable for continuous, linear, and normally distributed data, which is not the case in the present study. Mutual information can be regarded as the uncertainty of one random variable being reduced because another random variable is known [32]. Proposed in 2011, the maximum information coefficient (MIC) [36] is the latest method for detecting nonlinear correlations between variables. Applying the concepts of information theory and probability to continuous data, the MIC algorithm can represent various linear and nonlinear relationships and has been used widely. Its range is between 0 and 1, and the higher the value, the stronger the correlation.

To study the relationship between needle displacement and the various influencing factors, it is necessary to provide batch sample data for the feature selection algorithm. The range and iteration step of each variable in equation (31) are given in Table 3.

Set F_{max} = 15 N and E = 2.11 × 10^{11} Pa and use Maple 2019 to calculate equation (31) numerically to obtain 97,510 data. The MIC calculation process discussed in [32] can be described aswhere (x, y) represents the sample value and (X, Y) represents the sample set. The results calculated using Python 3.7.3 and minepy 1.2.3 are given in Table 4, where MIC is the maximum information coefficient, MAS is the maximum asymmetry score, MEV is the maximum edge value, GMIC is the generalized MIC, and TIC is the total information coefficient.

To analyze the positive and negative effects of the parameter characteristics on the objective function, the Kendall coefficient is used. The results calculated by using Python 3.7.4 and scikitlearn 0.21.3 are given in Table 5.

From Tables 4 and 5, the order of influence of each parameter on the needle displacement is d_{ned} > L_{2} > m > L_{1} > x, where m, d_{ned}, and x are negative correlations and L_{1} and L_{2} are positive correlations.
3.2. Optimization Example of Needle Diameter and Bar Length
Based on the conclusion reached in Section 3.1, this study optimizes the two most important characteristics that affect needle displacement, namely, d_{ned} and L_{2}.
3.2.1. Constraint Analysis of Displacement Conditions on Key Morphological Parameters
The displacement constraint is the result of equation (31) must be less than to ensure that the needle tip does not collide with the hole wall, where β is the safety factor.
(1) Needle Diameter. Given the process conditions of L_{2} = 7 mm, d_{k} = 2.3 mm, and β = 0.7, constraint conditions become as follows:
For F_{max} = 0 to 50 N and d_{ned} = 0.8 to 2.0 mm, the contour relationship among , F_{max}, and d_{ned} is as shown in Figure 11.
The constraint condition in equation (33) shows that the region above and to the left of the 8.05 × 10^{−4} displacement contour in Figure 11 is safe and feasible. From equation (33), this can be obtained by solving the boundary condition 0.805 mm. For F_{max} ≤23 N, d_{ned} satisfies the displacement safety requirements of the current example in the range of 0.8 to 2.0 mm. For F_{max} = 50 N, d_{ned} ≥1.035 mm meets the safety requirements of the current example of displacement.
(2) Needle Rod Length. Given the engineering conditions of d_{ned} = 1.1 mm, d_{k} = 2.3 mm, and β = 0.7, equation (31) becomes as follows:
For F_{max} = 0 to 50 N and L_{2} = 5 to 20 mm, the contour relationship among , F_{max}, and L_{2} is as shown in Figure 12.
The constraint conditions of equation (35) show that the area below and to the left of the 8.05 × 10^{−4} m displacement contour in Figure 12 is safe and feasible. Equation (34) can be obtained by solving the constraint boundary condition 8.05 × 10^{−4} m. For L_{2} = 20 mm, F_{max} meets the displacement safety requirements of the current example within the range of 0 to 11.9 N. For F_{max} = 50 N, L_{2} ≤7.95 mm meets the displacement safety requirements of the current example.
3.2.2. Analysis of the Strength Constraint of the Needle under Stress
(1) Needle Diameter. Given that the commonly used conditions of the project are L_{1} = 5 mm, L_{2} = 7 mm, σ = 1.514 × 10^{8} Pa, and a safety factor of n = 2.5, equations (17)–(19) become as follows:
For F_{max} = 0 to 50 N and d_{ned} = 0.8 to 2.0 mm, the contour relationships among σ, F_{max}, and d_{ned} are as shown in Figure 13.
(a)
(b)
(c)
(d)
Figures 13(a) and 13(b) show that the areas (i) above and to the right of the 6.056 × 10^{7} strength contour at point A′ for 0 to 10 N and (ii) above and to the left of the −6.056 × 10^{7} Pa strength contour at point A′ for 10 to 50 N constitute the allowable safety range of needle strength. Figures 13(c) and 13(d) show that the areas (i) above and to the right of the −6.056 × 10^{7} strength contour at point A″ for 0–10 N and (ii) above and to the left of the 6.056 × 10^{7} Pa strength contour at point A″ for 10 to 50 N constitute the allowable safety range of needle strength. The results obtained by solving equation (36) subject to the boundary conditions of 6.056 × 10^{7} Pa and −6.056 × 10^{7} Pa are given in Table 6.

(2) Needle Rod Length. Given that the commonly used conditions of the project are L_{1} = 5 mm, d_{ned} = 1.1 mm, σ = 1.514 × 10^{8} Pa, and a safety factor of n = 2.5, equations (17) and (19) become
For F_{max} = 0 to 50 N and L_{2} = 5 to 20 mm, the contour relationships among σ, F_{max}, and L_{2} are as shown in Figure 14.
(a)
(b)
(c)
(d)
Figure 14(a) shows that the area below and to the right of the 6.056 × 10^{7} strength contour at point A′ for 0 to 10 N is the safe range allowed by the needle strength. However, the −6.056 × 10^{7} strength contour at point A′ does not appear in the range of 10 to 50 N in Figure 14(b), so all that can be inferred is that most of the areas in Figure 15(b) satisfy the requirements of the needle strength. Figure 14(c) shows that the area below and to the right of the 6.056 × 10^{7} strength contour at point A″ in the range of 0 to 10 N is the safe range allowed by the needle strength. Similar to Figure 13(b), the 6.056 × 10^{7} strength contour at point A″ does not appear in the range of 10 to 50 N in Figure 14(d). Figures 14(b) and 14(d) show that the solution corresponding to the boundary condition at this time is not within the scope of the previous domain. The results obtained by solving equation (37) subject to the boundary conditions of 6.056 × 10^{7} and −6.056 × 10^{7} are given in Table 7.
(a)
(b)
(c)

3.2.3. Optimization Results of d_{ned} and L_{2} Example
The paper gives the range of parameters allowed by the process: F_{max} = 0 to 50 N, d_{ned} = 0.8 to 20 mm, L_{2} = 5 to 20 mm, and the feasible domain that satisfies the displacement and strength constraint mentioned above is d_{ned} = 1.159 to 2.0 mm and L_{2} = 7.07 to 7.96 mm.
The calculation parameters in Table 8 are substituted into equation (31) to obtain the partial derivatives of d_{ned} and L_{2} as follows:

Equation (37) shows that displacement is decremented monotonically with d_{ned}, and increases monotonically with L_{2}. However, the thicker the needle bar is, the more the fiber is squeezed, which reduces the fiber volume fraction of the 3D fabric, and the shorter the needle bar is, the fewer the carbon cloth layers that can be pierced at once, which reduces the work efficiency. Therefore, the optimization results of this example are d_{ned} = 1.159 × 10^{−3} m and L_{2} = 7.95 × 10^{−3} m.
Obtained numerically with four groups of parameters, the results given in Table 9 show that the theoretical displacement of this example varies from 0.206 to 0.45 mm, and the effect of parameter optimization on reducing displacement is obvious.

4. Experiment and Results
The randomness of the yarn around the needle bending direction means that the direction of the resultant force of the needle tip during IP is unknown. Therefore, it is difficult to compare directly the calculation results of the needle displacement in the theoretical model with the displacement of each steel needle in the experiment. This paper designs a verification scheme based on machine vision to measure tip displacement. Four sets of results were obtained using a combination of two needle diameters (1.0 mm and 1.159 mm) and two puncture needle lengths (7.95 mm and 10 mm). After image processing, four sample data sets for tip displacement were obtained, and the mean and variance were calculated to assess how optimizing the needle shape parameters affected the displacement of the needle array.
4.1. Experimental Parameters
The system has three main parts, namely, (i) integral puncture tooling, (ii) an industrial camera, and (iii) a light source. The tooling and material parameters of the experiment are given in Table 10, and the measurement requirements and camera parameters are given in Table 11.


4.2. Experimental Operation
In the experiment, the tip states of four groups of parameters were photographed. Two photographs were taken for each group of states to reduce the error caused by lighting and equipment stability. Figure 15(a) shows the machine vision imageacquisition experimental platform. Figure 15(b) shows four situations, namely, (i) a damaged needle, (ii) an intact needle, (iii) a dangerous position, and (iv) a safe position. Figure 15(c) shows the designed experimental process.
4.3. Image Processing
Computer vision technology is used widely in various fields, such as industrial automation, inspection and monitoring, visual navigation, humancomputer interaction, and virtual reality. The application of image processing technology to the positioning of steel needle tips generally requires three main stages, namely, (i) preprocessing, (ii) image segmentation, and (iii) feature extraction. As shown in Figure 15(b), the unsafe and dangerous positions are divided according to the distance of the needle tip from the hole wall. An example of the detailed image processing algorithmic flow and image processing is shown in Figure 16.
(a)
(b)
The data acquired in the actual project will inevitably contain missing values, abnormal points, and noisy data. The data are cleaned as follows: first, this paper observed the basic distribution of the data through histograms, density maps, and box plots. Then, because the missing data in this data set accounted for only 0.05% in the statistical analysis of the sample set, the direct deletion was used for simple processing. Finally, a box chart was used to analyze the outliers in the sample data, and these were replaced by the average value calculated from the data around the outliers.
4.4. Analysis of Results
As mentioned earlier, the present research aims to reduce the collision damage caused by the excessive displacement of the steel needles in contact with the hole walls during the overall puncture process. Table 12 shows that the recognition rate for the needle tip displacement is above 99%, and the experimental design is safe and reliable; at the same time, the proportion of needle tip displacement danger areas has been reduced from 43.13% to 5.54%. The displacement safety rate and dangerous rate indicate that the needle tip displacement safety rate is higher when the needle is wider and the needle bar is shorter and is lower when the needle is narrower and the needle bar is longer.

The analytical results of the data in Table 13 show that the mean and variance of needle array displacement are smaller when the needle tip is wider and the needle bar is shorter, and the mean and variance of needle array displacement are larger when the needle tip is narrower and the needle bar is longer. Compared with the fourth group of experiments, in the other three groups with optimized parameters, the average tip displacement of the needle array is reduced by 18.16–31.42%, and the variance is reduced by 22.51–27.71%. Comparing the results of the two experiments in each group, the reliability of the experimental sample collection can be seen from Tables 12 and 13.

Figure 17 shows the distributions of needle tip displacement under all four sets of parameters. Figures 17(a) and 17(b) show the average and variance of the steel needle array with small displacements, but the curve fluctuates more obviously. The effect of the thread on the complex tension of the steel needle is relatively obvious. Figures 17(c) and 17(d) show the large displacement mean and variance of the steel needle array, and the curve fluctuations have weakened significantly. The distribution law in Figure 17(d) is close to a normal distribution, which indicates that the effect of fiber tension decreases when the tip displacement is large.
(a)
(b)
(c)
(d)
The statistics and analysis results of the image processing results are consistent with the laws of the theoretical model, and although they deviate from the calculations of the theoretical model, this is acceptable. During the experiment, the displacement unit used pixels without camera calibration and distance conversion. This is because the purpose of this experiment is not to calculate the actual physical displacement of a single steel needle but to verify the correctness of the laws in the theoretical model based on statistical analysis. Without coordinate transformation, the experimental steps are simplified and the experimental purpose is satisfied.
5. Conclusions
Previous studies have modeled the yarn tension around the needle bending stage during the IP of woven carbon cloth, and the curve of the needle tip has been optimized for the different IP stages, but most of those studies did not consider comprehensively how the displacement is influenced by the fabric structure, the piercingneedle shape parameters, and the piercing template parameters. However, the combined variablesection cantilever beam deflection model based on a variable concentration load proposed herein does consider the above three factors comprehensively. In addition, previous studies have mainly optimized the tip shape curve. However, the present research results based on machine learning feature selection indicate that the order of influence of the needle shape parameters on the displacement should be d_{ned} > L_{2} > m > L_{1} > x. Therefore, the optimization of d_{ned} and L_{2} herein is very important and provides a theoretical basis for problems in engineering practice to determine the key influencing factors in specific feasible domains.
The yarn tension model herein compensates the yarn structure parameters for the yarn stretch and the Zdirection slip of the yarn during the piercing process, thereby improving the compatibility and applicability of the tension model. The mean and variance results from the needle displacement experiment based on machine vision show that the representative singleneedle displacement model for woven carbon cloth has obvious importance for reducing the displacement after being applied to the needle array. Therefore, the theoretical model is effective for reducing the collision between the steel needle and the hole wall, thereby indirectly improving the uniformity of the Zdirection fiber distribution.
Using an RPC fabric structure in the present study simplifies the IP model, thereby making it feasible to perform finiteelement simulations on the mechanics and motion mechanisms of the IP process. However, some limitations should be noted: the minimum structural elements corresponding to various structures of woven fabric are different, and the minimum structural elements of the needle array are determined by different engineering backgrounds. Future work should include (i) the piercing mechanism for a combination of multiple fabric structures and (ii) the structure with the smallest needle array unit.
Data Availability
The data used to support the findings of this study are included in the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors thank Hikvision for providing camera equipment and visual measurement technical support. This research was funded by the Tianjin Science and Technology Support Key Project (15ZCZDGX00840) (Research and development of the key technology and equipment for the whole piercing and forming automation of carbon fiber threedimensional fabric).
References
 L. Sun, C. Wang, X. F. Li, H. Li, M. Y. Ye, and C. An, “Research progress of C/C composite preforms,” Journal of Aeronautical Materials, vol. 38, no. 2, pp. 86–95, 2018. View at: Google Scholar
 P. G. Roincik, “146 Properties and application of mod3 pierced fabric reinforced carbon/carbon composites,” Carbon, vol. 11, no. 6, p. 690, 1973. View at: Google Scholar
 J. X. Zhu, “The sturetuarl charaeteirsties and properties of fine weave piecred fabire,” Aerospace Materials & Technology, vol. 1998, no. 1, pp. 41–43, 1998. View at: Google Scholar
 J. X. Zhu, Weaving Technology and Structural Parameter Optimization of Carbon Cloth with Whole Piercing Fabric, Southeast University, Nanjing, China, 2003.
 J. X. Zhu, J. M. He, and H. Y. Wang, “The mechanism of fiber bending and elongation in the integrated piercing process of orthogonal laminated woven fabrics,” Strategic Study of CAE, vol. 5, no. 5, pp. 59–69, 2003. View at: Google Scholar
 J. X. Zhu, J. M. He, Z. G. Zhou, and Y. Huang, “Optimizing of the form of steel needlepoint based on bending of woven in the integrated piercing process,” Journal of Tianjin University, vol. 37, no. 8, pp. 690–694, 2004. View at: Google Scholar
 J. X. Zhu, J. M. He, H. Y. Wang, and Z. G. Zhou, “Optimizing Zdirectional steely needle point based on fiber bending and elongation,” Strategic Study of CAE, vol. 5, no. 9, pp. 18–31, 2003. View at: Google Scholar
 J. Z. Dong, Z. Y. Tan, X. M. Jiang, and Z. J. Song, “Optimizing of puncture steely needle point Morphology duringthe stage of fiber bending and elongation needles,” Journal of Harbin Engineering University, vol. 40, no. 2, pp. 387–392, 2019. View at: Google Scholar
 Z. W. Qiao and C. Li, “Principle of optimization design for integral stiffness of needle matrix based on interated piercing,” in Proceedings of the 17th National Conference on Composite Materials, Beijing, China, 2012. View at: Google Scholar
 J. M. Hu and F. T. Hu, “Piercing damage of carbon fabrics containing metal powder versus needle tip shape,” Fiber Glass, vol. 3, no. 29, pp. 29–32, 2012. View at: Google Scholar
 Z. Y. Tan, Study on the Optimization of the Needle Shape and the Puncture Length of the Solid Fabric, Tiangong University, Tianjin, China, 2019.
 M. Salar, M. R. Ghasemi, and B. Dizangian, “A fast Gabased method for solving truss optimization problems,” International Journal of Optimization in Civil Engineering, vol. 6, pp. 101–114, 2015. View at: Google Scholar
 H. Ghasemi, H. Park, and T. Rabczuk, “A levelset based iga formulation for topology optimization of flexoelectric materials,” Computer Methods in Applied Mechnics and Engineering, vol. 313, pp. 239–258, 2016. View at: Publisher Site  Google Scholar
 S. K. Azad, “Seeding the initial population with feasible solutions in Metaheuristic optimization of steel trusses,” Engineering Optimization, vol. 50, no. 1, pp. 89–105, 2017. View at: Publisher Site  Google Scholar
 P. Zakian, “Metaheuristic design optimization of steel moment resisting frames subjected to natural frequency constraints,” Advances in Engineering Software, vol. 135, pp. 1–18, 2019. View at: Publisher Site  Google Scholar
 M. Bybordiani and S. K. Azad, “Optimum design of steel braced frames considering dynamic soilstructure interaction,” Structural and Multidisciplinary Optimization, vol. 60, pp. 1123–1137, 2019. View at: Publisher Site  Google Scholar
 P. Shoaei and M. Mahsuli, “Reliabilitybased design of steel moment frame structures isolated by leadrubber bearing systems,” Structures, vol. 20, pp. 765–778, 2019. View at: Publisher Site  Google Scholar
 A. Naderi, M. R. Sohrabi, M. R. Ghasemi, and B. Dizangian, “A swift technique for damage detection of determinate truss structures,” Engineering with Computers, vol. 2020, 2020. View at: Publisher Site  Google Scholar
 A. Kaveh, K. Biabani Hamedani, S. Milad Hosseini, and T. Bakhshpoori, “Optimal design of planar steel frame structures utilizing metaheuristic optimization algorithms,” Structures, vol. 25, pp. 335–346, 2020. View at: Publisher Site  Google Scholar
 P. Gou and J. Yu, “A nonlinear ann equalizer with minibatch gradient descent in 40gbaud pam8 im/dd system,” Optical Fiber Technology, vol. 46, pp. 113–117, 2018. View at: Publisher Site  Google Scholar
 H. Y. Pan, X. Niu, R. C. Li, Y. Dou, and H. Jiang, “Annealed gradient descent for deep learning,” Neurocomputing, vol. 380, pp. 201–211, 2019. View at: Publisher Site  Google Scholar
 Z. Yang, C. Wang, Z. Zhang, and J. Li, “Accelerated stochastic gradient descent with step size selection rules,” Signal Processing, vol. 159, pp. 171–186, 2019. View at: Publisher Site  Google Scholar
 A. Naderi, M. R. Sohrabi, M. R. Ghasemi, and B. Dizangian, “Total and partial updating technique: a swift approach for crosssection and geometry optimization of truss structures,” KSCE Journal of Civil Engineering, vol. 24, no. 4, pp. 1219–1227, 2020. View at: Publisher Site  Google Scholar
 J. Z. Dong, X. M. Jiang, J. C. Yang, S. H. Zhao, and R. W. Yuan, “Experimental study and development of integrated piercing steel needles array laying device of 3D fabric,” Journal of Textile Research, vol. 36, no. 3, pp. 115–120, 2015. View at: Google Scholar
 J. Z. Dong, B. L. Mei, X. M. Jiang, and J. Z. Yang, “Design and experimental study on steel needle gripper of replacementof Z directional steel needles in threedimensional fabric,” Journal of Textile Research, vol. 39, no. 12, pp. 101–106, 2018. View at: Google Scholar
 J. X. Zhu, “Analysis of the movement and mechanical behavior of carbon cloth in the whole puncture process,” Composite Materials: Innovation and Sustainable Development, vol. 2021, 2010. View at: Google Scholar
 J. B. Xie, X. M. Chen, Y. F. Zhang, G. D. Fang, and L. Chen, “Experimental and numerical investigation of the needling process for quartz fibers,” Composites Science and Technology, vol. 165, pp. 115–123, 2018. View at: Publisher Site  Google Scholar
 Y. Ma and Y. Zhang, “Deflection calculation of two section cantilever rectangular beam with variable cross section,” Structural Engineers, no. 1, pp. 1–5, 1994. View at: Google Scholar
 X. Z. Yuan and P. J. Liu, “Stiffness superposition of variable section cantilever beam,” Mechanics in Engineering, vol. 33, no. 2, pp. 8990, 2009. View at: Google Scholar
 L. P. Wang, “Method of deflection calculation of cantilever beam with variable cross section,” Henan Science and Technology, vol. 8, no. 22, pp. 122–124, 2018. View at: Google Scholar
 U. M. Khaire and R. Dhanalakshmi, “Stability of feature selection algorithm: a review,” Journal of King Saud UniversityComputer and Information Sciences, vol. 2019, pp. 1–14, 2019. View at: Google Scholar
 N. R. David, R. Yakir, K. F. Hilary, and R. G. Sharon, “Detecting novel associations in large data sets,” Science (New York, N.Y.), vol. 334, no. 6062, pp. 1518–1524, 2011. View at: Publisher Site  Google Scholar
 T. H. G. Megson, Structural and Stress Analysis, Elsevier, Amsterdam, Netherland, 2014.
 N. AlNuaimi, M. M. Masud, M. A. Serhani, and N. Zaki, “Streaming feature selection algorithms for big data: a survey,” Applied Computing and Informatics, pp. 1–13, 2019. View at: Google Scholar
 V. BolónCanedo and A. AlonsoBetanzos, “Ensembles for feature selection: a review and future trends,” Information Fusion, vol. 52, pp. 1–12, 2019. View at: Publisher Site  Google Scholar
 X. Qi, X. Fan, Y. Gao, and Y. Liu, “Learning bayesian network structures using weakest Mutualinformationfirst strategy,” International Journal of Approximate Reasoning, vol. 114, pp. 84–98, 2019. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Jingzhao Yang et al. 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.