Development of a Delamination Fatigue Testing Machine for Composite Materials

: Experimental studies on fatigue behavior are usually conducted on servo-hydraulic testing machines that are expensive and have high maintenance costs. In this work, a much simpler testing machine was developed, intended mainly for delamination fatigue tests on composite materials. After a literature review on the methods and parameters of such tests, the machine was designed, and its parts manufactured and assembled into a fully operational testing machine. Additionally, the electrical components and the control and data acquisition software were also developed and implemented. Finally, several mode II delamination fatigue tests were conducted using the end-notched flexure test. The results were consistent with the well-known Paris law, which for composite materials relates the crack propagation rate to the strain-energy release rate range. Therefore, the developed machine seems to be an excellent alternative to the highly costly testing machines.


Introduction
High-performance composite materials are nowadays used in a wide range of structural applications, being a consolidated alternative to metals [1]. In addition to their excellent mechanical properties and low specific weight, polymer matrix composites enable a more flexible design stage relative to conventional materials. It is also possible to manufacture complex shapes in a wide range of dimensions. These advantages have led to increased applications in various industries, such as aeronautical, automotive, and sports [2]. However, in high-performance laminated composites, layers may separate, a failure mode that is known as delamination.
Delamination occurs mainly due to the low strength of the thin resin layer that binds the laminate layers [1,2]. In addition, delamination occurs preferentially in regions with significant interlaminar stresses, especially between layers with different orientations. Delamination may also occur due to transverse cracks in the layers. Delamination resistance is analyzed in the scope of fracture mechanics. Accordingly, there are three pure modes of delamination (j = I, II, III), which are characterized by critical strain-energy release rates Gjc [3] measured in different tests [4].
Fatigue testing involves the application of cyclic loads below the strength of materials. Such loads cause the formation and propagation of damage, which, in the end, can promote failure. In homogeneous and isotropic materials, fatigue failure is usually dictated by the initiation and propagation of a single crack, while in composite materials there is an accumulation of different types of damage, e.g., transverse cracking, delamination, and longitudinal cracking. Nevertheless, delamination under fatigue has been considered relevant enough to warrant specific testing procedures. In fact, delamination under fatigue loading has been the subject of several studies [4][5][6][7][8][9][10][11][12][13]. In those studies, cyclic loads have been applied to specimens adopted by the standards for fracture toughness testing.
Fatigue testing is a time-consuming procedure, which often requires long periods of continuous testing, and gives rise to costly experimental programs. Therefore, it is of considerable interest to have dedicated low-cost machines for fatigue testing. For example, McAlorum et al. [13] recently developed a low-cost small-scale machine (25 kN) for the custom fatigue testing of structural beams. The authors carried out an analysis comparing the costs of a dedicated machine with those of outsourcing. The results indicated that the custom machine was the preferred option for long-term testing of at least 373 h.
The development of several fatigue testing machines (FTMs) has been developed and reported in the literature, e.g., the closed loop servo-hydraulic FTM developed by Lee [14]. This testing machine consists of a structural bending load frame, a structural torsion load frame, a hydraulic system, and a control system. The hydraulic system allows bending loadings up to 490 kN. Such value makes this machine appropriate for large-scale testing, but there are high costs involved. Therefore, FTMs are typically designed for lower loading capacities, unless the specificity of the test and the lack of standard equipment justify custom-made solutions. More recently, Pach et al. [15] developed an FTM in order to study the fatigue behavior of fibre-reinforced polymer (FRP) composites under constant amplitude load-cycles and a wide range of load ratios. Pach et al. [15] highlighted the need for conducting the tests at low frequencies, in order to avoid temperature rises resulting from viscoelastic dissipation of the polymer matrix. Consequently, the test programs demand much more time than those involving metals. In contrast to Lee [14], Pach et al. [15] also adopted a low-cost solution instead of conventional servo-hydraulic testing machines. The approach consisted of applying axial loads (10 kN max) by using an "seesaw" motion created by a set of helical springs compressed by an electric motor driven mechanism [16]. Nevertheless, some design problems were identified, such as lack of correct grip [16]. Vincent et al. [16] designed a much simpler FTM based on an electric motor drive. However, this simple and low-capacity solution was very limited in terms of maximum load (300 N). This is a significant constraint for fatigue testing of composites.
This paper describes a fatigue testing machine (FTM) developed for fatigue testing of composite materials. Its main advantages are the low costs and easy adaptability to the standard delamination test specimens.

Fatigue Testing Machine
The machine overall 3D drawing is shown in Figure 1. Basically, it consists of a rigid frame and a movable cross-head. The rigid frame has a top and a bottom table, which are connected to each other by four support rods, while the movable cross-head is driven by two sets of nut/spindle (TR24 × 5) along a maximum 600 mm course. The two spindles that drive the cross-head are mounted and connected to the machine tables by means of bearings. One of the spindles is coupled to the stepper motor. The other spindle is simultaneously driven by a timing-belt transmission. To adjust the belt tension, there is also a stretcher mounted on the machine's top table. The gear ratio between the motor and the spindles is 1:1. The load cell bolted to the movable crosshead allows the attachment of loading rigs needed for various types of tests. The lower table has a 300 mm wide open workspace and can accommodate various types of fixtures. In the design of the FTM, a BSHB3913-H hybrid stepper motor was chosen [17]. Hybrid motors have high precision and high torque, and can be configured for small angles. Stepper motors are useful in applications with accurate incremental actions [18].
Regarding the load cell, a Ningbo Xinlan XL 8111 with 10 kN capacity was used. This is an Stype tensile/compression load cell. The load cell support system made of a nonconductive thermoplastic was machined in a Mikron VCE 500 CNC machine centre. Other parts of the FTM were also machined in this centre.
A DM3522M driver [17] was used to power the stepper motor of the FTM. The driver's main purpose is to amplify the current received from the microcontroller (Arduino Uno Rev3 based on the ATmega328) that controls the motor. In addition, the driver's characteristics make it easier to control the rotor [17]. The driver has a built-in panel to program parameters, such as micro or single step, number of steps per rotation, and direction of rotation, among others.
The power and amplification of the signal from the load cell were carried out by an integrated circuit, model INA 125 [19]. This model is a low-power, high-accuracy instrumentation amplifier with a precision voltage reference.
To correct any possible errors in the number of steps of each cycle, a Sharp GP1S53V000F sensor was used to evaluate and correct any loss of steps. This sensor was employed instead of other systems such as encoders, because it is a cheaper and easier solution to implement. This sensor has an infrared light-emitting diode (LED) and a phototransistor detector. It was placed in a housing attached to the movable beam, as shown in Figure 2. This adjustable system is connected to an aluminum beam, which is rigidly connected to the tables.  The control can be divided into two main sectors: a first sector with Arduino, driver, stepper motor, and sensor; and a second sector with the power supply, load cell, and amplifier. In the first sector, the computer with the control software sends via USB input to Arduino, which processes the commands and transmits to the driver the operation sequence of the stepper motor. The zero sensor sends to the Arduino the maximum and minimum values, adjusting them in the case of missing steps. In the second sector, the load cell is powered with 10 V. The data sent from the load cell to the Arduino are amplified by the INA125 instrumentation amplifier. These data are sent to the computer by the Arduino, displaying it graphically by another routine developed in Processing 1.5.1 and recording it in a text file for later processing. All cables from the load cell, sensor, and Arduino were connected to an electrical panel, which includes buttons to rotate the engine in a specific direction, a switch for manual/automatic control, a speed control potentiometer, and an emergency button. Additionally, ground connections were created for the machine's structure and electric circuit. In the first measurements, significant noise was detected. A GoldStar Oscilloscope OS-9020A 20 MHz showed that the noise was caused essentially by the stepper motor and, to a lesser extent, by the driver. The noise was of the order of 1 V and had a low frequency. In order to reduce noise to acceptable levels, shielded cables were used to insulate against electromagnetic interference. Additionally, a low-pass filter was used between the load cell and the INA125 instrumentation amplifier.
Finally, the load cell was successfully calibrated in a Shimadzu AG 50 kN testing machine. A maximum load of 1 kN was used for calibration, since this value was above the maximum load expected in mode II fatigue tests. Figure 4 shows the completely assembled FTM. Table 1 presents the key performance parameters of the FTM developed in this work, comparing it with traditional ones.  1 A 10 kN load cell was employed, but it is only used up to 5 kN in order to guarantee durability in fatigue cycles.
As seen in Table 1, the FTM has the high accuracy levels needed for fatigue tests on composite materials. Its results were validated by using simultaneously with Arduino the NI CDAQ-9178 data acquisition board, manufactured by National Instruments (mV precision). This equipment is typically used for highly precise measurements, and thus was employed for calibration of the FTM measurement systems. It is also worth mentioning that the closed circuit used for the stepper motor and the zero sensor enables precise control of the displacement. In addition, the machine compliance was measured and taken into account when imposing or computing displacements.
The fatigue strength of the FTM components experiencing significant stress levels was carefully approached during the design process. Nevertheless, bearings and timing belt will have to be replaced in due time, which depends on the recorded usage.
After the aforementioned calibrations and validations, the mode II fatigue delamination tests described below provided further assessment of the FTM performance and practical usefulness.

Mode II Fatigue Delamination
The first application of the developed testing machine concerned the mode II fatigue delamination of a carbon/epoxy material. Considerable research has been reported on this issue, especially on measurement and modeling of crack propagation rates [8,20]. The main goal of the experimental work conducted was precisely to measure crack propagation rates within the framework of the well-known Paris law [21], which relates the crack propagation rate to the stress intensity factor range ΔK. Actually, strain-energy release rates G are used for characterizing delamination resistance of composites, and thus, in the scope of mode II loadings, Paris law has been expressed as where a is the crack length, N the number of load cycles, ΔGII the strain-energy release rate range, and D and m empirical parameters. Nevertheless, it is well known that GII = ψ(KII) 2 , where ψ is a function of the material elastic constants. Therefore, in order to ensure the correctness of loading superposition and similitude with ΔK of the original Paris law, one should define ΔGII = (√GIImax− √GIImin) 2 [20]. The comprehensive review by Pascoe et al. [21] showed that such definition of ΔG provides more consistent results. Unidirectional specimens were prepared from a high-strength carbon fiber (T300)/toughened epoxy Pre-preg (reference HS160REM) supplied by Seal. The specimen preparation procedures adopted are detailed in previous studies [22], in which the main elastic properties were found to have the following average values: longitudinal Young's modulus E1 = 130 GPa, transverse Young's modulus E2 = 8.2 GPa, in-plane shear modulus G12 = 4.1 Gpa, and longitudinal Poisson's ratio ν12 = 0.27. The mode II test selected was the well-known end-notched flexure (ENF) ( Figure 5). The specimens had nominal thickness 2h = 3 mm, width b = 20 mm, and a starter delamination created by a 13 μm thick polytetrafluoroethylene (PTFE) insert. In order to check that no spurious dynamic, friction, or specimen sliding effects were relevant, two combinations of starter crack length (a0) and span (2L) were used: a0 = 40 mm, 2L = 125 mm; and a0 = 32 mm, 2L = 100 mm. The total length of the specimens was 200 mm, so that several tests could be carried out on each specimen: a first quasi-static test to measure the critical strain-energy release rate GIIc, fatigue tests, and a final quasi-static test. The data reduction method employed was the effective crack method (ECM) in the form presented in [22,23]. It consists of computing an effective crack length a e = √ 8E 1 bh 3 C f -2L 3 3 3 (2) from experimental values of the longitudinal flexural modulus E1 and specimen flexural compliance Cf. The latter is obtained by discounting the contribution of transverse shear to the total experimental compliance C by where G13 is the through-thickness shear modulus. The strain-energy release rate is then obtained from II = 9 2 2 16 2 1 ℎ 3 (4) As discussed in [22], the main advantages of the ECM are twofold: simplicity, since it avoids the crack length monitoring known to be particularly difficult in mode II tests; and accuracy, as it intrinsically takes into account the extensive Fracture Process Zone (FPZ) associated to mode II delamination [24]. Under steady-state crack propagation conditions, the FPZ remains constant and thus dae/dN ≈ da/dN. On the other hand, during the initial stages of quasi-static and fatigue tests, the ECM is unable to distinguish FPZ formation from actual crack propagation. Therefore, in quasi-static tests, GIIc was defined at the plateau of the GII-ae curves. Similarly, the GII-dae/dN data for the first Δae values were discarded, which also avoided the effects of the previous test the specimens had undergone. Evidently, this approach demanded that fatigue damage was confined to the close vicinity of the crack tip. As seen below, the results proved that this hypothesis was valid.
Fatigue tests were carried out at a low 2 Hz test frequency in order to avoid viscoelastic effects. Displacement control was adopted rather than load control, which in the ENF specimen would lead to very fast crack propagation and produce very few data points. The stress ratio R = δmin/δmax was varied from 0.26-0.36, while GIImax/GIIc = 0.32-0.66, so that tests were not excessively long, e.g., typically around 10 5 cycles. Another advantage of displacement control was that GIImin and GIImax varied mildly with crack propagation, i.e., less than 3% for a Δae = 1 mm increment, except very close (2 mm) to the load-point, a region that the crack should not be allowed to reach. This enabled a very robust evaluation of dae/dN, which is particularly important, given the unavoidable noise in raw data.

Results
In order to verify the correct functioning and accuracy of the developed FTM, an experimental study was carried out on mode II delamination. Figure 6 depicts typical GII-ae curves measured in quasi-static tests. Tests conducted from the insert did not reveal any R-curve effect. In turn, because of fatigue damage, a very mild R-curve was seen in tests conducted from the fatigue pre-crack. Nevertheless, there was no significant difference in the maximum GIIc of specimens tested from the insert (average 1.21 kJ/m 2 , 6.9% standard deviation) or from the fatigue pre-crack (average 1.12 kJ/m 2 , 4.5% standard deviation). This shows that fatigue damage was limited to the FPZ and did not affect the other regions of the specimen, thereby supporting the data reduction approach adopted (Section 2.2). In order to take into account the scatter in fracture toughness values, the maximum R-curve plateau GIIc from the fatigue pre-crack was used to normalize fatigue results. The Paris law form adopted was thus As expected from the aforementioned specimen characteristics, ae versus N plots presented several quasi-linear regions, despite some noise and scatter (Figure 7). In other words, dae/dN remained nearly constant throughout considerable parts of the test, hence favoring the consistency and representativeness of GII-dae/dN data under steady-state propagation. It is worth reminding that, as discussed above, GII-dae/dN data for the first Δae values (typically Δae = 1-1.5 mm) were discarded.   (5) with D = 16.7 mm/cycle and m = 6.9, given the high 0.967 correlation coefficient. Moreover, the m = 6.9 of Equation (5) obtained was similar to the values reported in other studies [11,25]. In fact, Matsubara et al. [11] obtained m = 6.5 for a unidirectional glass/epoxy material tested under a stress ratio R = 0.3, which was within the range here adopted. Matsubara et al. [11] employed a form of the Paris law similar to Equation (5), but used the threshold GII range, ΔGIIth, instead of GIIc. Using the ΔGIIth, GIIc, and the Paris law parameters reported in [11], one concludes that, if processed with Equation (5), the results by Matsubara et al. [11] give D = 7.1 mm/cycle. This value is lower than the one measured here, with the obvious explanation being the different materials employed. Tanaka et at. [26] studied the mode II fatigue delamination of a unidirectional carbon/epoxy material under different R-ratios. They [26] adopted the original form of the Paris law, i.e., they used Equation (1) with ΔKII instead of ΔGII. Considering that ΔGII = ψ(ΔKII) 2 when ΔGII = (√GIImax − √GIImin) 2 , we conclude that Equation (5) applied to the R = 0.2-0.5 results of Tanaka et al. [26] yields m = 6-9.5, an interval where the present m = 6.9 fits quite well. Moreover, using the ΔKIIth and KIIc measured in [26], one finds that Equation (5) gives D = 2.8-393 mm/cycle, an interval which, again, includes the present D value.

Figure 8.
Effective crack propagation rate versus the strain-energy release rate range experimental data points. The fitting with Paris law-type Equation (5) is also shown.
A broader scope study of the mode II fatigue delamination of the carbon/epoxy composites can be found in [27]. Here, it is important to stress out that the testing machine developed proved to be suitable for mode II delamination fatigue testing of composites.

Conclusions
The present work presents the entire development of a fully operational FTM, including its control and data acquisition software. The design, assembly, control system, and software for test data acquisition are also described. The solution here presented is simple, robust, and has low manufacturing and maintenance cost, especially when compared to servo-hydraulic testing machines. Moreover, it has a significant loading capacity relative to other similar machines described in the literature.
Additionally, in order to assess the performance of the FTM, several mode II delamination fatigue tests were conducted. The results were consistent with the well-known Paris law and compared favorably with results available in the literature, thereby confirming the suitability and usefulness of the FTM. Preliminary work devoted to mode I fatigue has already been carried out with the FTM, confirming its advantageous characteristics.