International Journal of Polymer Science
Volume 2017 (2017), Article ID 2790198, 11 pages
https://doi.org/10.1155/2017/2790198

Research Article

Effects of Hole Perpendicularity Error on Mechanical Performance of Single-Lap Double-Bolt Composite Joints

Xueshu Liu,¹ Yuxing Yang,² Yiqi Wang,² Yongjie Bao,² and Hang Gao²

¹School of Automotive Engineering, Dalian University of Technology, Dalian 116024, China
²School of Mechanical Engineering, Dalian University of Technology, Dalian 116024, China

Correspondence should be addressed to Xueshu Liu; liuxs@dlut.edu.cn

Received 9 August 2017; Revised 10 October 2017; Accepted 7 November 2017; Published 27 November 2017

Academic Editor: Peng He

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.

Abstract.

Static tensile experiments and numerical simulations were carried out to study effects of hole perpendicularity error on mechanical performance of single-lap double-bolt 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 Chang-Lessard 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 load-carrying capacity is concerned, bolted joints are the most important methods of composite assembling especially in aerospace industry because of their high load-carrying 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 load-carrying 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 physical-mechanical 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 single-lap single-bolt 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 single-lap double-bolt 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)

Figure 1: Characterization of perpendicularity error of hole from [9].

Only tilting directions and are discussed to reduce the parameters of the study. The influence of tilting direction for double-bolt 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

Figure 2: Patterns of single-lap double-bolt joints.

3. Methodology

3.1. Specimens and Experimental Procedures

The specimens used in this study were manufactured from IMS194/CYCOM977-2 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.

Table 1: Mechanical properties of IMS194/CYCOM977-2.

Elastic property
Value
Strength property
Value

165

3150

1450

4.27

2.52

0.33

108

108

0.55

108

value is set according to the transversely isotropic assumption and engineering experience.

Table 2: Plastic behavior of SCM35 alloy steel.

Stress/GPa
Plastic strain

0.63
0

0.76
0.01

1.08
0.05

1.13
0.08

1.2
0.1

Single-lap double-bolt composite joints were prepared according to the ASTM standard D5961/D5661 M-13 with width-to-hole diameter ratio and edge distance-to-hole diameter ratio [13]. Standard specimen configurations using fixed values of test parameters are shown in Figure 3.

Figure 3: Configuration of specimens (unit: mm).

According to the ASTM D5961 standard, the static tensile tests were performed on a WDW-100 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 YSJ50-5-ZC 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 single-lap double-bolt composite joint.

Table 3: Test groups and parameters.

Group
Pattern
β/°
Repeat

P0-0
P0
0
3

P0-0-30
P0
0
1

P1-4
P1
4
3

P1-4-30
P1
4
1

P2-2
P2
2
3

P2-4
P2
4
3

P2-4-30
P2
4
1

Figure 4: Experimental setup.

3.2. Numerical Model

To investigate the effects of hole perpendicularity error on mechanical behavior of single-lap double-bolt composite joints, a nonlinear three-dimensional finite element model was constructed with 8-node solid element C3D8R using the commercial finite element software, ABAQUS [14], as illustrated in Figure 5.

Figure 5: FEA model and interactions of bolted joint.

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 width-to-thickness ratio of the element [7].

Contact interactions between all surfaces were modelled by surface-based 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 two-dimensional model into three-dimensional model. There are diverse failure criteria to predict the damage initiation including Chang criteria [15], Tsai-Wu criteria [16], Hashin type criteria [17], and Puck criteria [18]. In this paper, a damage model considering the effect of nonlinear shear stress-strain behavior was established based on extended Chang-Lessard 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, :

Fiber-matrix shear-out 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⁻⁸.

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.

Table 4: Material property degradation for each failure criterion.

Model of failure
Material properties degradation rules

Matrix cranking failure
  ,

Matrix crushing failure

Tensile fiber failure

Compressive fiber failure

Fiber-matrix shear-out failure

Delamination onset in tension

Delamination onset in compression

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.

Figure 6: Comparison of numerical and experimental load–displacement curves.

As shown in Table 5, the average ultimate loads of P0, P1-4, and P2-4 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 single-lap double-bolt 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.

Table 5: The experimental and numerical bearing behaviors of all specimens.

Group
Stiffness 
(N/mm)
Mean 
(N/mm)
Ultimate load 
(N)
Mean 
(N)

P0-0
33289
30711
38206
38047

P0-0
30424
38568

P0-0
28421
37367

P1-4
31889
31181
37800
37734

P1-4
30964
37668

P1-4
30691
37976

P2-2
26500
28083
38894
37450

P2-2
28673
37783

P2-2
29076
35671

P2-4
27171
26540
37964
37745

P2-4
25804
37526

P2-4
26645
37535

FEA-P0
—
36690
—
35900

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

Figure 7: Failure morphology of the top laminate under failure load.

(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) Fiber-matrix shear-out failure (SDV7)

Figure 8: Damage envelope of various failure modes under ultimate load.

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 fiber-matrix shear-out failure occur along with delamination between layers. Gradually, as fiber and matrix failures continue to expand, the load-carrying capacity of the composite laminates decreases.

4.2. Parameters Study

In this section, how the load-carrying capacity of a single-lap double-bolt 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.

Table 6: The stiffness and ultimate load of FEA results.

Group
Stiffness/N/mm
Decline/%
Ultimate load/N
Decline/%

P0
26515
—
35901
—

P1-1
26325
0.72
35552
0.97

P1-2
26061
1.71
35336
1.58

P1-3
25751
2.88
34901
2.79

P1-4
25187
5.01
33394
6.98

P2-1
25154
5.13
35077
2.30

P2-2
23389
11.79
35052
2.36

P2-3
21805
17.76
34419
4.13

P2-4
21482
18.98
33302
7.24

P3-1
25327
4.48
35506
1.10

P3-2
23453
11.55
35353
1.53

P3-3
22015
16.97
34088
4.28

P3-4
22005
17.01
33182
7.57

P4-1
25090
5.4
35558
0.96

P4-2
20203
23.8
35429
1.31

P4-3
16926
36.2
34675
3.41

P4-4
15501
41.5
33441
6.85

P5-1
26161
1.3
35388
1.43

P5-2
23667
10.7
35453
1.25

P5-3
23185
12.5
35013
2.47

P5-4
21140
20.3
32830
8.55

(a) P1

(b) P2

(c) P3

(d) P4

(e) P5

(f) Comparison of different joint patterns

Figure 9: Numerical results of hole perpendicularity error on joint performance.

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.

Figure 10: The influence of hole tilting direction on joint stiffness.

4.3. Damage Comparison

The appearance of hole perpendicularity error will alter the load transfer path and affect loads transferred by laminates in a single-lap double-bolt 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).

Figure 11: Load transfer path and damage of matrix crushing failure.

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 , A-1 indicates the left hole of top laminate, and B-2 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.

Figure 12: Hole damage pictures for specimen loaded up to 30 kN.

5. Conclusions

A progressive damage model based on Chang-Lessard criterion was implemented using UMAT subroutine of ABAQUS/standard to predict damage initiation, evolution, and ultimate failure of single-lap double-bolt 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 single-lap double-bolt 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.

T. E. A. Ribeiro, R. D. S. G. Campilho, L. F. M. da Silva, and L. Goglio, “Damage analysis of composite-aluminium adhesively-bonded single-lap joints,” Composite Structures, vol. 136, pp. 25–33, 2016.

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.

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.

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. 1-4, pp. 217–226, 1995.

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 single-lap, single-bolt countersunk composite joints,” Composite Structures, vol. 107, no. 1, pp. 315–324, 2014.

A. Riccio and L. Marciano, “Effects of geometrical and material features on damage onset and propagation in single-lap bolted composite joints under tensile load: Part I - Experimental studies,” Journal of Composite Materials, vol. 39, no. 23, p. 2071, 2005.

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.

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.

C. T. McCarthy, M. A. McCarthy, and V. P. Lawlor, “Progressive damage analysis of multi-bolt composite joints with variable bolt-hole clearances,” Composites Part B: Engineering, vol. 36, no. 4, pp. 290–305, 2005.

C. T. McCarthy and M. A. McCarthy, “Three-dimensional finite element analysis of single-bolt, single-lap composite bolted joints. Part II. Effects of bolt-hole clearance,” Composite Structures, vol. 71, no. 2, pp. 159–175, 2005.

H. Gao, J. Wang, Y. Yang et al., “Effect of hole perpendicularity error on mechanical behavior of single-lap composite joints,” Astronautica Sinica, vol. 02, pp. 285–293, 2017.

ASTM standard D 5961/D 5961M–96, Standard Test Method for Bearing Response of Polymer Matrix Composite Laminates, 1996.

Abaqus-Inc., 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.

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.

Z. Hashin, “Failure criteria for unidirectional fiber composites,” Journal of Applied Mechanics, vol. 47, no. 2, pp. 329–334, 1980.

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.

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.

Á. Olmedo and C. Santiuste, “On the prediction of bolted single-lap composite joints,” Composite Structures, vol. 94, no. 6, pp. 2110–2117, 2012.

S. C. Tan, “Progressive failure model for composite laminates containing openings,” Journal of Composite Materials, vol. 25, no. 5, pp. 556–577, 1991.

F.-K. Chang and K.-Y. Chang, “Post-Failure Analysis of Bolted Composite Joints in Tension or Shear-Out Mode Failure,” Journal of Composite Materials, vol. 21, no. 9, pp. 809–833, 1987.

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.

C. T. McCarthy, R. M. O'Higgins, and R. M. Frizzell, “A cubic spline implementation of non-linear shear behaviour in three-dimensional progressive damage models for composite laminates,” Composite Structures, vol. 92, no. 1, pp. 173–181, 2010.

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.