The concept of lightweight design is driving future aircraft to exploit all the available strength of materials and further reduce the weight of an aircraft leading to lower fuel consumption and more sustainable aviation. Current designs introduce rivets and bolts to join the structure and simultaneously create holes in the pristine material, which is not wanted due to the local stress concentration. If the design wants to exploit all the strength of the material, stress concentration should be avoided using suitable joining technologies that don't require holes in the structure. This is possible with the adhesive bonding technology. However, many factors, such as the effects of manufacturing and impact-induced damage are still not fully understood. Thus, the damage tolerance of the design cannot be guaranteed. This results in conservatory safety factors being prescribed in the design process, which possibly reduces the exploitation of all available material strength.

The occurrence of barely visible impact damage (BVID) in aircraft composite structures, especially in adhesive joints, is a serious issue that can jeopardise an aircraft's structural safety during operation. However, virtual testing and numerous numerical approaches allow high-fidelity simulation of damage initiation and propagation in adhesives and composite adherends to predict the residual strength of the bonding accurately. Although the development is fast, accurate impact simulation is still computationally very expensive and thus not applicable in industrial cases where behaviour needs to be analysed to assess the structures' design, maintenance and repair.

This thesis investigated two topics that allow robust and accurate simulation methods to assess the effect of manufacturing and impact damage on the residual strength of bonded joints. The first topic is related to the development of the composite material damage model, where the multiscale material model is created by the use of a representative volume element (RVE). The second topic of simulation methodology is the development of the quasi-static simulation approach for the damage tolerance assessment of bonded joints. Modelling approaches to simulate residual strength are investigated through the finite element modelling of three groups (i) pristine, (ii) artificially damaged during manufacturing, and (iii) impacted bonded joints. In the numerical models, an approach based on the observation of the fracture surface of single-lap joints in different geometry and layup configuration is proposed in which the damage assessment focuses on the bondline area up to the first 0° ply. For the modelling of the damage in joints, two modelling techniques were studied, first removing elements in the damaged area and second detachment of the elements. Utilising those techniques, simplified approaches to model damage resulting from the impact were studied, with the modelling of impact damage as a hole, where all through-thickness elements are deleted and as delamination, where interlaminar cohesive zone elements are deleted.

The developed multiscale material model allowed for an accurate representation of failure modes occurring in the composite material by the characterisation of elastic, damage and plasticity parameters of the fibre and matrix constituents in the homogenisation and inverse characterisation process. The results showed that the deletion of the elements can be used to represent defects and damage of different kinds in the composite bonded joints, both in the adhesive and adherend. The comparison of different representations of the impact damage in the single lap joint configuration revealed that the best prediction in terms of the ultimate load, failure mode and size of the model is obtained with the simplified representation as a single delamination positioned before the first 0° ply in the layup and in the studied adherend configuration that was between the first (45°) and second (0°) ply in the layup. The study ends with the conclusion that in different geometries of single lap joints and different damage types studied, a numerical analysis should focus on the region of the overlap edge and in the thickness direction from the bondline up to the first 0° ply in the lay-up. The results and numerical method developed during this thesis create a base for the further investigation of a variety of impact cases on bonded joints.