Modeling of Precise Tension with Passive Dancers for Automated Fiber Placement

: This paper modeled the tension fluctuation during automated fiber placement (AFP), which depicted the tension variations under different operating conditions. The stability and validity of the model were demonstrated using Bode plots and experiments, respectively. According to the model, the tension fluctuations of AFP at different stages were obtained. Additionally, the passive dancer parameters with the better system performance were determined using the evaluation methodology presented in this paper. Moreover, it was discovered that the damping coefficient affects the tension variation more significantly than the elasticity coefficient. Finally, the placement experiments showed that the determined passive dancer parameters improved the laying quality significantly.


Introduction
Resin-based carbon-fiber-reinforced composites have been widely used in aerospace applications due to their lightweight, high strength, and superior specific strength and modulus [1][2][3][4].Automated fiber placement (AFP) is one of the most significant processes for manufacturing composite parts, and its manufacturing process is shown in Figure 1a.During AFP, each tow has separate transport and cutting channels.After the tows have been transferred onto the surface of the mold, they are placed side by side on the mold under the combined action of the compacting roller and the heating device [5,6].
Tension is one of the most critical parameters affecting placement quality.Applying appropriate tension to the tows during the laying process helps to maintain the smoothness of the tows and reduces the probability of folding and bridging defects to some extent.Moreover, it can moderately enhance the mechanical properties of the components after curing [7][8][9][10][11].Thus, the tension should be stabilized as much as possible within a suitable range during AFP.Firstly, frequent tension fluctuations could lead to the bridging, folding, and excessive slipping distances of the tow during the laying process [12,13].Secondly, repeated tension fluctuations increase the probability of wire blocking, which can lead to placement failure.A stable tension can reduce laying defects and the likelihood of laying failures, thus enhancing the quality of the laying [7,[14][15][16][17].
Many factors affect the laying quality, such as laying path, laying pressure, laying tension, and laying speed [18][19][20].However, there have been few studies on tension fluctuations during layup to.Most studies on layup treat tension as a stable value and ignore its fluctuations, but this is incorrect.During the clamping, retransfering, and layup stages, the tow tension fluctuates dramatically due to changes in speed, mold curvature, Many factors affect the laying quality, such as laying path, laying pressure, laying tension, and laying speed [18][19][20].However, there have been few studies on tension fluctuations during layup to.Most studies on layup treat tension as a stable value and ignore its fluctuations, but this is incorrect.During the clamping, retransfering, and layup stages, the tow tension fluctuates dramatically due to changes in speed, mold curvature, etc.And according to the previous paragraph, such fluctuations could affect the placement quality [21].However, the pattern of tension fluctuations is not yet clear.Therefore, it is necessary to establish a precise tension fluctuation model for the AFP system, which has not been extensively investigated yet.Dodangeh et al. [19] developed two setups for layup and measured tow tackiness under different laying situations.At the same time, they optimized the parameters of the layup process using the Taguchi method and experimentally validated these optimal values.The optimal combination of the process parameters predicted by the Taguchi method was a heat gun temperature of 400 °C, a compaction force of 178 N, and a feed speed of 50 mm/sec using a roller of 60 durometers.The study results indicated that the interplay between the process parameters affected the tow viscosity more than the process parameters acting alone.By analyzing the contact between the tow and compaction roller at any laying path point, Wang et al. [22] proposed a mechanical model of compression force distribution between the compaction roller and the tow for laying arbitrary curvature molds.Based on this model, a numerical calculation algorithm was presented to calculate the pressure distribution during AFP.Meanwhile, they verified the model's validity by measuring the pressure distribution.Budelmann et al. [23] used a single-stage peel method to identify the correlation between key processing parameters and tow viscosity at varying degrees of room temperature over time.Qu et al. [24] performed the required analysis for laying suitability before path planning.The results showed that choosing the appropriate path parameters can decrease the risk of unreasonable path planning and enhance the effectiveness.Through tow feeding and laying tests, Xiao et al. [21] investigated the tension during the laying process.The results of tows feeding trials showed that the tension greatly affected the tows' feeding reliability and precision.In addition, the laying tests indicated that an unsuitable tension reduced the layup quality.However, this paper still did not consider tension fluctuations.Jiang et al. [18] employed segmented compaction rollers in their study.By analyzing the contact relationship between the segmented compaction rollers and the tow, they established a theoretical model for arbitrary layer pressure distribution and conducted experimental validation.Based on this model, they proposed the concepts of single-point and whole-path pressure uniformity to evaluate the pressure distribution quantitatively.At present, the majority of studies assume that the tension remains constant during layup.Few papers have focused on the tension fluctua- Dodangeh et al. [19] developed two setups for layup and measured tow tackiness under different laying situations.At the same time, they optimized the parameters of the layup process using the Taguchi method and experimentally validated these optimal values.The optimal combination of the process parameters predicted by the Taguchi method was a heat gun temperature of 400 • C, a compaction force of 178 N, and a feed speed of 50 mm/sec using a roller of 60 durometers.The study results indicated that the interplay between the process parameters affected the tow viscosity more than the process parameters acting alone.By analyzing the contact between the tow and compaction roller at any laying path point, Wang et al. [22] proposed a mechanical model of compression force distribution between the compaction roller and the tow for laying arbitrary curvature molds.Based on this model, a numerical calculation algorithm was presented to calculate the pressure distribution during AFP.Meanwhile, they verified the model's validity by measuring the pressure distribution.Budelmann et al. [23] used a single-stage peel method to identify the correlation between key processing parameters and tow viscosity at varying degrees of room temperature over time.Qu et al. [24] performed the required analysis for laying suitability before path planning.The results showed that choosing the appropriate path parameters can decrease the risk of unreasonable path planning and enhance the effectiveness.Through tow feeding and laying tests, Xiao et al. [21] investigated the tension during the laying process.The results of tows feeding trials showed that the tension greatly affected the tows' feeding reliability and precision.In addition, the laying tests indicated that an unsuitable tension reduced the layup quality.However, this paper still did not consider tension fluctuations.Jiang et al. [18] employed segmented compaction rollers in their study.By analyzing the contact relationship between the segmented compaction rollers and the tow, they established a theoretical model for arbitrary layer pressure distribution and conducted experimental validation.Based on this model, they proposed the concepts of single-point and whole-path pressure uniformity to evaluate the pressure distribution quantitatively.At present, the majority of studies assume that the tension remains constant during layup.Few papers have focused on the tension fluctuations in AFP.
Numerical simulation is one of the most critical methods for scientific research [25,26].In this paper, an accurate tension fluctuation model is developed to depict the tension variations during the layup process.The stability and validity of the model are demonstrated by the Bode plot and experiment (Sections 3.1 and 3.2, respectively).After that, a MATLAB/Simulink simulation analysis was carried out to illustrate the superiority of the tension control system with passive rollers (Section 3.2).At the same time, the tension fluctuations under different elasticity coefficients and damping are obtained for the clamping, retransmission, and layup stages (Section 3.3).The evaluation method established in this paper determines the passive dancer parameters with a better system performance, which is apply to the tow material and laying mold selected in this article (Section 3.3).In addition, the effect of the damping coefficient on the tension fluctuation is more significant than that of the elasticity coefficient by simulation analysis (Section 3.3).Finally, the placement experiments were carried out in this paper using the passive dancer parameters determined by the simulation (Section 3.4).The experiment results show that the selected parameters improve the laying quality significantly (Section 3.4).
The tension fluctuation model developed in this paper accurately reflected the tension variation at each stage in AFP.Through this model, suitable passive dancer parameters were selected for different layup molds and tow materials, minimizing tension fluctuations and improving the layup quality.

Model of the Tension Fluctuation for AFP
The range of tension fluctuations significantly impacts the stability and reliability of laying.Therefore, it is crucial to model the tension fluctuations during AFP.Considering the working conditions, several issues must be explored when modeling tension fluctuations.Firstly, the clamping module, the retransmission module, and the cut module will frequently work alternately during the laying process, which will cause frequent abrupt changes in tension.Secondly, the laying speed and tension frequently change when the large curvature mold is laid.Thus, to address the above problems, this paper established a tension fluctuation model in the laying process, which reflected the tension variation rule under different working conditions.The schematic diagram of the tension control system is shown in Figure 1b.

The Model for the Pay-off Axis
Figure 2a shows the schematic diagram of the pay-off axis working.The motor used for the pay-off axis was a permanent magnet synchronous servo motor (PMSM), which directly drove the pay-off axis.Therefore, the pay-off axis and servo motor were modeled as a whole in this paper.The space magnetic field distribution was sine, and the magnetic circuit's core magnetic saturation and nonlinear characteristics were ignored.Thus, according to the torque balance equation, the equation of motion of the pay-off axis was obtained as follows: where M t is the resistance torque of motor output, T(t) is the tow tension, r(t) is the radius of the pay-off axis, M f is the dry friction torque, and β is the viscous friction coefficient.J(t) is the rotational inertia and ω(t) is the rotating angular velocity of the pay-off axis.
According to Equation (1), the equation of motion of the pay-off axis is strongly timevarying.Hence, to simplify the problem, the following simplification was considered for Equation (1): (1) Since within the experimental scope of this study, the variation in the radius of the pay-off axis was not significant and the change in radius had a negligible impact on the model established in this study, the radius of the pay-off axis was considered a constant in this study; (2) Based on the simplification assumption (1), the change in the rotational inertia could also be neglected; (3) Due to the characteristics of the PMSM, the dry friction torque of the system was not considered.According to Equation (1), the equation of motion of the pay-off axis is strongly timevarying.Hence, to simplify the problem, the following simplification was considered for Equation ( 1): (1) Since within the experimental scope of this study, the variation in the radius of the pay-off axis was not significant and the change in radius had a negligible impact on the model established in this study, the radius of the pay-off axis was considered a constant in this study; (2) Based on the simplification assumption (1), the change in the rotational inertia could also be neglected; (3) Due to the characteristics of the PMSM, the dry friction torque of the system was not considered.
Therefore, Equation (1) was simplified as follows: The PMSM operated in the torque mode.Hence, the current vector control was adopted, and the relationship between current and voltage was expressed as [27]: where  and  are the current and voltage, respectively. is the inductance; , , and  are the motor winding resistance, the number of motor poles, and rotor flux, respectively.
Meanwhile, based on the electromagnetic torque equation, the relationship between the current and output torque was expressed as follows [27]: The back-EMF coefficient and torque constant of the motor were defined as [27]: where  and  are the back-EMF coefficient and torque constant in turn, respectively.Therefore, Equation (1) was simplified as follows: The PMSM operated in the torque mode.Hence, the current vector control was adopted, and the relationship between current and voltage was expressed as [27]: where i and u are the current and voltage, respectively.L is the inductance; R, p, and ψ f are the motor winding resistance, the number of motor poles, and rotor flux, respectively.Meanwhile, based on the electromagnetic torque equation, the relationship between the current and output torque was expressed as follows [27]: The back-EMF coefficient and torque constant of the motor were defined as [27]: where K e and K t are the back-EMF coefficient and torque constant in turn, respectively.By joining Equation (2) to (5), the differential equation between the rotating angular velocity and the voltage was acquired:
We considered the small inductance L of the PMSM and Rβ ≪ K t K e .Thus, L and Rβ were ignored in this paper.Equation ( 6) was simplified as follows: where τ m = RJ K t K e , K c = R K t , and K m = 1 K e .The Laplace transform was applied to both ends of Equation (7).Hence, the transfer function of the motor was obtained as follows: If the output of the system is velocity, Equation ( 8) can be written as:

The Model for Velocity Difference-Tension
From the previous section, the speed of the pay-off axis will directly affect the tow tension.Meanwhile, for the AFP system, there is a strong coupling between the tow tension and the difference in speed.Therefore, the velocity difference-tension model is established in this section to achieve tension control.
To simplify the calculation, the AFP process shown in Figure 1a was simplified to the model shown in Figure 2b.In Figure 2b, index "1" means the tow is on the pay-off axis, "2" indicates that the tow is in the area where tension control is required, and "1-2" is the length of the laying area (the same below).In addition, the following laws need to be considered when establishing the velocity difference-tension model: (1) Hooke's law [28], which is used to define the elastic properties of the tow; (2) The law of conservation of mass [29], which represents the cross-coupling between the velocity and tension of the tow; (3) Newton's laws of motion [30] described the dynamics of the passive dancers.We considered that there was no mass transfer between the substances.Thus, by the law of conservation of mass, the following equation was obtained: where ρ, A, and v are the density, cross-sectional area, and velocity of the tow, respectively.dx is the length of the tow in the unit area.
Since the tow is a transversely isotropic viscoelastic material [31], in the direction of motion, the tow is linearly elastic.The deformation of the tow can be obtained as follows: where ε is the tow strain.u and n are the tensioned and non-tensioned states of the tow, respectively.
Assuming the mass of the tow remains constant, then Combining Equations ( 10)-( 12) produces the following: Considering the strain of the tow is equal everywhere in the control area, Equation ( 13) can be written as: During the tensile state, the tow strain on the pay-off axis is considered constant and much less than 1.From this, Equation ( 14) can be simplified to: where L 0 is the total length of the tow in the control area.The simplification of Equation ( 15) produces: According to the ε ≪ 1, we obtained: The loss between the tow and each mechanism was ignored.Substituting Equation ( 17) into ( 16) and combined with Hooke's law (T(t) = EAε 2 (t)), the following equation was ob- tained: where E is the tow elastic modulus.

The Model of Passive Dancers
The schematic diagram of the passive dancers is shown in Figure 2c.We obtained from the geometric relationship that where x(t) and δL 1 are the passive dancer movement distance and upstream variation in the tow length, respectively.Since δθ 1 (t) is small during AFP, Equation ( 19) can be simplified as: At the same time, The passive dancer is subjected to tow tension, spring, and damping forces.Hence, the following equation is acquired by Newton's second law of motion: where m, D, and K are the mass, damping, and elasticity coefficient of the passive dancer, respectively.
According to the geometric relationship shown in Figure 2c, it could also be obtained that: From the above, it is evident that x(t) is small.Thus, Equation ( 23) can be simplified and written as: Similarly, it can be obtained that: where , so that Equation ( 22) can be expressed as: The Laplace transform was applied to both ends of Equation (26).Therefore, the transfer function of passive dancers could be written as follows:

The Model of the Tension Control with Passive Dancers
During AFP, passive dancers are in the left and right movements.Combining Equations ( 20) and ( 21), the total length of the tow can be defined as follows: Let 1 sin θ 1 + 1 sin θ 2 = η.Taking Equation ( 28) into ( 15), the following equation can be obtained: In addition, Equation ( 27) can be written as: Similar to Equation ( 18), using Hooke's law, it can be obtained as follows: Linearizing Equation (31), let V n = v n − v nr , T l = T − t r , and X = x − x r .Substituting V n , T l , and X into Equation (31), the linearized tensorial dynamics can be derived as: Defining L 0 − ηx r = γ and applying the Laplace transform to Equation (32) produces: Substituting Equation (34) into Equation (33), the following transfer function is derived The tension sensor can be treated as an amplifier link due to the short delay.Therefore, the transfer function of the sensor is When tension is controlled during AFP, nonlinear and time-variant problems are insignificant.Hence, the integral separated PID was chosen as the system's control strategy.The control law is shown below: where u(t) is the controlled voltage.T t , λ, and e 0 is the sampling time, control coefficient of integral terms, and set threshold, respectively.Combining Equations ( 9), ( 27), ( 35), ( 36) and ( 38), the control framework of the system can be obtained, as shown in Figure 3b.If the passive dancer is disregarded in this system, the control diagram is as shown in Figure 3a.
The tension sensor can be treated as an amplifier link due to the short delay.Therefore, the transfer function of the sensor is When tension is controlled during AFP, nonlinear and time-variant problems are insignificant.Hence, the integral separated PID was chosen as the system's control strategy.The control law is shown below: where () is the controlled voltage.  , , and  0 is the sampling time, control coefficient of integral terms, and set threshold, respectively.Combining Equations ( 9), ( 27), ( 35), (36), and (38), the control framework of the system can be obtained, as shown in Figure 3b.If the passive dancer is disregarded in this system, the control diagram is as shown in Figure 3a.
The transfer function derivation method of the system that neglects the passive rollers is similar to the above.Thus, this paper's derivation of this transfer function was not carried out.The system equations are as follows:  The transfer function derivation method of the system that neglects the passive rollers is similar to the above.Thus, this paper's derivation of this transfer function was not carried out.The system equations are as follows:

Analysis of the Precise Tension Control Model
In this section, Simulink simulates the spring damping and elasticity coefficient effects on the tension fluctuations at different stages of the AFP.The passive dancer parameters that make the system perform better were selected by analyzing the results.The specific simulated stages were clamping, retransmission, and laying.This paper simplified these processes as tension step signals, velocity perturbations, and the combination of tension step and velocity step, respectively.The required parameters for the simulation are shown in Table 1, and the step function is illustrated in Figure 4.In this section, Simulink simulates the spring damping and elasticity coefficient effects on the tension fluctuations at different stages of the AFP.The passive dancer parameters that make the system perform better were selected by analyzing the results.The specific simulated stages were clamping, retransmission, and laying.This paper simplified these processes as tension step signals, velocity perturbations, and the combination of tension step and velocity step, respectively.The required parameters for the simulation are shown in Table 1, and the step function is illustrated in Figure 4.  Figure 5 shows the experimental platform of the automatic wire-feeding used to carry

Model Validation
Figure 5 shows the experimental platform of the automatic wire-feeding used to carry out the experiments for verifying the model's validity.The tow used for experiments was a unidirectional carbon/epoxy tow with a width of 6.35 mm (Type: USN12500/7901).The parameters of the tow are shown in Table 1.In Figure 5, the experimental platform is equipped with two rotating axes, the pay-off axis and the backing paper axis, driven by torque servo motors.The pay-off axis holds the filament roll, while the backing paper axis holds the backing paper roll.The filament roll is the source of tows.As the filament pay-off axis rotates, the tows within the filament roll unwind, allowing the tows to enter the filament bundle conveying mechanism.Meanwhile, the backing paper that wraps around the tows separates from the tows and is collected back onto the backing paper roll.After leaving the filament roll, the tows pass through the guide roller, several guide wheels, and tension sensing device and smoothly enter the filament bundle action mechanism.The bundle action mechanism consists of a clamp module, a tow feeding system, and a cut module.The clamping module is responsible for fixing the position of the tow when the feeding system is not operational.The feeding system is accountable for transporting the tow forward, while the cut module is responsible for cutting the tow.
guide rail until a state of force equilibrium is achieved.The tension-sensing device can convert tension into pressure on the guide wheel, indirectly measuring the tension of the bundles.The clamping module is driven by a cylinder, where the extension and retraction of the cylinder's piston rod drive the closure of the clamping block, thereby achieving the function of fixing the tow.The tow feeding system consists of a pair of mutually movable rollers, with the lower roller driven by a motor, rotating in real-time.The tows are positioned between the two rollers.When the upper roller clamps the tows between the rollers under the action of the cylinder, the tows are driven forward by the frictional force of the rollers, following the linear velocity of the rollers.The cut module is propelled by a cylinder to cut the tows using the cutting blade.
The IPC, motor, and sensor models used in this experiment platform were C6130, AKM33E, and FS127, respectively.Furthermore, the system was designed as a function block and utilized the EtherCAT Fieldbus protocol for communication.The tension control flow chart is shown in Figure 6.It can be obtained from the flow chart that the system essentially controls the tension by controlling the motor torque.Thus, the motor torque needs to be calibrated before the experiment, and the calibration result is: where  is the digital output of the PLC and  is the actual output torque of the motor.Equation (40) calculates the torque in the flow chart.The operating principle of the experimental platform shown in Figure 5 is as follows: The unwound tows from the filament roll smoothly enter the spring roller mechanism through the reversing device.The spring roller consists of a guide wheel, a translating guide rail, and a spring.The tension of the tows drives the guide wheel to move along the guide rail until a state of force equilibrium is achieved.The tension-sensing device can convert tension into pressure on the guide wheel, indirectly measuring the tension of the bundles.The clamping module is driven by a cylinder, where the extension and retraction of the cylinder's piston rod drive the closure of the clamping block, thereby achieving the function of fixing the tow.The tow feeding system consists of a pair of mutually movable rollers, with the lower roller driven by a motor, rotating in real-time.The tows are positioned between the two rollers.When the upper roller clamps the tows between the rollers under the action of the cylinder, the tows are driven forward by the frictional force of the rollers, following the linear velocity of the rollers.The cut module is propelled by a cylinder to cut the tows using the cutting blade.
The IPC, motor, and sensor models used in this experiment platform were C6130, AKM33E, and FS127, respectively.Furthermore, the system was designed as a function block and utilized the EtherCAT Fieldbus protocol for communication.
The tension control flow chart is shown in Figure 6.It can be obtained from the flow chart that the system essentially controls the tension by controlling the motor torque.Thus, the motor torque needs to be calibrated before the experiment, and the calibration result is: where x is the digital output of the PLC and y torque is the actual output torque of the motor.Equation (40) calculates the torque in the flow chart.The pressure, target tension, and retransmission speed were set as 0.6 MPa, 5 N, and 60 mm/s, respectively.The relative humidity was 65%, and the ambient temperature was 25 °C.Based on the automatic wire-feeding experimental platform established in Figure 5, the following three experiments were carried out to verify the validity of the tension system model established in this paper: (1) The tow was clamped in the stationary state, after which a tension step of 5 N was applied.This process is depicted in Figure 7a.It can also be observed from Figure 7a that the actual tension values overlap very well with the simulated values.
(2) The lower roller of the feeding module was set to 60 mm/s.At the moment of feeding module initiation, the upper roller was pressed down, and the clamping block was lifted, causing the tow to be driven forward by the friction of the roller.Due to the tow's significant friction and low mass, it rapidly accelerated from zero to 60 mm/s with a very high acceleration.Therefore, the moment of feeding module initiation can be considered as a velocity step process.After a brief yet significant disturbance, the tension quickly stabilized and fluctuated within a small range near 5 N.After maintaining a steady delivery speed of 60 mm/s for a while, the lower roller speed was reduced to 0, resulting in another significant fluctuation in tension.However, the tension quickly returned to 5 N. Figure 7b records the tension fluctuations throughout this process.The experimental curve shows a high degree of agreement with the simulated curve.In particular, as revealed in Figure 7b, small tension fluctuations occur during wire feeding, which is unavoidable and acceptable.
(3) The tow feeding speed was set as a sine wave.For the tension system, the tension value is outputted as a sinusoidal wave with the same frequency, as shown in Figure 7c.It can also be observed from Figure 7c that the tension amplitude and the phase difference with the speed signal exhibit a high degree of agreement between the experimental and simulated results.
In conclusion, the tension system model established in this study has a high credibility.The pressure, target tension, and retransmission speed were set as 0.6 MPa, 5 N, and 60 mm/s, respectively.The relative humidity was 65%, and the ambient temperature was 25 • C. Based on the automatic wire-feeding experimental platform established in Figure 5, the following three experiments were carried out to verify the validity of the tension system model established in this paper: (1) The tow was clamped in the stationary state, after which a tension step of 5 N was applied.This process is depicted in Figure 7a.It can also be observed from Figure 7a that the actual tension values overlap very well with the simulated values.

Analysis of the System Stability and Superiority
Illustration: In classical control theory, the dynamic performance of a model is typically described using step response, and the most important parameters in the step response are the maximum overshoot (  ), rise time (  ), and settling time (  ).Therefore, in this study, we considered using these three parameters to assess the quality of the established tension system.Specific analyses are provided below.
Before performing the simulation, it is necessary to determine the stability and superiority of the model with the passive dancer developed in this paper.Figure 8 shows the Bode diagram of the system (without a controller) with different damping and elasticity coefficients.In Figure 8a, both   (gain margin) and   (phase margin) are infinite, which indicates that the established model is steady.Moreover, the larger the  , the (2) The lower roller of the feeding module was set to 60 mm/s.At the moment of feeding module initiation, the upper roller was pressed down, and the clamping block was lifted, causing the tow to be driven forward by the friction of the roller.Due to the tow's significant friction and low mass, it rapidly accelerated from zero to 60 mm/s with a very high acceleration.Therefore, the moment of feeding module initiation can be considered as a velocity step process.After a brief yet significant disturbance, the tension quickly stabilized and fluctuated within a small range near 5 N.After maintaining a steady delivery speed of 60 mm/s for a while, the lower roller speed was reduced to 0, resulting in another significant fluctuation in tension.However, the tension quickly returned to 5 N. Figure 7b records the tension fluctuations throughout this process.The experimental curve shows a high degree of agreement with the simulated curve.In particular, as revealed in Figure 7b, small tension fluctuations occur during wire feeding, which is unavoidable and acceptable.
(3) The tow feeding speed was set as a sine wave.For the tension system, the tension value is outputted as a sinusoidal wave with the same frequency, as shown in Figure 7c.It can also be observed from Figure 7c that the tension amplitude and the phase difference with the speed signal exhibit a high degree of agreement between the experimental and simulated results.
In conclusion, the tension system model established in this study has a high credibility.

Analysis of the System Stability and Superiority
In classical control theory, the dynamic performance of a model is typically described using step response, and the most important parameters in the step response are the maximum overshoot (M p ), rise time (T r ), and settling time (T s ).Therefore, in this study, we considered using these three parameters to assess the quality of the established tension system.Specific analyses are provided below.
Before performing the simulation, it is necessary to determine the stability and superiority of the model with the passive dancer developed in this paper.Figure 8 shows the Bode diagram of the system (without a controller) with different damping and elasticity coefficients.In Figure 8a, both G m (gain margin) and P m (phase margin) are infinite, which indicates that the established model is steady.Moreover, the larger the D, the wider the medium frequency, and the smoother the system in the medium frequency.Meanwhile, in the high band, the system decays faster with a larger D, which demonstrates the ability of the system to resist high-frequency interference.The pattern in Figure 8b is similar to that in Figure 8a.Thus, Figure 8b was not analyzed in this paper.

Analysis of the System Stability and Superiority
Illustration: In classical control theory, the dynamic performance of a model is typically described using step response, and the most important parameters in the step response are the maximum overshoot (  ), rise time (  ), and settling time (  ).Therefore, in this study, we considered using these three parameters to assess the quality of the established tension system.Specific analyses are provided below.
Before performing the simulation, it is necessary to determine the stability and superiority of the model with the passive dancer developed in this paper.Figure 8 shows the Bode diagram of the system (without a controller) with different damping and elasticity coefficients.In Figure 8a, both   (gain margin) and   (phase margin) are infinite, which indicates that the established model is steady.Moreover, the larger the  , the wider the medium frequency, and the smoother the system in the medium frequency.Meanwhile, in the high band, the system decays faster with a larger , which demonstrates the ability of the system to resist high-frequency interference.The pattern in Figure 8b is similar to that in Figure 8a.Thus, Figure 8b was not analyzed in this paper.As shown in Figure 9, the tension fluctuations under clamping, retransmission, and laying processes were analyzed, respectively.In Figure 9, compared to the other two systems, the system without passive dancers has a larger   and a shorter   and   .However, for the tension control system in AFP, if the   and   are excessively short, the motor does not have enough time to respond, and the tension varies rapidly.In this case, the probability of laying defects increases.Meanwhile, for the system without the controller, all indicators are inferior to those with a passive roller and controller.In contrast to the other two systems, the tension control system with the passive dancers and the As shown in Figure 9, the tension fluctuations under clamping, retransmission, and laying processes were analyzed, respectively.In Figure 9, compared to the other two systems, the system without passive dancers has a larger M p and a shorter T r and T s .However, for the tension control system in AFP, if the T r and T s are excessively short, the motor does not have enough time to respond, and the tension varies rapidly.In this case, the probability of laying defects increases.Meanwhile, for the system without the controller, all indicators are inferior to those with a passive roller and controller.In contrast to the other two systems, the tension control system with the passive dancers and the controller provides a better buffer for the motor and reduces the maximum overshoot significantly, which is better for the AFP.
controller provides a better buffer for the motor and reduces the maximum overshoot significantly, which is better for the AFP.

Selection of the Passive Dancer Parameters during Different Stages
For AFP, the tension control system needs an excellent buffering mechanism, small tension fluctuations, and faster stabilization.In other words, the smaller the   ,   , and   , the better the system performance.But in practice, perfect systems are almost nonexistent.This paper proposed a new evaluative method to evaluate the performance of the established tension control system and selected the appropriate passive dancer parameters.The principle of this method is as follows: Firstly, a total score is defined for a perfect system.Secondly, the defined total score is assigned according to the different importance of each performance indicator (  ,  and   ) to the system, and then the value of each performance indicator is multiplied by its assigned score.Finally, the resulting values are added to obtain a value, and the magnitude of this value evaluates the system.In this paper, a total score of 10 was defined.And by analyzing the AFP,   ,   , and   were defined to be 5, 3, and 2 scores, respectively.It should be noted that allocating 10 points in a 5-3-2 ratio is more of a demonstration and an example.While this specific allocation of 5-3-2 may not be perfect, it reflects the requirements of a certain type of layup process.These requirements can be dynamically adjusted based on the actual needs of the process personnel.The calculation formula for evaluating the tension control system is as follows: where  is the system score, which is used to judge the system performance.   ,    , and    are the values of   ,   , and   for different systems, respectively.As for AFP, the smaller the , the better the system performance.This section analyzes the effect of different damping and elasticity coefficients on the system during different working stages in AFP.
The tension step signal mainly causes tension fluctuation during the clamping stage, and the result is shown in Figure 10. Figure 10a,b show the tension fluctuations and performance analysis of the control system under different damping conditions, respectively.From Figure 10a, it can be seen that, as  increases,   and   decrease gradually.This is due to the fact that the greater the , the lower the passive dancer deformation, and the better the ability to absorb shock and buffer.Therefore, less time is required for the system to stabilize.In particular, when  increases from 5 N • s/m to 100 N • s/m,   progressively decreases.However, as  > 100 N • s/m ,   starts to rise gradually.This is because  is excessively large, and the passive dancer moves very slowly.Thus, to make the system reach stability faster,   becomes greater.In addition, combining Figure 10a,b, it could be obtained that, when  = 5 N • s/m, 50 N • s/m, and 100 N • s/m, the system performance is poor.However, the system performs well when  = 200 N • s/m and

Selection of the Passive Dancer Parameters during Different Stages
For AFP, the tension control system needs an excellent buffering mechanism, small tension fluctuations, and faster stabilization.In other words, the smaller the T r , M p , and T s , the better the system performance.But in practice, perfect systems are almost nonexistent.This paper proposed a new evaluative method to evaluate the performance of the established tension control system and selected the appropriate passive dancer parameters.The principle of this method is as follows: Firstly, a total score is defined for a perfect system.Secondly, the defined total score is assigned according to the different importance of each performance indicator (M p , T s , and T r ) to the system, and then the value of each performance indicator is multiplied by its assigned score.Finally, the resulting values are added to obtain a value, and the magnitude of this value evaluates the system.In this paper, a total score of 10 was defined.And by analyzing the AFP, M p , T s , and T r were defined to be 5, 3, and 2 scores, respectively.It should be noted that allocating 10 points in a 5-3-2 ratio is more of a demonstration and an example.While this specific allocation of 5-3-2 may not be perfect, it reflects the requirements of a certain type of layup process.These requirements can be dynamically adjusted based on the actual needs of the process personnel.The calculation formula for evaluating the tension control system is as follows: where δ is the system score, which is used to judge the system performance.V M p , V T s , and V T r are the values of M p , T s , and T r for different systems, respectively.As for AFP, the smaller the δ, the better the system performance.This section analyzes the effect of different damping and elasticity coefficients on the system during different working stages in AFP.
The tension step signal mainly causes tension fluctuation during the clamping stage, and the result is shown in Figure 10. Figure 10a,b show the tension fluctuations and performance analysis of the control system under different damping conditions, respectively.From Figure 10a, it can be seen that, as D increases, T s and T r decrease gradually.This is due to the fact that the greater the D, the lower the passive dancer deformation, and the better the ability to absorb shock and buffer.Therefore, less time is required for the system to stabilize.In particular, when D increases from 5 N•s/m to 100 N•s/m, M p progressively decreases.However, as D > 100 N•s/m, M p starts to rise gradually.This is because D is excessively large, and the passive dancer moves very slowly.Thus, to make the system reach stability faster, M p becomes greater.In addition, combining Figure 10a,b, it could be obtained that, when D = 5 N•s/m, 50 N•s/m, and 100 N•s/m, the system performance is poor.However, the system performs well when D = 200 N•s/m and 400 N•s/m.Therefore, both D = 200 N•s/m and 400 N•s/m can be considered as the damping coefficients for the clamping stage.
progressively increases, while   and   diminish.Furthermore, combining Figure 10c,d, it could be concluded that, compared to  = 100 N/m and 300 N/m, the system performance is better when  = 600 N/m, 900 N/m, and 1800 N/m.Therefore, the elasticity factor of the clamping stage was chosen as  = 600 N/m , 900 N/m, and 1800 N/m.In particular, the damping coefficient has a more significant effect on the system than the elasticity coefficient, as shown in Figure 10.During the retransmission phase of the AFP, the tow suddenly accelerates or decelerates at the beginning and end of the transmission.Thus, in this study, the cause of tension fluctuation in the retransmission stage was simplified to a speed step.The effect of different passive dancer parameters on tension fluctuation during the retransmission stage was analyzed, and the results are shown in Figure 11.It should be noted that the tow had already an applied tension of 5 N at the beginning of the retransmission period.Hence, when analyzing the tension control system in the retransmission phase,   was not considered.Comparing Figure 11 and Figure 10, all the performance indicators of the tension system become more significant with the velocity step, which indicates that the velocity step has a greater effect on tension fluctuations than the tension step.
Figure 11a,b show the effect of different damping coefficients on the tension fluctuation during the retransfer process.From Figure 11a, it can be concluded that, as  increases,   gradually decreases, while   first reduces and then rises, which is similar to what occurs during the clamping stage.In particular, compared to the clamping phase, the damping value that makes   change from a downward to an upward trend is In the clamping phase, the elasticity coefficient also impacts the control system performance, and the results are shown in Figure 10c,d.As K increases, the distance traveled by the passive rollers and the response time of the system decrease, but the force required to reach equilibrium increases.Hence, as shown in Figure 10a, when K increases, M p progressively increases, while T r and T s diminish.Furthermore, combining Figure 10c,d, it could be concluded that, compared to K = 100 N/m and 300 N/m, the system performance is better when K = 600 N/m, 900 N/m, and 1800 N/m.Therefore, the elasticity factor of the clamping stage was chosen as K = 600 N/m, 900 N/m, and 1800 N/m.In particular, the damping coefficient has a more significant effect on the system than the elasticity coefficient, as shown in Figure 10.
During the retransmission phase of the AFP, the tow suddenly accelerates or decelerates at the beginning and end of the transmission.Thus, in this study, the cause of tension fluctuation in the retransmission stage was simplified to a speed step.The effect of different passive dancer parameters on tension fluctuation during the retransmission stage was analyzed, and the results are shown in Figure 11.It should be noted that the tow had already an applied tension of 5 N at the beginning of the retransmission period.Hence, when analyzing the tension control system in the retransmission phase, T r was not considered.Comparing Figures 10 and 11, all the performance indicators of the tension system become more significant with the velocity step, which indicates that the velocity step has a greater effect on tension fluctuations than the tension step.
sion fluctuation in the layup phase is much larger.This difference is because the tension in the layup stage suffers from tension step and velocity step effects.Thus, the discussion of tension fluctuations in the layup stage was not carried out in this paper.Only the corresponding conclusions are provided: The values of100 N • s/m, 200 N • s/m, and 400 N • s/m were considered as the damping coefficients in the layup phase.Meanwhile, 600 N/m, 900 N/m, and 1800 N/m were selected as the elasticity coefficients.The above analysis shows that the passive dancer parameters significantly affect the tension fluctuations during the clamping, retransmission, and layup processes.At the same time, compared to the retransmission and layup stages, the clamping process does   11a, it can be concluded that, as D increases, T s gradually decreases, while M p first reduces and then rises, which is similar to what occurs during the clamping stage.In particular, compared to the clamping phase, the damping value that makes M p change from a downward to an upward trend is greater in the retransmission phase.This is because the speed step has more influence on the tension fluctuation than the tension step.Therefore, more damping force must be provided for buffering to reduce the tension swing during retransmission.In addition, combining Figure 8a,b Figure 11c,d show the control system's tension fluctuations and performance analysis for the retransmission phase with different elasticity coefficients, respectively.Similar to what occurs during the clamping stage, M p gradually increases with the increase in K, while T s decreases.However, in the retransmission stage, the change in performance indicator parameters is insignificant with the K change, which indicates that the K variation could not greatly suppress the tension fluctuation during the retransmission stage.Meanwhile, by analyzing Figure 11c,d, K = 600 N/m, 900 N/m, and 1800 N/m were all considered as the elasticity coefficients of the retransmission phase.
Since the tension fluctuation law in the layup phase is similar to that in the retransmission stage, the only difference is that, compared to the retransmission phase, the tension fluctuation in the layup phase is much larger.This difference is because the tension in the layup stage suffers from tension step and velocity step effects.Thus, the discussion of tension fluctuations in the layup stage was not carried out in this paper.Only the corresponding conclusions are provided: The values of 100 N•s/m, 200 N•s/m, and 400 N•s/m were considered as the damping coefficients in the layup phase.Meanwhile, 600 N/m, 900 N/m, and 1800 N/m were selected as the elasticity coefficients.
The above analysis shows that the passive dancer parameters significantly affect the tension fluctuations during the clamping, retransmission, and layup processes.At the same time, compared to the retransmission and layup stages, the clamping process does not require a high system performance, especially T s .In summary, this section considered the damping and elasticity coefficients determined during the retransmission and layup phases as the passive dancer parameters for the AFP process.

Layup Experiments
Placement experiments were carried out using the passive dancer parameters selected by the simulation.The test setup is shown in Figure 12a.This study focused solely on investigating the influence of tension on the layup process.In this experiment, only tension was treated as a variable, while the remaining parameters were kept constant.Therefore, we assumed that any observed effects were attributed to the tension system itself rather than changes in other parameters.The experimental results show fewer layup defects are generated when the selected parameters are used for a layup, and the quality is high.As depicted in Figure 12b, the tension fluctuates excessively when D = 400 N/(m/s).Thus, it is easier to obtain pucker and edge unevenness during laying.As for D = 100 N/(m/s) and 200 N/(m/s), there is a missing tow in the laying process, which could be ignored.It can be seen from Figure 12c that the variation in the elasticity coefficient does not have much effect on the layup.When K = 600 N/m, sticking occurs during laying.And when K = 900 N/m and K = 1800 N/m, there are small wrinkles during laying.Since defects are unavoidable, the minor defects caused by changing the elasticity coefficient during the laying process were neglected.In addition, comparing the influence of the damping coefficient and elasticity coefficient on the laying process, it could be obtained that the damping coefficient has a more significant effect.(Note: The results are only for the type of prepreg selected in this paper, the AFP parameters, and the mold used for a layup.If the above conditions change, the simulation and experiments must be reconducted to determine the passive dancer parameters).
Actuators 2024, 13, 70 16 of 18 not require a high system performance, especially  .In summary, this section considered the damping and elasticity coefficients determined during the retransmission and layup phases as the passive dancer parameters for the AFP process.

Layup Experiments
Placement experiments were carried out using the passive dancer parameters selected by the simulation.The test setup is shown in Figure 12a.This study focused solely on investigating the influence of tension on the layup process.In this experiment, only tension was treated as a variable, while the remaining parameters were kept constant.Therefore, we assumed that any observed effects were attributed to the tension system itself rather than changes in other parameters.The experimental results show fewer layup defects are generated when the selected parameters are used for a layup, and the quality is high.As depicted in Figure 12b, the tension fluctuates excessively when D = 400 N/(m/s).Thus, it is easier to obtain pucker and edge unevenness during laying.As for D = 100 N/(m/s) and 200 N/(m/s), there is a missing tow in the laying process, which could be ignored.It can be seen from Figure 12c that the variation in the elasticity coefficient does not have much effect on the layup.When K = 600 N/m, sticking occurs during laying.And when K = 900 N/m and K = 1800 N/m, there are small wrinkles during laying.Since defects are unavoidable, the minor defects caused by changing the elasticity coefficient during the laying process were neglected.In addition, comparing the influence of the damping coefficient and elasticity coefficient on the laying process, it could be obtained that the damping coefficient has a more significant effect.(Note: The results are only for the type of prepreg selected in this paper, the AFP parameters, and the mold used for a layup.If the above conditions change, the simulation and experiments must be reconducted to determine the passive dancer parameters).

Conclusions
This paper establishes a precise tension fluctuation model to describe the pattern of tension variation during AFP.The tests confirm the effectiveness of this model, and the Bode diagram illustrates the system's stability.Moreover, the superiority of the tension control system with passive rollers is demonstrated by simulation analysis.After that, the tension fluctuations under different elasticity coefficients and damping are obtained for each phase in the AFP.The evaluation method developed in this article determines the

Conclusions
This paper establishes a precise tension fluctuation model to describe the pattern of tension variation during AFP.The tests confirm the effectiveness of this model, and the Bode diagram illustrates the system's stability.Moreover, the superiority of the tension control system with passive rollers is demonstrated by simulation analysis.After that, the tension fluctuations under different elasticity coefficients and damping are obtained for each phase in the AFP.The evaluation method developed in this article determines the passive dancer parameters with a better system performance.Furthermore, placement experiments are carried out using the determined passive dancer parameters.
(1) The model's validity was verified by developing wire-feeding and clamping experiments.Also, the stability and superiority of the control system with a passive dancer were illustrated by simulation analysis.The system provides a better buffer for the motor and reduces the maximum overshoot significantly.(2) The tension fluctuation at each stage of AFP can be derived from the precise tension control model.By analyzing the tension fluctuation with different passive dancer parameters, it can be concluded that the damping has more influence on the tension fluctuation than the elasticity coefficient.Meanwhile, by appropriately increasing the damping and elasticity coefficients, the system performance is improved, and the tension fluctuation is reduced, thus enhancing the layup quality.(3) By analyzing the tension fluctuations at each stage, the elastic coefficients of the passive dancer in the AFP process were determined to be 600 N/m, 900 N/m, and 1800 N/m, and the damping was 100 N•s/m, 200 N•s/m, and 400 N•s/m.Furthermore, it was demonstrated that using the determined passive dancer parameters for laying could improve the laying quality significantly.

Figure 1 .
Figure 1.(a) Schematic diagram of the automatic fiber placement.(b) Schematic diagram of the tension control system.

Figure 1 .
Figure 1.(a) Schematic diagram of the automatic fiber placement.(b) Schematic diagram of the tension control system.

Figure 2 .
Figure 2. (a) Schematic diagram of the pay-off axis model.(b) Schematic diagram of the simplified velocity difference-tension model.(c) Schematic diagram of the passive dancer model.

Figure 2 .
Figure 2. (a) Schematic diagram of the pay-off axis model.(b) Schematic diagram of the simplified velocity difference-tension model.(c) Schematic diagram of the passive dancer model.

Figure 3 .
Figure 3. (a) Flow chart of the system control without passive dancers.(b) Flow chart of the system control with passive dancers.

Figure 3 .
Figure 3. (a) Flow chart of the system control without passive dancers.(b) Flow chart of the system control with passive dancers.

Figure 4 .
Figure 4. (a,b): Schematic diagram of the tension step signal.(c,d): Schematic diagram of the velocity perturbations.

Figure 4 .
Figure 4. (a,b): Schematic diagram of the tension step signal.(c,d): Schematic diagram of the velocity perturbations.

Figure 6 .
Figure 6.Flow chart of the algorithm for the tension control.

Figure 6 .
Figure 6.Flow chart of the algorithm for the tension control.

Figure 7 .
Figure 7.The validity verification of the model.(a) Tension step model validation; (b) Velocity step model validation; (c) Velocity sinusoidal input model validation.

Figure 7 .
Figure 7.The validity verification of the model.(a) Tension step model validation; (b) Velocity step model validation; (c) Velocity sinusoidal input model validation.

Figure 9 .
Figure 9. (a) Tension fluctuations during the tension step.(b) Tension fluctuations during the velocity step.(c) Tension fluctuations during both steps.

Figure 9 .
Figure 9. (a) Tension fluctuations during the tension step.(b) Tension fluctuations during the velocity step.(c) Tension fluctuations during both steps.

Figure 10 .
Figure 10.(a) Tension fluctuations in the clamping phase with different damping coefficients ( = 300 N/m).(b) Performance analysis of the control system with different damping coefficients.(c) Tension fluctuations in the clamping phase with different elasticity coefficients ( = 200 N • s/m).(d) Performance analysis of the control system with different elasticity coefficients.

Figure 10 .
Figure 10.(a) Tension fluctuations in the clamping phase with different damping coefficients (K = 300 N/m).(b) Performance analysis of the control system with different damping coefficients.(c) Tension fluctuations in the clamping phase with different elasticity coefficients (D = 200 N•s/m).(d) Performance analysis of the control system with different elasticity coefficients.

Figure 11 .
Figure 11.(a) Tension fluctuations in the retransmission phase with different damping coefficients ( = 300).(b) Performance analysis of the control system with different damping coefficients.(c) Tension fluctuations in the retransmission phase with different elasticity coefficients ( = 100).(d) Performance analysis of the control system with different elasticity coefficients.

Figure 11 .
Figure 11.(a) Tension fluctuations in the retransmission phase with different damping coefficients (K = 300).(b) Performance analysis of the control system with different damping coefficients.(c) Tension fluctuations in the retransmission phase with different elasticity coefficients (D = 100).(d) Performance analysis of the control system with different elasticity coefficients.

Figure
Figure11a,b show the effect of different damping coefficients on the tension fluctuation during the retransfer process.From Figure11a, it can be concluded that, as D increases, T s gradually decreases, while M p first reduces and then rises, which is similar to what occurs during the clamping stage.In particular, compared to the clamping phase, the damping value that makes M p change from a downward to an upward trend is greater in the retransmission phase.This is because the speed step has more influence on the tension fluctuation than the tension step.Therefore, more damping force must be provided for buffering to reduce the tension swing during retransmission.In addition, combining Figure8a,b, it also can be obtained that the system has a low resistance to tension fluctuation when D = 5 N•s/m and 50 N•s/m.And when D = 100 N•s/m, 200 N•s/m, and 400 N•s/m, the overall performance of the system is better.Therefore, D = 100 N•s/m, 200 N•s/m, and 400 N•s/m were all considered as the damping values for the retransmission phase.Figure11c,d show the control system's tension fluctuations and performance analysis for the retransmission phase with different elasticity coefficients, respectively.Similar to what occurs during the clamping stage, M p gradually increases with the increase in K, while T s decreases.However, in the retransmission stage, the change in performance indicator parameters is insignificant with the K change, which indicates that the K variation could not greatly suppress the tension fluctuation during the retransmission stage.Meanwhile, by analyzing Figure11c,d, K= 600 N/m, 900 N/m, and 1800 N/m were all considered as the elasticity coefficients of the retransmission phase.Since the tension fluctuation law in the layup phase is similar to that in the retransmission stage, the only difference is that, compared to the retransmission phase, the tension fluctuation in the layup phase is much larger.This difference is because the tension in the layup stage suffers from tension step and velocity step effects.Thus, the discussion of tension fluctuations in the layup stage was not carried out in this paper.Only the corresponding conclusions are provided: The values of 100 N•s/m, 200 N•s/m, and 400 N•s/m were considered as the damping coefficients in the layup phase.Meanwhile, 600 N/m, 900 N/m, and 1800 N/m were selected as the elasticity coefficients.
Figure11a,b show the effect of different damping coefficients on the tension fluctuation during the retransfer process.From Figure11a, it can be concluded that, as D increases, T s gradually decreases, while M p first reduces and then rises, which is similar to what occurs during the clamping stage.In particular, compared to the clamping phase, the damping value that makes M p change from a downward to an upward trend is greater in the retransmission phase.This is because the speed step has more influence on the tension fluctuation than the tension step.Therefore, more damping force must be provided for buffering to reduce the tension swing during retransmission.In addition, combining Figure8a,b, it also can be obtained that the system has a low resistance to tension fluctuation when D = 5 N•s/m and 50 N•s/m.And when D = 100 N•s/m, 200 N•s/m, and 400 N•s/m, the overall performance of the system is better.Therefore, D = 100 N•s/m, 200 N•s/m, and 400 N•s/m were all considered as the damping values for the retransmission phase.Figure11c,d show the control system's tension fluctuations and performance analysis for the retransmission phase with different elasticity coefficients, respectively.Similar to what occurs during the clamping stage, M p gradually increases with the increase in K, while T s decreases.However, in the retransmission stage, the change in performance indicator parameters is insignificant with the K change, which indicates that the K variation could not greatly suppress the tension fluctuation during the retransmission stage.Meanwhile, by analyzing Figure11c,d, K= 600 N/m, 900 N/m, and 1800 N/m were all considered as the elasticity coefficients of the retransmission phase.Since the tension fluctuation law in the layup phase is similar to that in the retransmission stage, the only difference is that, compared to the retransmission phase, the tension fluctuation in the layup phase is much larger.This difference is because the tension in the layup stage suffers from tension step and velocity step effects.Thus, the discussion of tension fluctuations in the layup stage was not carried out in this paper.Only the corresponding conclusions are provided: The values of 100 N•s/m, 200 N•s/m, and 400 N•s/m were considered as the damping coefficients in the layup phase.Meanwhile, 600 N/m, 900 N/m, and 1800 N/m were selected as the elasticity coefficients.

Table 1 .
The required parameters for the simulation.

Table 1 .
The required parameters for the simulation.