Skip to content
Jad El Badaoui

Composites · FEA · Experiment · 2025/26 · Final-year research project, University of Bristol

Thin-Ply Composite Analysis

Does the Double-Double strength penalty survive when the plies get thin? Theory, manufacture, experiment and FE, combined to answer one question.

Repository on GitHub
  • Abaqus
  • MATLAB
  • Python
  • eLamX
  • Instron 1342
  • SEM

Undergraduate research project (AENG30017), supervised by Prof. Michael R. Wisnom

Three stages of composite manufacture: a laid-up laminate panel, the vacuum-bagged assembly, and the autoclave used for curing
The panels behind every number on this page: laid-up laminate, vacuum-bagged assembly, and the autoclave curing system.
4 × 32-ply
Laminates manufactured

TC33/K51 spread-tow prepreg, 0.03 mm plies, autoclave-cured

496 to 538 MPa
Failure stress range

All four architectures, failure strains 1.41 to 1.53 %

45 to 254 J/m²
O'Brien energy release rates

Within or below the Gc ≈ 200 to 500 J/m² critical range

65
Automated tests

Asserting the analysis pipeline's core properties

01 / Overview

The research question

Conventional-thickness Double-Double (DD) laminates carry a well-documented unnotched tensile penalty relative to stiffness-equivalent quasi-isotropic (QI) laminates, usually attributed to elevated free-edge stresses, earlier matrix cracking and delamination.

This project asks whether that penalty is intrinsic to the DD architecture, or whether it is only a consequence of matrix- and interface-dominated damage mechanisms that thin plies suppress. It asks the same question of the classical strength penalty associated with ply blocking.

The four architectures. Blocking groups four plies of the same angle together while preserving the in-plane stiffness of the dispersed equivalent, so any strength difference is a stacking effect rather than a stiffness effect.

Each laminate is 32 plies at 0.03 mm, giving a 0.96 mm laminate. The blocked stacks contain the same proportion of each orientation as their dispersed counterparts, so their extensional stiffness (A) matrices are identical. That equivalence is the premise the whole comparison rests on, and it is asserted by an automated test rather than assumed.

  1. 01

    Theory

    CLT, ABD matrices, Hashin first-ply failure, O'Brien energy release rates

  2. 02

    Manufacture

    Hand layup, vacuum bagging and autoclave cure of four 32-ply panels

  3. 03

    Experiment

    Unnotched tension on an Instron 1342 with video-gauge strain

  4. 04

    FE model

    3D ply-by-ply Abaqus free-edge stress extraction

  5. 05

    Statistics

    Welch tests with Holm correction for multiple comparisons

  6. 06

    Fractography

    SEM identification of the governing damage mechanisms

02 / Process

Manufacturing the panels

Panels were laid up by hand from TC33/K51 spread-tow thin-ply prepreg, vacuum-bagged and autoclave-cured, then cut into untabbed coupons.

Hand layup of the spread-tow thin-ply prepreg.
Panel bagged for debulking and cure.
Autoclave loaded, before cure.
Cured panel, with surface wrinkling visible across the plate.

Working at 0.03 mm ply thickness makes layup defects a live concern, and the repository documents them rather than editing them out. They are the reason results are ultimately interpreted comparatively rather than as intrinsic material strengths.

The defect record below is part of the engineering argument. It covers wrinkling after debulking, wrinkling concentrated on the positive-angle plies, fibre separation in the spread-tow material, and a ply repair where a second ply was laid over a broken one.

Defect record

Wrinkling after debulking.
Wrinkling on positive-angle plies.
Fibre separation in the spread-tow material.
Ply repair: a second ply laid over a broken one.

03 / Experiment

Tensile testing

A cut coupon, ready for test.
The Instron 1342 rig, with hydraulic grips and the video gauge on its tripod.
Coupon mounted for test.

The non-contact video gauge matters on a 0.96 mm coupon. A bonded extensometer or strain gauge would introduce exactly the local stiffening and stress concentration the experiment is trying to avoid. Coupons were run untabbed after preliminary tabbed tests failed prematurely at the tab ends.

04 / Numerical

The 3D free-edge finite-element model

The FE work is the numerical half of the argument. If DD really carries an intrinsic penalty, it should show up as a more severe free-edge stress state than QI at matched in-plane stiffness.

Every one of the 32 plies is meshed discretely through the thickness. The free-edge problem is a ply-scale effect, so smearing the laminate into an equivalent shell would delete the very thing being measured.

Model & boundary conditions

One end fully fixed, the other pulled to Ux = 1 mm. Interlaminar stresses are read along a path across the width at mid-length, far from the grip constraints.

Through-thickness mesh

All 32 plies individually meshed through the thickness.
The three interlaminar components extracted at the central same-angle interface: S33 (peel), S13 and S23 (interlaminar shear). The interface of interest is the 90°/90° pair in QI and the 67.5°/67.5° pair in DD.

All four laminates

|S33| across the coupon width for all four laminates. The interior is essentially stress-free, and everything happens in the last 1 mm or so at each edge.

QI blocked, signed

The same picture in signed form for QI blocked: flat through the interior, then a steep excursion to roughly −175 MPa right at the free edge.

Two results come out of this. QI and DD develop genuinely different free-edge stress states despite matched in-plane stiffness, so the architectures are not interchangeable at the ply scale even when their ABD matrices say they are. Blocking also raises the local magnitudes in both families, which is the conventional explanation for the blocking strength penalty.

Feeding these fields into the O'Brien energy release rate gives per-interface values of 45 to 254 J/m², sitting within or below the expected critical range for this material system (Gc ≈ 200 to 500 J/m²) for three of the four laminates. The free-edge stresses are real and architecture-dependent, but mostly not large enough to drive delamination first. That is the crux of the whole result.

05 / Results

What the experiment showed

All four laminates failed at similar global strains (1.41 to 1.53 %) and stresses (496 to 538 MPa). Blocked DD gave the highest mean, but the question worth answering is which of the differences are real.

Mean stress-strain response. The two architectures track each other almost exactly to failure, with DD sitting marginally higher. Shaded bands are the scatter across specimens.

Local ply strain at failure

Critical local ply strain at failure. QI reaches 1.42 % against a 1.5 % reference, DD only 1.16 %.

Failure stress vs reference

Failure stress: QI 13 % below its reference value, DD 25 % below.

Classical laminate theory converts the laminate-level measurements into something comparable against the material's own fibre allowable. Both architectures underperform their references and DD underperforms by more, which taken alone looks like the classic DD penalty. It is not. The shortfall tracks the grip-failure caveat below, and the architectures are statistically indistinguishable when compared against each other rather than against literature values.

Raw p-values

Raw pairwise Welch p-values for failure strain. Several comparisons appear significant.

After Holm correction

The same comparisons after Holm correction. Only one pairing still clears 0.05.

That pair is what decides the question. Raw pairwise tests throw up differences, but once Holm correction is applied for the number of comparisons actually being made, almost all of them disappear. The headline QI-versus-DD difference does not survive it.

06 / Mechanism

Fractography closes the loop

SEM of the fracture surfaces makes the mechanism visible directly. QI shows matrix-dominated separation localised at the 90° ply, so a single weak orientation is doing the failing. DD shows fibre fracture across all four orientations, with damage distributed rather than concentrated at one ply.

That is what the thin-ply suppression hypothesis predicts. With no weak matrix-dominated ply to fail first, the DD architecture loses the mechanism that conventionally penalises it. The statistics say the penalty is gone, and the fractography says why.

QI (top) and DD (bottom), with full fracture surfaces on the left and magnified cross-sections across the ply orientations on the right. Matrix cracking is annotated on the DD surface.

07 / Methods

Reproducibility as an engineering deliverable

The reported numbers come from the MATLAB pipeline, which covers ABD construction, Hashin first-ply failure across all four modes, O'Brien energy release rates, and the Welch and Holm statistics. Because MATLAB needs a licence, the central calculation is also implemented in Python and runs against a clearly labelled synthetic sample dataset with no licensed software at all.

A 65-test suite asserts the properties the analysis depends on rather than just that the code runs. It checks that a symmetric laminate really produces a zero B matrix, that the A matrix is genuinely unchanged by blocking (the stiffness-matching premise of the whole comparison), that a unidirectional laminate recovers its engineering constants exactly, that Hashin returns unity at each allowable, and that the Holm implementation matches hand-computed adjusted p-values.

The complete write-up lives in the repository: methods, data, code, provenance and limitations.

github.com/Jadbadawi/thin-ply-composite-analysis