Research Article  Open Access
Xueshu Liu, Yuxing Yang, Yiqi Wang, Yongjie Bao, Hang Gao, "Effects of Hole Perpendicularity Error on Mechanical Performance of SingleLap DoubleBolt Composite Joints", International Journal of Polymer Science, vol. 2017, Article ID 2790198, 11 pages, 2017. https://doi.org/10.1155/2017/2790198
Effects of Hole Perpendicularity Error on Mechanical Performance of SingleLap DoubleBolt Composite Joints
Abstract
Static tensile experiments and numerical simulations were carried out to study effects of hole perpendicularity error on mechanical performance of singlelap doublebolt composite joints. Hole tilting angle, varying from 1 to 4 degrees, and hole tilting direction to account for anisotropic properties of composite material were investigated. Progressive damage model (PDM) based on ChangLessard type criterion with an extension by Olmedo was created, which involved seven failure modes and made material properties be a direct function of predefined field variables. The model was implemented in ABAQUS/Standard using a UMAT subroutine. Good agreement was found when comparing numerical simulation results with experimental outcomes. In addition, the results demonstrate that, with increasing of hole tilting angle, damage is prone to arise and the load path of the composite bolted joints alters with changing of hole tilting direction, which result in severe stress concentration around the edge of hole and joint failure in advance.
1. Introduction
Designable combinations of fiber and matrix make the composite a versatile material. Their high strength coupled with light weight leads to their use wherever structural efficiency is at a premium. However, joints are always necessary to manufacture a large composite structure, among which bonded joints become increasingly common in automotive and marine field due to their ability to improve structural performance and help weight reduction [1–3]. Nevertheless, when loadcarrying capacity is concerned, bolted joints are the most important methods of composite assembling especially in aerospace industry because of their high loadcarrying capacity and facility to assemble, disassemble, and repair. For such a joint, it is essential to machine the composite laminates with the utilization of, for example, drills, and thus the structure integrity of composite laminate is broken, which results in stress concentration around the edge of hole and a significant reduction of the loadcarrying capacity. Therefore, the bolted joint is a key point of the composite structural design.
The main purpose of studies on composite bolted joints is to determine the effect of various parameters on joint strength. Theoretical results are also used to serve the improvement of bolted joint design and exploit the full potential of CFRP. Unsustainable and conservative design can therefore be avoided with methods that have been benchmarked against satisfactory data. The parameters that have been studied include geometrical factors [4, 5], fastener parameters [6], material properties [7, 8], application factors [9], and manufacturing process [10–12].
The process factors determine the structural parameters of composites, the physicalmechanical properties, and the residual stresses produced by assembly process. Typical process factors include the bolt hole clearance effects, the accuracy and interchangeability of the joints, and the machining quality of the holes. Many factors have been studied by scholars except for hole perpendicularity error. In this paper, hole perpendicularity error is taken into consideration to enrich the study of composite bolted joints. Finite element analysis (FEA) and experiments were carried out and compared to assess the mechanical performance of composite joints.
2. Problem Description
Hole perpendicularity error is a kind of manufacturing defects, which results in stress concentration around hole edge area. In the case of singlelap singlebolt composite joint, the appearance of hole perpendicularity error will lead to bearing capability decrease of the bolted joint [12]. However, composite parts are often assembled with multiple bolts in engineering applications. Bolt holes with perpendicularity errors will affect each other, which eventually influences the mechanical performance of the joint. In this paper, by designing and carrying out tensile test for singlelap doublebolt composite joints, the bearing capacity and failure mode of the joints were obtained, and the effects of hole perpendicularity error on the joint performance were investigated.
Considering the anisotropic properties of CFRP, two parameters are used to characterize the hole perpendicularity error, tilting direction and tilting angle β, as shown in Figure 1. Project the hole on a plane, say , where is the longitudinal direction of the specimen. The initial side of is on the positive axis and the terminal one is on OC′, namely, the projection of bolt hole center axis (Figure 1(a)), with positive angles representing rotations toward the positive axis. In the case where the plate thickness is constant, the angle between the bolt hole center axis, OC′, and the theoretical one, OC, is used to characterize the magnitude of hole perpendicularity error (Figure 1(b)).
(a)
(b)
Only tilting directions and are discussed to reduce the parameters of the study. The influence of tilting direction for doublebolt composite joints can be classified into 6 different patterns after taking into account the symmetry of the structure. The diagrammatic sketches are demonstrated in Figure 2, where red lines indicate bolts and the black arrows represent the laminates and their loading directions.
(a) P0
(b) P1
(c) P2
(d) P3
(e) P4
(f) P5
3. Methodology
3.1. Specimens and Experimental Procedures
The specimens used in this study were manufactured from IMS194/CYCOM9772 carbon fiber/epoxy prepreg (nominal ply thickness is 0.188 mm) and designed to fail initially in the bearing mode with symmetric stacking sequence of [−45/0/45/90/0/90/45/0/−45/45. The material properties provided by the manufacturer are shown in Table 1. Protruding head bolts made of SCM35 alloy steel are used as fasteners, whose elastic modulus and Poisson ratio are, respectively, 190 GPa and 0.3. The elasticity and Poisson’s ratio of nuts and gaskets are the same as those of bolts. The plastic behavior of the SCM35 alloy steel is given in Table 2.
 
value is set according to the transversely isotropic assumption and engineering experience. 

Singlelap doublebolt composite joints were prepared according to the ASTM standard D5961/D5661 M13 with widthtohole diameter ratio and edge distancetohole diameter ratio [13]. Standard specimen configurations using fixed values of test parameters are shown in Figure 3.
According to the ASTM D5961 standard, the static tensile tests were performed on a WDW100 testing machine in displacement control at a rate of 1.0 mm/min. The experimental setup is shown in Figure 4. The displacement was recorded using a YSJ505ZC extensometer with a gauge length of 50 mm and measure range of 12.5 mm. The loading procedure is not stopped until the load has dropped off about 30% after it reaches the maximum, or displacement exceeds the extensometer range. A preload torque of 7.1 Nm was applied to both bolts. The configurations of test are shown in Table 3 and the group is marked as “Pi  β ,” where “Pi” indicates the pattern as shown in Figure 2, means hole tilting angle, and , if there is any, represents the maximum load applied to the joint during test. Otherwise, the joint will be loaded until its failure. Note that 30 kN is about 80% of the failure load of the singlelap doublebolt composite joint.

3.2. Numerical Model
To investigate the effects of hole perpendicularity error on mechanical behavior of singlelap doublebolt composite joints, a nonlinear threedimensional finite element model was constructed with 8node solid element C3D8R using the commercial finite element software, ABAQUS [14], as illustrated in Figure 5.
To reduce the computational cost, the grip area of the composite laminate was neglected and the laminate was modelled with 5 elements in thickness direction, which avoids excessive widthtothickness ratio of the element [7].
Contact interactions between all surfaces were modelled by surfacebased contact considering the linear penalty method. Seventeen contact interactions with finite sliding formulation were defined to simulate the interfaces between any two parts. The stiffness of master surface was relatively larger and the meshes were designed to be coarser than those of slave surfaces. Default contact control and contact stabilization were defined for each contact pair to improve convergence. Constant value of friction coefficient was set to 0.3 [8].
Boundary conditions were applied to reference points, which defined coupling constraints to node sections in the FEA model as shown in Figure 5. Section A was limited by three displacement degrees of freedom, while displacement load was applied to section B in direction with the other two displacement degrees of freedom fixed. A preload force of 6 kN is applied to the bolt shank by “bolt load” function.
There were 42,272 elements in total in the finite element model. Considering that the region around the hole is the key area of stress concentration, the size of refined meshes was approximately 0.36 × 0.77 mm in the vicinity of the hole in order to ensure the accuracy. Coarser meshes were adopted further to reduce computational cost. The element size of bolts and washers was set up corresponding to the size of the contact areas of the laminates.
The progressive damage method has attracted extensive attention in analysis and design of composite structures, which effectively simulates damage formation, evolution, and ultimate failure process of complex composite structures. Unlike isotropic materials, the destruction of composite materials is a gradual process. The progressive damage method is based on damage mechanics, which regards the failure of composite materials as damage. The material degradation model was used to describe the mechanical behavior of damaged material.
The progressive damage analysis includes the stress solving of the nonlinear finite element equilibrium equation, the prediction of damage initiation, damage evolution, and the stiffness degradation based on the material failure criterion. Stress analysis is mainly carried out by finite element software and has evolved from twodimensional model into threedimensional model. There are diverse failure criteria to predict the damage initiation including Chang criteria [15], TsaiWu criteria [16], Hashin type criteria [17], and Puck criteria [18]. In this paper, a damage model considering the effect of nonlinear shear stressstrain behavior was established based on extended ChangLessard type criteria [19, 20]. Seven basic microscopic failure modes in composite were detected and distinguished by utilizing the failure criterion, as illustrated in (1)–(7).
Tensile fiber failure, : Compressive fiber failure, :Fibermatrix shearout failure, : Matrix cracking failure, :Matrix crushing failure, :Delamination onset in tension, :Delamination onset in compression, :where are the scalar components of the stress tensor; , , , , , and are normal strength, in which subscripts “” and “” denote tensile and compressive status, respectively; and are shear strength. is a constant parameter which was determined experimentally and set as 2.44 × 10^{−8}.
The composite laminates are subjected to gradual increasing load. Once the stress matches the equations of failure criteria, damage occurs. The mechanical properties of the material elements in the damaged area will be degenerated with a sudden drop and the degradation rules are given in Table 4.

There are no acknowledged degradation criteria of material property. Most of them are determined based on experiments and experience. The sudden drop degradation model is relatively easy to implement and thus is widely used in the progressive damage analysis of composites. Tan [21], F.K. Chang and K.Y. Chang [22], Camanho and Matthews [23], McCarthy et al. [24], and Tserpes et al. [25] put forward different degradation models, which match the specific failure criteria, respectively.
4. Results
4.1. Model Validation
Figure 6 shows the load–displacement curves of experimental results and FEA one of P0. It can be found that all results show the same tendency. In addition, the ultimate bearing load is about 38 kN, which proves that the numerical result is reliable when bearing load is concerned.
As shown in Table 5, the average ultimate loads of P0, P14, and P24 are 38387 N, 37737 N, and 37745 N, respectively. The difference is less than 2%, which indicates that when hole perpendicularity error is less than 4°, it has limited influence on ultimate load of singlelap doublebolt joints. The numerical prediction for ultimate load of the joint is 35900 N. The errors of the ultimate load and stiffness between the experimental outcomes and FEA results are about 6.5% and 12.5%, which demonstrates that the PDM can be used to predict the failure loads of composite joints.

Figure 7 shows the experimental damage morphology of the top laminate of P1. It can be seen that the laminate mainly experienced bearing damage. Fiber breakage mainly happened around the bolt hole and in the transvers direction of the laminate. In addition, the bearing failure of fiber and matrix mainly occurred on the contact surface of the bolt shank and the hole. Both mentioned damage occurrences were perfectly predicted by numerical simulation. As shown in Figure 8, tensile fiber failure (SDV1), compressive fiber failure (SDV2), tensile matrix failure (SDV3), and matrix crushing failure (SDV4) appeared perfectly at the corresponding zones. Moreover, the region of bearing damage could be identified by the uplift of the material around the holes and it extended about half of the hole diameter from the wall of the hole (Figure 7, light blue curves). Similar results were obtained as given in Figures 8(b) and 8(d). Finally, delamination caused by secondary bending of the singlelap joint occurred around the hole and close to the loading side (green curves in Figure 7), which was also successfully predicted by the proposed model (Figure 8(f)).
(a) Tensile fiber failure (SDV1)
(b) Compressive fiber failure (SDV2)
(c) Tensile matrix failure (SDV3)
(d) Compressive matrix failure (SDV4)
(e) Delamination onset in tension (SDV5)
(f) Delamination onset in compression (SDV6)
(g) Fibermatrix shearout failure (SDV7)
Bearing failure of composite laminate is a combination of multiple damage occurrences. When subjected to tensile loading, the bolt shank is pressed against the bolt hole. Then, the fibers firstly bear the load due to their larger elastic modulus and are gradually buckling as the load increases until fracture. After that, matrix compression failure and fibermatrix shearout failure occur along with delamination between layers. Gradually, as fiber and matrix failures continue to expand, the loadcarrying capacity of the composite laminates decreases.
4.2. Parameters Study
In this section, how the loadcarrying capacity of a singlelap doublebolt joint is affected will be discussed. Numerical results are shown in Figure 9 and details are given in Table 6. Generally speaking, both stiffness and ultimate load of the joints have a tendency to decrease with increasing of tilting angle. However, the decline of stiffness is larger than that of ultimate load, which indicates that hole perpendicularity error has a great influence on joint stiffness.

(a) P1
(b) P2
(c) P3
(d) P4
(e) P5
(f) Comparison of different joint patterns
As shown in Figure 9(a), the linear elastic segments of load–displacement curves are substantially coincident for P1. The deviations of stiffness and ultimate load are less than 3% when β is no more than 3 degrees (see Table 6). The variations of all curves, from linear elastic segment to the nonlinear segment, are also consistent, which indicates that hole perpendicularity error has a minimal impact on mechanical performance of this kind of joints, while, for P2 and P3, there is a significant decline with increasing of β. When , the stiffness decreases by nearly 12%. The largest decline of joint stiffness happens in P4, which is about 23.8% even when . It is easy to understand that P4 experiences twice the influence of hole perpendicularity error of P2.
Figure 9(f) shows the comparison of different joint patterns with . The only difference between P1 and P2, as well as P4 and P5, is tilting direction α. However, the stiffness of P1 and P5 is higher than that of P2 and P4, respectively. It can be concluded that when hole tilts against the load direction , hole perpendicularity error has less influence on joint stiffness than that of . Detailed comparison results are shown in Figure 10, from which similar conclusions can be drawn.
4.3. Damage Comparison
The appearance of hole perpendicularity error will alter the load transfer path and affect loads transferred by laminates in a singlelap doublebolt joint. The expectant result is obviously that more elements are damaged around the bolt hole when it bears more load. Take tensile load of 30 kN as an example and numerical results of matrix crushing damage are given in Figure 11. Figure 11(a) is the damage nephogram for P0, which has almost the same scale damage around both holes for each laminate. But, when hole perpendicularity error is involved, element damage becomes severer and uneven, which can be seen by comparing Figures 11(b) and 11(c) with 11(a).
To prove the above conclusion, specimens loaded up to 30 kN, about 80% of the ultimate load, were disassembled to observe the damage around the holes. Pictures captured by optical microscope are shown in Figure 12, where first row is the results of P0, second row is the results of P1 with , the last row is the results of P2 with , A1 indicates the left hole of top laminate, and B2 represents the hole of bottom laminate for each configuration. Obviously, the damage of P2 is severer than that of P1 and P0. The damage of P0 and P1 has only minor differences, which proves that the parameter has less influence on the joint than that of . Taking the damage location and degree into consideration, numerical results show a good agreement with the experimental ones.
5. Conclusions
A progressive damage model based on ChangLessard criterion was implemented using UMAT subroutine of ABAQUS/standard to predict damage initiation, evolution, and ultimate failure of singlelap doublebolt composite joint. Seven modes of failures are captured to forecast the internal damage of the composite laminates. The prediction of failure for composite bolt joints shows good agreement with experimental results.
The bolt hole tilting angle up to 4° has little impact on ultimate bearing load of singlelap doublebolt composite joints. However, the stiffness is apparently affected. The stiffness of P2 and P3 decreases rapidly with increase of tilting angle.
By loading up to 30 kN, the damage of P2 is severer than that of P1, which proves that the initial damage of P2 will occur under relatively less amount of load. When composite bolt joints are subjected to tensile load, hole tilting against loading direction should be more concerned and avoided.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (Grant nos. 51375068 and 51475073), the Major State Basic Research Development Program (Grant no. 2014CB046504), the Fundamental Research Funds for the Central Universities (Grant no. DUT17JC19), and the Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry.
References
 R. Hazimeh, R. Othman, K. Khalil, and G. Challita, “Influence of plies' orientations on the stress distribution in adhesively bonded laminate composite joints subjected to impact loadings,” Composite Structures, vol. 152, pp. 654–664, 2016. View at: Publisher Site  Google Scholar
 T. E. A. Ribeiro, R. D. S. G. Campilho, L. F. M. da Silva, and L. Goglio, “Damage analysis of compositealuminium adhesivelybonded singlelap joints,” Composite Structures, vol. 136, pp. 25–33, 2016. View at: Publisher Site  Google Scholar
 S. Budhe, M. D. Banea, S. de Barros, and L. F. M. da Silva, “An updated review of adhesively bonded joints in composite materials,” International Journal of Adhesion and Adhesives, vol. 72, pp. 30–42, 2017. View at: Publisher Site  Google Scholar
 A. Valenza, V. Fiore, C. Borsellino, L. Calabrese, and G. Di Bella, “Failure map of composite laminate mechanical joint,” Journal of Composite Materials, vol. 41, no. 8, pp. 951–964, 2007. View at: Publisher Site  Google Scholar
 C. Cooper and G. J. Turvey, “Effects of joint geometry and bolt torque on the structural performance of single bolt tension joints in pultruded GRP sheet material,” Composite Structures, vol. 32, no. 14, pp. 217–226, 1995. View at: Publisher Site  Google Scholar
 P. J. Gray, R. M. O'Higgins, and C. T. McCarthy, “Effect of thickness and laminate taper on the stiffness, strength and secondary bending of singlelap, singlebolt countersunk composite joints,” Composite Structures, vol. 107, no. 1, pp. 315–324, 2014. View at: Publisher Site  Google Scholar
 A. Riccio and L. Marciano, “Effects of geometrical and material features on damage onset and propagation in singlelap bolted composite joints under tensile load: Part I  Experimental studies,” Journal of Composite Materials, vol. 39, no. 23, p. 2071, 2005. View at: Publisher Site  Google Scholar
 G. Kolks and K. I. Tserpes, “Efficient progressive damage modeling of hybrid composite/titanium bolted joints,” Composites Part A: Applied Science and Manufacturing, vol. 56, pp. 51–63, 2014. View at: Publisher Site  Google Scholar
 S. Heimbs, S. Schmeer, J. Blaurock, and S. Steeger, “Static and dynamic failure behaviour of bolted joints in carbon fibre composites,” Composites Part A: Applied Science and Manufacturing, vol. 47, no. 1, pp. 91–101, 2013. View at: Publisher Site  Google Scholar
 C. T. McCarthy, M. A. McCarthy, and V. P. Lawlor, “Progressive damage analysis of multibolt composite joints with variable bolthole clearances,” Composites Part B: Engineering, vol. 36, no. 4, pp. 290–305, 2005. View at: Publisher Site  Google Scholar
 C. T. McCarthy and M. A. McCarthy, “Threedimensional finite element analysis of singlebolt, singlelap composite bolted joints. Part II. Effects of bolthole clearance,” Composite Structures, vol. 71, no. 2, pp. 159–175, 2005. View at: Publisher Site  Google Scholar
 H. Gao, J. Wang, Y. Yang et al., “Effect of hole perpendicularity error on mechanical behavior of singlelap composite joints,” Astronautica Sinica, vol. 02, pp. 285–293, 2017. View at: Google Scholar
 ASTM standard D 5961/D 5961M–96, Standard Test Method for Bearing Response of Polymer Matrix Composite Laminates, 1996.
 AbaqusInc., Abaqus analysis user’s guide, Version 6.14, 2014, http://50.16.225.63/v6.14/books/usb/default.htm?startat=pt06ch28.html.
 F.K. Chang and K.Y. Chang, “A Progressive Damage Model for Laminated Composites Containing Stress Concentrations,” Journal of Composite Materials, vol. 21, no. 9, pp. 834–855, 1987. View at: Publisher Site  Google Scholar
 S. W. Tsai and E. M. Wu, “A general theory of strength for anisotropic materials,” Journal of Composite Materials, vol. 5, pp. 58–80, 1971. View at: Publisher Site  Google Scholar
 Z. Hashin, “Failure criteria for unidirectional fiber composites,” Journal of Applied Mechanics, vol. 47, no. 2, pp. 329–334, 1980. View at: Publisher Site  Google Scholar
 A. Puck and H. Schürmann, “Failure analysis of FRP laminates by means of physically based phenomenological models,” Composites Science and Technology, vol. 58, no. 7, pp. 1045–1067, 1998. View at: Publisher Site  Google Scholar
 F.K. Chang and L. B. Lessard, “Damage tolerance of laminated composites containing an open hole and subjected to compressive loadings. Part I. Analysis,” Journal of Composite Materials, vol. 25, no. 1, pp. 2–43, 1991. View at: Publisher Site  Google Scholar
 Á. Olmedo and C. Santiuste, “On the prediction of bolted singlelap composite joints,” Composite Structures, vol. 94, no. 6, pp. 2110–2117, 2012. View at: Publisher Site  Google Scholar
 S. C. Tan, “Progressive failure model for composite laminates containing openings,” Journal of Composite Materials, vol. 25, no. 5, pp. 556–577, 1991. View at: Publisher Site  Google Scholar
 F.K. Chang and K.Y. Chang, “PostFailure Analysis of Bolted Composite Joints in Tension or ShearOut Mode Failure,” Journal of Composite Materials, vol. 21, no. 9, pp. 809–833, 1987. View at: Publisher Site  Google Scholar
 P. P. Camanho and F. L. Matthews, “A progressive damage model for mechanically fastened joints in composite laminates,” Journal of Composite Materials, vol. 33, no. 24, pp. 2248–2280, 1999. View at: Publisher Site  Google Scholar
 C. T. McCarthy, R. M. O'Higgins, and R. M. Frizzell, “A cubic spline implementation of nonlinear shear behaviour in threedimensional progressive damage models for composite laminates,” Composite Structures, vol. 92, no. 1, pp. 173–181, 2010. View at: Publisher Site  Google Scholar
 K. I. Tserpes, G. Labeas, P. Papanikos, and T. Kermanidis, “Strength prediction of bolted joints in graphite/epoxy composite laminates,” Composites Part B: Engineering, vol. 33, no. 7, pp. 521–529, 2002. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Xueshu Liu 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.