Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2020 / Article

Research Article | Open Access

Volume 2020 |Article ID 5078590 |

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.

Reducing the Collision Damage Done to the Tips of Steel Needles during Integrated Piercing by Using Shape Optimization with Feature Selection

Academic Editor: Mohamed Shaat
Received30 Apr 2020
Revised12 Jul 2020
Accepted24 Jul 2020
Published24 Aug 2020


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 high-performance and heat-insulating C/C composite materials and are used widely in aircraft brake pads and rocket nozzles after curing C/C composites [15]. 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 [611]. 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 [1219]. 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 mini-batch gradient descent, annealing gradient descent, and accelerated stochastic gradient descent [2022]. 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 time-consuming process, and engineering applications need to combine experience and trial and error. Naderi et al. [23] proposed a multistage metaheuristic method called top-up 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 single-needle 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 Z-direction slippage, and (iii) the compensation of yarn tension caused by the Z-direction 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 variable-section 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 Unit-Cell 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 unit-cell (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.

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).

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 Z-direction displacements. Although in some ways similar to [5], the present study refines the model and considers the fiber slip in the vertical direction.

2.1.2. Bending Model of Fiber around Needle considering Z-Direction 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 i21, i22, and i23 have symmetry; therefore, only the length of the fiber B1B2 between i21 and i22 is analyzed in Figure 5, and the dotted line between C2B2′ is the projection of C2B2 on the plane.

The fiber length B1B2 after bending and stretching between two adjacent needles can be obtained from Figure 5 as follows:

The four-segment lengths that makeup B1B2 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, R2 is the same for tip 2, and dkj is the distance between the two holes adjacent to the center of the piercing template.where dk 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, plain-weave, twill-weave, or satin-weave 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 satin-structured 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, Lxj):where Z is the fiber path curve, C is the crimp height of the fabric structure phase, Lxj 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 Ec is the carbon cloth modulus of elasticity and Az is the carbon cloth fiber-bundle cross-sectional area. Satin-weave 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 dned 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, R1, R2, and dned 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 Fmax 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 Fmax, the elastic modulus Ec, and the cross-sectional area Az. It is known from engineering experience that yarn with a strain of less than 2 × 10−2 is stretched safely.


Ec2.3 × 1011PaModulus of elasticity
Dc9 × 10−4mmFiber diameter
η2 × 10−2Contact coefficient
k3 × 103Yarn number k
dned1.1mmNeedle bar diameter
dkj2.3mmDistance between holes
dk2.0mmHole diameter
Lxj5mmAdjacent yarn distance
C2.7 × 10−1mmFabric buckling wave height
τ1155Wave height adjustment factor
h0.5mmFiber Z-direction elongation

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 cross-sectional 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 cross-sectional 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.

Figure 7(a) shows that the needle tip curve passes through point B′ with coordinates (L1, dned), from which we can derivewhere L1 is the length of the needle tip, dned 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 F1 in the horizontal direction and F2 in the vertical direction. As shown in Figure 7(c), the F2 action point is translated to the intersection point C, which is the intersection point of the neutral layer of the steel needle and F2, and then, F1 can be translated to the tip vertex F1′. Next, a couple of M1 must be added, which is determined by the translation theorem of the force, and the couple moment of component F2 to point A is M2.

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 F2 is θ, then

From Figure 7(b), the tension components F1 and F2 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 AA″ is a dangerous section. To analyze the position of the dangerous point in the section AA″, the stress distributions of (a) the force component F2, (b) the additional force moment, and (c) the axial force component F1 are plotted in Figure 8.

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 cross-sectional 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 cross-sectional cantilever beam under concentrated load. The authors in [2830] consider how to calculate the deflection of various variable-section 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 cross-sectional needle bar, BO is a variable cross-sectional needle tip, and F2 acts at any point C. Figure 9(b) shows the additional force moment required for F2 to translate to point B, which is determined by the principle of force translation, namely,

Calculating the forced needle displacement is equivalent to calculating the deflection of the combined variable-section 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 cross-sectional cantilever beam has four parts, namely,where is the deflection under the concentrated load F2 after the tempered section AB (Figure 9(c)), is the deflection of the concentrated load F2 after the tempered segment BO is translated to point B and F2′ (Figure 9(d)), is deflection due to additional moment M of force F2 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 singular-point 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.

(1) Displacement v1and Rotation Angle θ1. The action of the fiber on the tip of the variable cross section is equivalent to F2 acting on any point C′ of the variable cross-sectional cantilever beam. In this case, the method of using the moment area is complicated, and the midsection moment of inertia of the straight-line cantilever beam after the second integral can be replaced by Ix 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 v2and 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 v3and 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 Mx, I, and AB to obtain

(4) Calculation of Cantilever Beam Deflection v. The displacement of L1 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 cross-validation.

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 up-to-date 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.

TypeScope of applicationComplexityRobustness

SpearmanLinear or simple monotonic nonlinearityLowMedium
KendallLinear or simple monotonic nonlinearityLowMedium
Maximum correlation coefficientLinear or nonlinearHighMedium
Distance dependentLinear or nonlinearMediumHigh
KNNLinear or nonlinearHighHigh
MICLinear or nonlinearLowLow

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 Fmax = 15 N and E = 2.11 × 1011 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 scikit-learn 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 dned > L2 > m > L1 > x, where m, dned, and x are negative correlations and L1 and L2 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, dned and L2.

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 L2 = 7 mm, dk = 2.3 mm, and β = 0.7, constraint conditions become as follows:

For Fmax = 0 to 50 N and dned = 0.8 to 2.0 mm, the contour relationship among , Fmax, and dned 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 Fmax ≤23 N, dned satisfies the displacement safety requirements of the current example in the range of 0.8 to 2.0 mm. For Fmax = 50 N, dned ≥1.035 mm meets the safety requirements of the current example of displacement.

(2) Needle Rod Length. Given the engineering conditions of dned = 1.1 mm, dk = 2.3 mm, and β = 0.7, equation (31) becomes as follows:

For Fmax = 0 to 50 N and L2 = 5 to 20 mm, the contour relationship among , Fmax, and L2 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 L2 = 20 mm, Fmax meets the displacement safety requirements of the current example within the range of 0 to 11.9 N. For Fmax = 50 N, L2 ≤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 L1 = 5 mm, L2 = 7 mm, σ = 1.514 × 108 Pa, and a safety factor of n = 2.5, equations (17)–(19) become as follows:

For Fmax = 0 to 50 N and dned = 0.8 to 2.0 mm, the contour relationships among σ, Fmax, and dned are as shown in Figure 13.

Figures 13(a) and 13(b) show that the areas (i) above and to the right of the 6.056 × 107 strength contour at point A′ for 0 to 10 N and (ii) above and to the left of the −6.056 × 107 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 × 107 strength contour at point A″ for 0–10 N and (ii) above and to the left of the 6.056 × 107 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 × 107 Pa and −6.056 × 107 Pa are given in Table 6.

PointConstraint planeFmaxdned

A′Z = 6.056 × 107 Pa≥3.18 N0.8 × 10−3 m
Z = 6.056 × 107 Pa1 N≥1.1 × 10−3 m
Z = −6.056 × 107 Pa≤27.62 N0.8 × 10−3 m
Z = −6.056 × 107 Pa50 N≥1.159 × 10−3 m
A″Z = 6.056 × 107 Pa≤44 N0.8 × 10−3 m
Z = 6.056 × 107 Pa50 N≥0.86 × 10−3 m
Z = −6.056 × 107 Pa≥3.33 N0.8 × 10−3 m
Z = −6.056 × 107 Pa1 N≥1.1 × 10−3 m

(2) Needle Rod Length. Given that the commonly used conditions of the project are L1 = 5 mm, dned = 1.1 mm, σ = 1.514 × 108 Pa, and a safety factor of n = 2.5, equations (17) and (19) become

For Fmax = 0 to 50 N and L2 = 5 to 20 mm, the contour relationships among σ, Fmax, and L2 are as shown in Figure 14.

Figure 14(a) shows that the area below and to the right of the 6.056 × 107 strength contour at point A′ for 0 to 10 N is the safe range allowed by the needle strength. However, the −6.056 × 107 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 × 107 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 × 107 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 × 107 and −6.056 × 107 are given in Table 7.

PointConstraint planeFmaxL2

A′Z= 6.056 × 107 Pa≥2.73 N20 mm
Z= 6.056 × 107 Pa1 N≥7.07 mm
Z= −6.056 × 107 Pa≤46.36 N20 mm
Z= −6.056 × 107 Pa50 N≤27.5 mm
A″Z= 6.056 × 107 Pa≥79.98 N (not in defined domain)20 mm
Z= 6.056 × 107 Pa50 N−121.94 mm (not in defined domain)
Z= −6.056 × 107 Pa≥2.79 N20 mm
Z= −6.056 × 107 Pa1 N≥7.01 mm

3.2.3. Optimization Results of dned and L2 Example

The paper gives the range of parameters allowed by the process: Fmax = 0 to 50 N, dned = 0.8 to 20 mm, L2 = 5 to 20 mm, and the feasible domain that satisfies the displacement and strength constraint mentioned above is dned = 1.159 to 2.0 mm and L2 = 7.07 to 7.96 mm.

The calculation parameters in Table 8 are substituted into equation (31) to obtain the partial derivatives of dned and L2 as follows:


EElastic modulus of steel needle2.11 × 1011 Pa
mExponent in power function of tip shape curve2
L1Tip length5 × 10−3 m
xDistance between fiber and tipL1

Equation (37) shows that displacement is decremented monotonically with dned, and increases monotonically with L2. 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 dned = 1.159 × 10−3 m and L2 = 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.

Group no.L1 (mm)mFmax (N)E (Pa)L2 (mm)dned (mm) (mm)

152152.11 × 1097.951.1590.206
252152.11 × 109101.1590.289
352152.11 × 1097.951.00.321
452152.11 × 109101.00.450

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.


Piercing templateOrifice size10 cm × 10 cm
Hole diameter2.3 mm
Distance between holes2.2 mm
Orifice thickness10 mm
Steel needleRod diameter1.0 mm/1.159 mm
Rod length7.95 mm/10 mm
Exponent in power function of tip shape curve2
Tip length5 mm
Woven carbon clothStructureWoven satin warp 8/3
Laying modeOrthogonal laminated


Measurement requirementsFOV10 cm × 10 cm
Accuracy0.1 mm
Camera parametersCamera typeMV-CE200-10 GM
Sensor modelIMX183
Resolution5472 × 3036
CCD chip size1″
Phase size2.4 μm
Lens parametersLens typeMVL-KF2528M-12MP
Optical distortion0.40%
Focal length25 mm
Light sourceLight typeMV-LGSS-175-W
Light source controllerMV-LEVS-24-1-SY

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 image-acquisition 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, human-computer 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.

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.

Group no.Picture no.L2 (mm)dned (mm)Safety holesDangerous holesTotal number of holesRecognition rate (%)Safety rate (%)Dangerous (%)


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.

Group no.Picture no.L2 (mm)dned (mm)Sum of safety distanceSum of hazard distanceMeanVariancePercentage reduction of meanPercentage reduction of variance


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.

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 piercing-needle shape parameters, and the piercing template parameters. However, the combined variable-section 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 dned > L2 > m > L1 > x. Therefore, the optimization of dned and L2 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 Z-direction 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 single-needle 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 Z-direction fiber distribution.

Using an RPC fabric structure in the present study simplifies the IP model, thereby making it feasible to perform finite-element 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.


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 three-dimensional fabric).


  1. 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
  2. P. G. Roincik, “146 Properties and application of mod-3 pierced fabric reinforced carbon/carbon composites,” Carbon, vol. 11, no. 6, p. 690, 1973. View at: Google Scholar
  3. 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
  4. J. X. Zhu, Weaving Technology and Structural Parameter Optimization of Carbon Cloth with Whole Piercing Fabric, Southeast University, Nanjing, China, 2003.
  5. 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
  6. 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
  7. J. X. Zhu, J. M. He, H. Y. Wang, and Z. G. Zhou, “Optimizing Z-directional 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
  8. 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
  9. 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
  10. 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
  11. Z. Y. Tan, Study on the Optimization of the Needle Shape and the Puncture Length of the Solid Fabric, Tiangong University, Tianjin, China, 2019.
  12. M. Salar, M. R. Ghasemi, and B. Dizangian, “A fast Ga-based method for solving truss optimization problems,” International Journal of Optimization in Civil Engineering, vol. 6, pp. 101–114, 2015. View at: Google Scholar
  13. H. Ghasemi, H. Park, and T. Rabczuk, “A level-set 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
  14. 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
  15. P. Zakian, “Meta-heuristic 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
  16. M. Bybordiani and S. K. Azad, “Optimum design of steel braced frames considering dynamic soil-structure interaction,” Structural and Multidisciplinary Optimization, vol. 60, pp. 1123–1137, 2019. View at: Publisher Site | Google Scholar
  17. P. Shoaei and M. Mahsuli, “Reliability-based design of steel moment frame structures isolated by lead-rubber bearing systems,” Structures, vol. 20, pp. 765–778, 2019. View at: Publisher Site | Google Scholar
  18. 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
  19. A. Kaveh, K. Biabani Hamedani, S. Milad Hosseini, and T. Bakhshpoori, “Optimal design of planar steel frame structures utilizing meta-heuristic optimization algorithms,” Structures, vol. 25, pp. 335–346, 2020. View at: Publisher Site | Google Scholar
  20. P. Gou and J. Yu, “A nonlinear ann equalizer with mini-batch gradient descent in 40gbaud pam-8 im/dd system,” Optical Fiber Technology, vol. 46, pp. 113–117, 2018. View at: Publisher Site | Google Scholar
  21. 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
  22. 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
  23. A. Naderi, M. R. Sohrabi, M. R. Ghasemi, and B. Dizangian, “Total and partial updating technique: a swift approach for cross-section 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
  24. 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 3-D fabric,” Journal of Textile Research, vol. 36, no. 3, pp. 115–120, 2015. View at: Google Scholar
  25. 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 three-dimensional fabric,” Journal of Textile Research, vol. 39, no. 12, pp. 101–106, 2018. View at: Google Scholar
  26. 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
  27. 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
  28. 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
  29. X. Z. Yuan and P. J. Liu, “Stiffness superposition of variable section cantilever beam,” Mechanics in Engineering, vol. 33, no. 2, pp. 89-90, 2009. View at: Google Scholar
  30. 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
  31. U. M. Khaire and R. Dhanalakshmi, “Stability of feature selection algorithm: a review,” Journal of King Saud University-Computer and Information Sciences, vol. 2019, pp. 1–14, 2019. View at: Google Scholar
  32. 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
  33. T. H. G. Megson, Structural and Stress Analysis, Elsevier, Amsterdam, Netherland, 2014.
  34. 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
  35. V. Bolón-Canedo and A. Alonso-Betanzos, “Ensembles for feature selection: a review and future trends,” Information Fusion, vol. 52, pp. 1–12, 2019. View at: Publisher Site | Google Scholar
  36. X. Qi, X. Fan, Y. Gao, and Y. Liu, “Learning bayesian network structures using weakest Mutual-information-first strategy,” International Journal of Approximate Reasoning, vol. 114, pp. 84–98, 2019. View at: Publisher Site | Google Scholar

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.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.