Preliminary Experimental and Numerical Study of Metal Element with Notches Reinforced by Composite Materials

: The presented study is related to the application of the composite overlays used in order to decrease the effect of the stress concentrations around the cut-outs in structural metal elements. The proposed approach with the application of the digital image correlation extends the recently presented studies. Such structural elements with openings of various shapes have been accommodated for a wide range of industrial applications. These structures exhibit certain stress concentrations which decrease their durability and strength. To restore their strength, various reinforcing overlays can be used. In the present paper, the ﬂat panel structure without and with the composite overlays made of HEXCEL TVR 380 M12/26%/R-glass/epoxy is under the experimental and the numerical study. Particular attention is paid to the investigation of the samples with the rectangular holes, which for smooth rounded corners offer a higher durability than the samples with the circular hole of the same size. The experimental results are obtained for the bare element and are reinforced with composite overlay samples. The experimental results are obtained with the use of the Digital Image Correlation method, while the numerical results are the product of the Finite Element Analysis. In the numerical analysis, the study of the shape, size and ﬁber orientation in applied overlays is done. The reduction of the stress concentration observed in opening notches has conﬁrmed the effectiveness of the overlay application. In the investigated example, the application of the square composite overlay increased the structure strength even by 25%.


Introduction
The structural elements with voids or openings of different shapes are commonly used in design practice. The main reason for this is that it is necessity to provide access to certain parts of the structure or the possibility to supply power or other media to perform the production process or correct operation. Such intentionally made holes are usually the source of the stress concentration, which may substantially decrease the structure durability or the fatigue endurance in cases of repetitive loading imposition. Another possible danger is the appearance of zones with plastic deformations, which are usually not allowed in standard structure operations. Such zones around the different holes have been recently analyzed using the Digital Image Correlation (DIC) [1], which is under rapid development and is successfully used for damage monitoring of different materials, such as isotropic metals [1], cementitious composites [2], and composites [3]. In the presented paper, special attention is focused on the reduction of such stress concentrations in notches reached through the application of the composite overlays around the cut-out hole.
The problem of stress concentration in constructional elements has attracted researchers for more than one hundred years and is still under thorough study due to the application of more and more advanced materials and the increasing shape complexity of the structural elements [1,. The first, basic mathematical formulas for the plane a single stiffened joint (from theoretical and numerical calculations) signalize significant stress concentrations at the ends of the adhesive joints. More information about such theoretical formulations can be found in Refs. [41,44,45]. There are also methods for analyzing materials at the macro-, meso-, and micro-levels [48,49]. Despite the many theoretical models as well as techniques (fracture mechanics [38,50], cohesive zone model [40,45], FEM [39,45,51,52], continuum mechanics), there is no universal approach for predicting the strength of adhesive joints [44].
Based on the Volkersen solution [46], the shear stress in the adhesive layer varies along the overlap length. Such shear stresses depend on the strain distribution in the adherents. Due to the differences in deformations, the highest shear stresses in the adhesive occur at both ends and are minimal at the center of the adhesive joint. It should be noted that stress concentrations of the shear stresses in the adhesive layer may achieve high values (typically 1-4, but in some cases even above 10) [40,41,45,51]. Because of this, the failure typically initiates at both ends of the adhesive joint [40,52]. Due to the high-stress concentrations at the ends of the adhesive joint and degradation of the adhesive layer over time, such connections require periodic inspections in order to control and repair damages. Such monitoring can be made with the use of the traditional non-destructive techniques [45,53], vibration-based methods [35], digital image correlation [1,3] or c-scan inspection [54].
Except for the problems with the adhesive joint choice, several advantages are present in a metal adhesive layer composite material combination. These are attributed to the nature of the composite material. The number of composite layers, choice of fiber material and matrix resin type, orientation of fibers in stacking layers and the size of the overlay can be freely shaped during the study of the overlays, and the above parameters can be the subject of the optimal choice for the analyzed problem.
The main aim of the presented paper is the experimental and numerical study of possibilities of reduction of the effect of the stress concentrations on the structure strength under tensile load. The proposed approach, with the application of the digital image correlation, extends the recently presented studies. The analyses are performed for metal samples with rectangular holes reinforced by the composite overlays. The reduction of the stress concentration observed in opening notches has confirmed the usefulness of the overlay application. In the investigated example, the application of the square composite overlay increased the structure strength even by 25%. The paper consists of five sections. The introduction to the studied problem and the literature review is given in Section 1. The description of the applied materials, samples and methods are described in Section 2. The results of the experimental investigations are presented and discussed in Section 3. The results of the numerical analyses are presented in Section 4. In these numerical studies, the influence of the geometry of the overlays and orientations of the layers on the structural strength is examined. The summary of the performed study and the final discussion are given in Section 5. The conclusions are given in Section 6.

Materials and Samples
The presented numerical and experimental studies are made for the samples made of structural steel with different cut-outs reinforced by composite overlays. The structural low-alloy steel S355J2+N was used as the base material. Such steel can be characterized by high strength and ductility, as well as good weldability, and is generally used in the manufacturing of machine parts, pressure vessels and pipes, and in civil engineering. The chemical composition of the material used in the experimental tests is given in Table 1. The maximal values of the particular chemical elements, as well as the minimal mechanical properties according to Standards [55], are also reported in Table 1. The reinforcement overlays were made of TVR 380 M12/26%/R-glass/epoxy composite material. The plate used in the preparation of the reinforcement overlays was made from prepregs with fiber and epoxy resin densities equal to 2.56 g·cm −3 and 1.24 g·cm −3 , respectively. More details about the preparation of the composite plate, including information about the autoclave operation, can be found in Ref. [56]. The mechanical properties of the composite material were determined with the use of the experimental tests for samples with orientations θ = [0 • ] 8 , θ = [90 • ] 8 , θ = [+45 • /−45 • ] 4 (θ = 0 • is the direction of the tension) according to the Standards [57]. Four samples were tested for each of the above cases with the loading rate equal to 3 mm/min. The detailed description and results of the tensile tests are presented in Ref [56], and the most important information is reported in Table 2. The obtained mechanical properties of composite material are in good agreement with the literature data [58].   8 ) rapid matrix cracking occurred. For this configuration, greater variance of the tensile strength is achieved. This may be caused by the geometry of the sample or matrix strength. However, the mean value from the tensile test is in good agreement with the literature data [58]. The average tensile strengths for investigated orientations are reported in Table 2.
The stress-strain curve of the base S355J2+N material was determined with the use of the sample with geometry given in Figure 2a. The dimensions of the sample were set according to the Standards [59]. The main experimental tests were carried for the samples with circular Ø15 mm ( Figure 2b) and rectangular 15 × 15 mm holes. In the case of rectangular cut-outs, two fillet radii were assumed-R2 ( Figure 2c) and R4 (Figure 2d). The hole was located at the center of the plate. The thickness of the base metal samples was 4 mm. The final test was made for the sample with reinforcement overlays on both sides (Figure 2e). The loading rate of tensile tests was equal 0.5 mm/min. More information about tests and samples are summarized in Table 3.
In the second case (Figure 1b-θ = [+45°/−45°]4) slow matrix cracking and separation of the fibers from the matrix were observed in the whole volume of the sample. In the last case (Figure 1c-θ = [90°]8) rapid matrix cracking occurred. For this configuration, greater variance of the tensile strength is achieved. This may be caused by the geometry of the sample or matrix strength. However, the mean value from the tensile test is in good agreement with the literature data [58]. The average tensile strengths for investigated orientations are reported in Table 2.  The stress-strain curve of the base S355J2+N material was determined with the use of the sample with geometry given in Figure 2a. The dimensions of the sample were set according to the Standards [59]. The main experimental tests were carried for the samples with circular Ø15 mm ( Figure 2b) and rectangular 15 × 15mm holes. In the case of rectangular cut-outs, two fillet radii were assumed-R2 ( Figure 2c) and R4 (Figure 2d). The hole was located at the center of the plate. The thickness of the base metal samples was 4 mm. The final test was made for the sample with reinforcement overlays on both sides (Figure 2e). The loading rate of tensile tests was equal 0.5 mm/min. More information about tests and samples are summarized in Table 3.    The overlays with a rectangular shape and dimensions 45 × 45 mm were made of the TVR 380 glass/epoxy composite material with the layer orientation θ = [+45°/−45°]4 (Figure 3). The thicknesses of the overlays were equal to about 2.12 mm. The microscopic image of the composite overlay through its thickness is also given in Figure 3c. The overlays were stuck to the metal base part with the use of S&P Resin 220 epoxy adhesive, which is dedicated to the reinforcement of steel structures by fiber composites. The compressive and shear strength of the used resin should not be less than 70 MPa and 26 MPa, respectively, and the minimum modulus of elasticity is equal to 7100 MPa. The epoxy adhesive was applied between the surfaces of the reinforcement overlays and the base metal part (red color in Figure 3a) and at both ends of composite overlays. The angle of the epoxy adhesive triangular spew fillets at the ends of the overlays was 45° for load transfer over a larger area.
The prepared samples with holes made of S355J2+N are presented in Figure 4 and the sample with composite reinforcement overlays is presented in Figure 5. The surfaces   The overlays with a rectangular shape and dimensions 45 × 45 mm were made of the TVR 380 glass/epoxy composite material with the layer orientation θ = [+45 • /−45 • ] 4 ( Figure 3). The thicknesses of the overlays were equal to about 2.12 mm. The microscopic image of the composite overlay through its thickness is also given in Figure 3c. The overlays were stuck to the metal base part with the use of S&P Resin 220 epoxy adhesive, which is dedicated to the reinforcement of steel structures by fiber composites. The compressive and shear strength of the used resin should not be less than 70 MPa and 26 MPa, respectively, and the minimum modulus of elasticity is equal to 7100 MPa. The epoxy adhesive was applied between the surfaces of the reinforcement overlays and the base metal part (red color in Figure 3a) and at both ends of composite overlays. The angle of the epoxy adhesive triangular spew fillets at the ends of the overlays was 45 • for load transfer over a larger area. terns was necessary for the application of the digital image correlations (DIC). This technique allows determining the distribution of the strain field on the tested surface. For this purpose, a series of photographs with constant time intervals (with the use of the Nikon D90 camera) were made during the tensile tests. The DIC analyses were made in the GOM Software [60]. The photographs were made in the constant time interval equal to 6 seconds. All analyses performed in GOM were made with the same configuration-the facet size and the point distance were equal to 66 and 15 pixels, respectively.    The prepared samples with holes made of S355J2+N are presented in Figure 4 and the sample with composite reinforcement overlays is presented in Figure 5. The surfaces of the samples were covered with random black speckle patterns with various shapes on a white background. This black/white coating with randomly distributed speckle patterns was necessary for the application of the digital image correlations (DIC). This technique allows determining the distribution of the strain field on the tested surface. For this purpose, a series of photographs with constant time intervals (with the use of the Nikon D90 camera) were made during the tensile tests. The DIC analyses were made in the GOM Software [60]. The photographs were made in the constant time interval equal to 6 seconds. All analyses performed in GOM were made with the same configuration-the facet size and the point distance were equal to 66 and 15 pixels, respectively. terns was necessary for the application of the digital image correlations (DIC). This technique allows determining the distribution of the strain field on the tested surface. For this purpose, a series of photographs with constant time intervals (with the use of the Nikon D90 camera) were made during the tensile tests. The DIC analyses were made in the GOM Software [60]. The photographs were made in the constant time interval equal to 6 seconds. All analyses performed in GOM were made with the same configuration-the facet size and the point distance were equal to 66 and 15 pixels, respectively.

Stress-Strain Behavior of S355J2+N
The engineering stress-strain curve for S355J2+N ( Figure 6) was designated with the use of the sample presented in Figure 2a. The determined upper Yield limit was higher than its minimum value [55] and was equal to 372 MPa (point 1 in Figure 6). After reaching the upper Yield limit, a slight reduction of the stress was observed. The S355J2+N material showed the Yield plateau (the range between points 1 and 2 in Figure 6) before strengthening. The maximal registered stress was equal to 543 MPa (p.3 in Figure 6). The elongation of the sample after fracture was equal to ε TOT = 31.5%.

Stress-Strain Behavior of S355J2+N
The engineering stress-strain curve for S355J2+N ( Figure 6) was designated with the use of the sample presented in Figure 2a. The determined upper Yield limit was higher than its minimum value [55] and was equal to 372 MPa (point 1 in Figure 6). After reaching the upper Yield limit, a slight reduction of the stress was observed. The S355J2+N material showed the Yield plateau (the range between points 1 and 2 in Figure  6) before strengthening. The maximal registered stress was equal to 543 MPa (p.3 in Figure 6). The elongation of the sample after fracture was equal to εTOT = 31.5%. The development of the strain field on the surface was studied with the use of DIC. Such analyses were carried in the range of strain from strain, corresponding to the upper Yield limit (εTOT = 0.85%- Figure 7a) to the material strengthening (εTOT = 3.39%- Figure  7f). Uniform stretching of the sample was observed up to the Yield limit ( Figure 7a). Above the Yield limit (point 1 in Figure 6), the strong concentration of the strain appeared in the weakest point and cross-section of the sample (Figure 7b). In the range of the Yield plateau (Figures 7b-d), the growth of the plastic strain in the surrounding of the weakest cross-section was observed. For the strain εTOT = 1.40% the second strain concentration appeared (Figure 7d). A further increase of tension leads to the plasticizing of the whole test part of the sample ( Figure 7e) and strengthening (Figure 7f). The determined stress-strain curve was used for material description in the FE simulations. The failure form is presented in Figure 8. The development of the strain field on the surface was studied with the use of DIC. Such analyses were carried in the range of strain from strain, corresponding to the upper Yield limit (ε TOT = 0.85%- Figure 7a) to the material strengthening (ε TOT = 3.39%- Figure 7f). Uniform stretching of the sample was observed up to the Yield limit ( Figure 7a). Above the Yield limit (point 1 in Figure 6), the strong concentration of the strain appeared in the weakest point and cross-section of the sample (Figure 7b). In the range of the Yield plateau (Figure 7b-d), the growth of the plastic strain in the surrounding of the weakest cross-section was observed. For the strain ε TOT = 1.40% the second strain concentration appeared (Figure 7d). A further increase of tension leads to the plasticizing of the whole test part of the sample ( Figure 7e) and strengthening (Figure 7f). The determined stress-strain curve was used for material description in the FE simulations. The failure form is presented in Figure 8.

Unreinforced Metal Samples with Holes
The influence of the stress concentrators in a form of circular and rectangular cut-outs in metal plates was investigated by the Authors in Ref. [1] Such holes or notches may significantly reduce the structure strength. Generally, the value of the stress concentration factor, denoted as K t , can be described by the given below formula: where: K t INF is the stress concentration factor for infinity plate, C W is the correction constant for a plate with finite width, C 0 , C 1 , C 2 , C 3 are the constant dependent on the cut-out shape, and f 1 , f 2 , f 3 are the functions dependent on the characteristic dimensions of the cut-out.
Such analyses were carried in the range of strain from strain, corresponding to the upper Yield limit (εTOT = 0.85%- Figure 7a) to the material strengthening (εTOT = 3.39%- Figure  7f). Uniform stretching of the sample was observed up to the Yield limit ( Figure 7a). Above the Yield limit (point 1 in Figure 6), the strong concentration of the strain appeared in the weakest point and cross-section of the sample (Figure 7b). In the range of the Yield plateau (Figures 7b-d), the growth of the plastic strain in the surrounding of the weakest cross-section was observed. For the strain εTOT = 1.40% the second strain concentration appeared (Figure 7d). A further increase of tension leads to the plasticizing of the whole test part of the sample ( Figure 7e) and strengthening (Figure 7f). The determined stress-strain curve was used for material description in the FE simulations. The failure form is presented in Figure 8.

Unreinforced Metal Samples with Holes
The influence of the stress concentrators in a form of circular and rectangular cut-outs in metal plates was investigated by the Authors in Ref. [1] Such holes or notches may significantly reduce the structure strength. Generally, the value of the stress concentration factor, denoted as Kt, can be described by the given below formula: where: Kt INF is the stress concentration factor for infinity plate, CW is the correction constant for a plate with finite width, C0, C1, C2, C3 are the constant dependent on the cut-out shape, and f1, f2, f3 are the functions dependent on the characteristic dimensions of the cut-out. The detailed formulae concerning the stress concentration factors for both rectangular and circular holes can be found in Ref. [1,[21][22][23][24]. The calculated theoretical values of the above constants, as well as stress concentration factors, are summarized in Table 4. Obviously, the largest value of Kt = 2.479 was obtained for the sample with the rectangular hole with the smallest fillet radius R2. However, based on the theoretical calculations, it is interesting that the circular hole is not the optimal one. It can be seen that the rectangular hole with a sufficiently large fillet radius shows smaller concentrations of stresses than the circular hole. Due to this fact, in the presented preliminary studies of the reinforcing of the metal plates with cut-outs with the use of composite overlays, the experimental studies were performed for sample 3 with the largest Kt. Table 4. Theoretical values of stress concentration factors for investigated samples. The detailed formulae concerning the stress concentration factors for both rectangular and circular holes can be found in Refs. [1,[21][22][23][24]. The calculated theoretical values of the above constants, as well as stress concentration factors, are summarized in Table 4. Obviously, the largest value of K t = 2.479 was obtained for the sample with the rectangular hole with the smallest fillet radius R2. However, based on the theoretical calculations, it is interesting that the circular hole is not the optimal one. It can be seen that the rectangular hole with a sufficiently large fillet radius shows smaller concentrations of stresses than the circular hole. Due to this fact, in the presented preliminary studies of the reinforcing of the metal plates with cut-outs with the use of composite overlays, the experimental studies were performed for sample 3 with the largest K t . Table 4. Theoretical values of stress concentration factors for investigated samples. The stress-strain characteristics for samples 2-4 with circular and rectangular holes are given in Figure 9. In the elastic and strain strengthening range till to strain ε TOT = 1.5%, all three samples show similar characteristics. A slightly greater force in this range (ε TOT = 0-1.5%) was registered for the sample with the circular hole. However, the maximal tensile forces (the ultimate tensile strength) were slightly greater in both samples with rectangular holes. The most significant difference in the behavior of these samples was the maximal elongation, which in the case of the sample with the circular hole was almost two times smaller than for samples with rectangular holes. The development of plastic strains on the surfaces of the tested samples, designated with the use of DIC, are presented in: Figure 10 for the sample with the circular hole; Figure 11 for the sample with the rectangular hole with fillet radius R2 and Figure 12 for the sample with the rectangular hole with fillet radius R4. In all the presented examples, the failure occurred in the weakest cross-section ( Figure 13). In the case of the sample with the circular hole, the crack direction was horizontal and was initiated in two points at which the highest strain concentrations were determined in DIC analysis ( Figure 10). In the case of the rectangular hole with rounded corners, the cracks were initiated in points on the vertical edge of the hole below the fillet radius. In both cases, such cracks appeared on the opposite corners of the rectangular hole. It can be seen that in the points at which cracks are initiated ( Figure 13), high stress (strain) concentrations are registered ( Figures 11 and 12).  The development of plastic strains on the surfaces of the tested samples, designated with the use of DIC, are presented in: Figure 10 for the sample with the circular hole; Figure 11 for the sample with the rectangular hole with fillet radius R2 and Figure 12 for the sample with the rectangular hole with fillet radius R4. In all the presented examples, the failure occurred in the weakest cross-section ( Figure 13). In the case of the sample with the circular hole, the crack direction was horizontal and was initiated in two points at which the highest strain concentrations were determined in DIC analysis ( Figure 10). In the case of the rectangular hole with rounded corners, the cracks were initiated in points on the vertical edge of the hole below the fillet radius. In both cases, such cracks appeared on the opposite corners of the rectangular hole. It can be seen that in the points at which cracks are initiated ( Figure 13), high stress (strain) concentrations are registered ( Figures 11 and 12). The development of plastic strains on the surfaces of the tested samples, designated with the use of DIC, are presented in: Figure 10 for the sample with the circular hole; Figure 11 for the sample with the rectangular hole with fillet radius R2 and Figure 12 for the sample with the rectangular hole with fillet radius R4. In all the presented examples, the failure occurred in the weakest cross-section ( Figure 13). In the case of the sample with the circular hole, the crack direction was horizontal and was initiated in two points at which the highest strain concentrations were determined in DIC analysis ( Figure 10). In the case of the rectangular hole with rounded corners, the cracks were initiated in points on the vertical edge of the hole below the fillet radius. In both cases, such cracks appeared on the opposite corners of the rectangular hole. It can be seen that in the points at which cracks are initiated ( Figure 13), high stress (strain) concentrations are registered ( Figures 11 and 12). (a) (b) (c) (d) Figure 11. Distribution of strains on the surface for sample no. 3 with rectangular holes with a fillet radius R2 for strain:  (a) (b) (c) (d) Figure 11. Distribution of strains on the surface for sample no. 3 with rectangular holes with a fillet radius R2 for strain: Figure 11. Distribution of strains on the surface for sample no. 3 with rectangular holes with a fillet radius R2 for strain:

Metal Samples with Holes Reinforced by Composite Materials
The main aim of the study was the analysis of the behavior of the composite reinforcement of the metal part with cut-out. The preliminary experimental test was performed for the sample with the rectangular hole with a fillet radius R2. The selection of the cut-out shape was based on the results of the experimental tests described in Section 3.2. The selected sample revealed the highest stress concentration factor with maintaining a high load capacity at the upper Yield limit of the material and ultimate tensile strength and high elongation at the fracture. The determined stress-strain curve for the metal with composite overlays (sample 5) is given in Figure 14. The stress in Figure 14 was calculated as a nominal stress in a metal plate with a cut-out with the given below formula: where: F is the measured force, W and t are the width and the thickness of the metal plate, respectively, and a is the width of the hole. Some portion of the loading is carried out by the composite overlays through the adhesive joint.
In the obtained curve ( Figure 14) a few interesting points can be observed. Till point 1, the curve is linear. The stress in point 1 corresponds to the plasticizing of the metal sample no. 3 in areas with stress concentrations. Point 1 corresponds to the strain in Figure 12. Distribution of strains on the surface for sample no. 4 with rectangular holes with a fillet radius R4 for strain:

Metal Samples with Holes Reinforced by Composite Materials
The main aim of the study was the analysis of the behavior of the composite reinforcement of the metal part with cut-out. The preliminary experimental test was performed for the sample with the rectangular hole with a fillet radius R2. The selection of the cut-out shape was based on the results of the experimental tests described in Section 3.2. The selected sample revealed the highest stress concentration factor with maintaining a high load capacity at the upper Yield limit of the material and ultimate tensile strength and high elongation at the fracture. The determined stress-strain curve for the metal with composite overlays (sample 5) is given in Figure 14. The stress in Figure 14 was calculated as a nominal stress in a metal plate with a cut-out with the given below formula: where: F is the measured force, W and t are the width and the thickness of the metal plate, respectively, and a is the width of the hole. Some portion of the loading is carried out by the composite overlays through the adhesive joint.

Metal Samples with Holes Reinforced by Composite Materials
The main aim of the study was the analysis of the behavior of the composite reinforcement of the metal part with cut-out. The preliminary experimental test was performed for the sample with the rectangular hole with a fillet radius R2. The selection of the cut-out shape was based on the results of the experimental tests described in Section 3.2. The selected sample revealed the highest stress concentration factor with maintaining a high load capacity at the upper Yield limit of the material and ultimate tensile strength and high elongation at the fracture. The determined stress-strain curve for the metal with composite overlays (sample 5) is given in Figure 14. The stress in Figure 14 was calculated as a nominal stress in a metal plate with a cut-out with the given below formula: where: F is the measured force, W and t are the width and the thickness of the metal plate, respectively, and a is the width of the hole. It can be observed that below this value the relationship between stress and strain is linear. The significant increase in the curvature of the σ-ε curve is visible for stresses above 400 MPa. This effect is clearly visible in Figure 15, where stress-strain curves for samples with square holes with a fillet radius R2 and with reinforcement (no. 5) and without composite overlays (no. 3) are compared. It should be noted that the Yield limit in sample no. 3 without reinforcement was achieved at εTOT = 0.32% (see Figure 15; it is corresponding to point 2 in Figure 14). The strain field on the composite overlays for such a loading level is visible in Figure 16c and can be compared with Figure 11a for the similar total strain for the sample without reinforcements. In the first part of the stress-strain curve ( Figure 15 for εTOT < 0.6%), which is the subject of the presented analysis, the benefits of the application of composite reinforcements are visible. The first differences in stress-strain curves can be observed for stresses greater than 325 MPa. The application of the overlays significantly reduces plastic stresses in a metal part. However, it is visible that above 350 MPa the plastic deformations occur in the metal part also in sample no. 5 with overlays. In the presented example, the influence of the composite overlays on the sample behavior begins to be visible, when plasticization of the metal part exceeds some critical level (above stress 300 MPa). The sudden drop of the stress and strain is observed after the separation of the Some portion of the loading is carried out by the composite overlays through the adhesive joint.
In the obtained curve ( Figure 14) a few interesting points can be observed. Till point 1, the curve is linear. The stress in point 1 corresponds to the plasticizing of the metal sample no. 3 in areas with stress concentrations. Point 1 corresponds to the strain in which the theoretical value of the maximal stress in the notch is equal to the Yield limit of the material. This stress was calculated as follows: It can be observed that below this value the relationship between stress and strain is linear. The significant increase in the curvature of the σ-ε curve is visible for stresses above 400 MPa. This effect is clearly visible in Figure 15, where stress-strain curves for samples with square holes with a fillet radius R2 and with reinforcement (no. 5) and without composite overlays (no. 3) are compared. It should be noted that the Yield limit in sample no. 3 without reinforcement was achieved at ε TOT = 0.32% (see Figure 15; it is corresponding to point 2 in Figure 14). The strain field on the composite overlays for such a loading level is visible in Figure 16c and can be compared with Figure 11a for the similar total strain for the sample without reinforcements.
In the first part of the stress-strain curve ( Figure 15 for ε TOT < 0.6%), which is the subject of the presented analysis, the benefits of the application of composite reinforcements are visible. The first differences in stress-strain curves can be observed for stresses greater than 325 MPa. The application of the overlays significantly reduces plastic stresses in a metal part. However, it is visible that above 350 MPa the plastic deformations occur in the metal part also in sample no. 5 with overlays. In the presented example, the influence of the composite overlays on the sample behavior begins to be visible, when plasticization of the metal part exceeds some critical level (above stress 300 MPa). The sudden drop of the stress and strain is observed after the separation of the particular overlays (see magnified zones in Figure 15). After failures of both overlays (caused by failure of the adhesive joint), a strong increase in the plastic strains in a metal part is observed in the σ-ε curve. This means that above 400 MPa, the overlays carry the greater part of the growing tension. The total plastic deformation of the metal part can be observed after the separation of both overlays (for strain ε TOT = 0.65%). Both failures of reinforcement overlays are caused by the exceeding of shear stresses in adhesive joints.
It can be observed that below this value the relationship between stress and strain is linear. The significant increase in the curvature of the σ-ε curve is visible for stresses above 400 MPa. This effect is clearly visible in Figure 15, where stress-strain curves for samples with square holes with a fillet radius R2 and with reinforcement (no. 5) and without composite overlays (no. 3) are compared. It should be noted that the Yield limit in sample no. 3 without reinforcement was achieved at εTOT = 0.32% (see Figure 15; it is corresponding to point 2 in Figure 14). The strain field on the composite overlays for such a loading level is visible in Figure 16c and can be compared with Figure 11a for the similar total strain for the sample without reinforcements. In the first part of the stress-strain curve ( Figure 15 for εTOT < 0.6%), which is the subject of the presented analysis, the benefits of the application of composite reinforcements are visible. The first differences in stress-strain curves can be observed for stresses greater than 325 MPa. The application of the overlays significantly reduces plastic stresses in a metal part. However, it is visible that above 350 MPa the plastic deformations occur in the metal part also in sample no. 5 with overlays. In the presented example, the influence of the composite overlays on the sample behavior begins to be visible, when plasticization of the metal part exceeds some critical level (above stress 300 MPa). The sudden drop of the stress and strain is observed after the separation of the  Figure 15). After failures of both overlays (caused by failure of the adhesive joint), a strong increase in the plastic strains in a metal part is observed in the σ-ε curve. This means that above 400 MPa, the overlays carry the greater part of the growing tension. The total plastic deformation of the metal part can be observed after the separation of both overlays (for strain εTOT = 0.65%). Both failures of reinforcement overlays are caused by the exceeding of shear stresses in adhesive joints. The first cracks in the adhesive joint are observed in the photographs for total strain εTOT = 0.5% (point 3 in Figure 14). This can be also observed in the results from DIC in Figures 17 and 18. The maximal stress for the sample with overlays is equal to 453 MPa (p. 4 in Figure 14). For such loading conditions, the first overlays are separated from the metal part. After the failure of the first overlay, the stress is reduced to 406 MPa. The second overlay is separated for stress equal to 422 MPa, and after that, the whole loading is carried only by the metal part. The first cracks in the adhesive joint are observed in the photographs for total strain ε TOT = 0.5% (point 3 in Figure 14). This can be also observed in the results from DIC in Figures 17 and 18. The maximal stress for the sample with overlays is equal to 453 MPa (p. 4 in Figure 14). For such loading conditions, the first overlays are separated from the metal part. After the failure of the first overlay, the stress is reduced to 406 MPa. The second overlay is separated for stress equal to 422 MPa, and after that, the whole loading is carried only by the metal part.
A significant role in the failure of the sample with composite reinforcements plays stress concentrations. In this case, the high strain concentrations occurred in: 1.
The metal part-four points (not visible in presented attached Figures, but can be assessed in the example of the plate without reinforcements- Figure 11); 2.
The composite overlays/adhesive joint-three points on each vertical edge of the rectangular hole; two points in the rounded corner and one in the middle length of the hole (Figure 16; Figure 17); 3.
The adhesive joint-four points at the corners of the reinforcement overlays (Figures 16-18).  Due to the fact that the concentrations in the metal part are not determined in the DIC analyses, the further comparison is limited to the strains in composite overlays and strains achieved in adhesive joints. For the lower tensile load (total strain εTOT < 0.4%- Figure 16), the higher strains occurred in the composite overlays around the rectangular hole. For higher strains, the strain concentrations are the highest in the adhesive joints ( Figure 17). Such concentrations grow from each end of rectangular reinforcement and the increase of these strains is significantly greater than the increase in the strains around the rectangular hole. This is probably caused by the decrease in the stiffness of the metal part (due to the plasticization in the range of the yield plateau and strain hardening) placed between the overlays. This leads to the situation in which the stiffness of the composite overlays is higher than the stiffness of the metal part. The further growth of the external load is redistributed between the steel panel and composite overlays and this redistribution depends on the tension stiffness of the respective parts. The straight assessment of this phenomenon is difficult to evaluate due to the spread of the plastic, weakened zone in steel plate, the size of which strongly influences the final tension stiffness of the steel part. It should be also noted that in this case, the adhesive joint is the weakest part of the sample.
On the basis of the taken photographs and the performed DIC analysis, it is difficult to determine what part of the overlay carried the load. However, it seems that both ends of the overlays (5 mm length from upper and bottom overlay end-see Figure 17 and areas on the overlays with dark blue color) are much less loaded. This is probably caused by the described-above effect of the reduction in the stiffness of the metal part. Let us compare distributions of strains in two samples: with composite overlays (Figure 17) and sample 3 without overlays (Figure 11). It can be seen that the area in which strains occurred on the overlay surface is similar to the area in which large plastic deformation Due to the fact that the concentrations in the metal part are not determined in the DIC analyses, the further comparison is limited to the strains in composite overlays and strains achieved in adhesive joints. For the lower tensile load (total strain ε TOT < 0.4%- Figure 16), the higher strains occurred in the composite overlays around the rectangular hole. For higher strains, the strain concentrations are the highest in the adhesive joints ( Figure 17). Such concentrations grow from each end of rectangular reinforcement and the increase of these strains is significantly greater than the increase in the strains around the rectangular hole. This is probably caused by the decrease in the stiffness of the metal part (due to the plasticization in the range of the yield plateau and strain hardening) placed between the overlays. This leads to the situation in which the stiffness of the composite overlays is higher than the stiffness of the metal part. The further growth of the external load is redistributed between the steel panel and composite overlays and this redistribution depends on the tension stiffness of the respective parts. The straight assessment of this phenomenon is difficult to evaluate due to the spread of the plastic, weakened zone in steel plate, the size of which strongly influences the final tension stiffness of the steel part. It should be also noted that in this case, the adhesive joint is the weakest part of the sample.
On the basis of the taken photographs and the performed DIC analysis, it is difficult to determine what part of the overlay carried the load. However, it seems that both ends of the overlays (5 mm length from upper and bottom overlay end-see Figure 17 and areas on the overlays with dark blue color) are much less loaded. This is probably caused by the described-above effect of the reduction in the stiffness of the metal part. Let us compare distributions of strains in two samples: with composite overlays ( Figure 17) and sample 3 without overlays ( Figure 11). It can be seen that the area in which strains occurred on the overlay surface is similar to the area in which large plastic deformation occurs in the plate without overlays (Figure 19). The final failure form of sample no. 5 with composite overlays is given in Figure 20.

Finite Element Analyses
The lack of the respective analytical solutions in the studied cases can be compensated by the numerical analysis. For that purpose, the finite element approach has been involved, and the ANSYS software [61] was used. The respective numerical analysis can be performed on various levels of approximation, in the general case the structure can be regarded as the 3D object, which demands the application of 3D solid finite elements. The main drawback of such an approach is the size of the numerical task, which should be solved to get reliable results. The reason for this is the need for the application of very dense mesh in the vicinity of notches to precisely predict and illustrate the stress concentration effect. The compromise in this is the application of the 2D shell or plane models, which are fully justified due to the small, in comparison with the remaining dimensions, thickness of the analyzed plate. In general, the ANSYS software offers different types of finite elements, which can be used for the analysis of orthotropic, layered laminates or their combinations with isotropic materials. Among them, there are the two finite strain shell elements, two layered structural solid elements and 3D layered structural shell elements, which allows for the inclusion of laminates with even 250 independently defined layers. The essential difference between the above-mentioned elements relies on the approximation level of the searched displacements. These can be modeled by polynomials with linear or quadratic shape functions. Here, in this paper, the finite strain shell elements with linear approximation (SHELL181) were used for the calculations. The main reason for this was because the thickness of the analyzed structure was several times smaller than the remaining dimension of the plate (W/t > 10). The detailed numerical calculations were performed for the formerly investigated experimentally flat samples with the circular hole or the rectangular hole with rounded corners with a fillet radius R2. The obtained results were compared between the plates without overlays and with overlays. The three types of overlays are of common use. These were: overlays in the shape of the circular patch, overlays of the rectangular shape and overlays in the form of

Finite Element Analyses
The lack of the respective analytical solutions in the studied cases can be compensated by the numerical analysis. For that purpose, the finite element approach has been involved, and the ANSYS software [61] was used. The respective numerical analysis can be performed on various levels of approximation, in the general case the structure can be regarded as the 3D object, which demands the application of 3D solid finite elements. The main drawback of such an approach is the size of the numerical task, which should be solved to get reliable results. The reason for this is the need for the application of very dense mesh in the vicinity of notches to precisely predict and illustrate the stress concentration effect. The compromise in this is the application of the 2D shell or plane models, which are fully justified due to the small, in comparison with the remaining dimensions, thickness of the analyzed plate. In general, the ANSYS software offers different types of finite elements, which can be used for the analysis of orthotropic, layered laminates or their combinations with isotropic materials. Among them, there are the two finite strain shell elements, two layered structural solid elements and 3D layered structural shell elements, which allows for the inclusion of laminates with even 250 independently defined layers. The essential difference between the above-mentioned elements relies on the approximation level of the searched displacements. These can be modeled by polynomials with linear or quadratic shape functions. Here, in this paper, the finite strain shell elements with linear approximation (SHELL181) were used for the calculations. The main reason for this was because the thickness of the analyzed structure was several times smaller than the remaining dimension of the plate (W/t > 10). The detailed numerical calculations were performed for the formerly investigated experimentally flat samples with the circular hole or the rectangular hole with rounded corners with a fillet radius R2. The obtained results were compared between the plates without overlays and with overlays. The three types of overlays are of common use. These were: overlays in the

Finite Element Analyses
The lack of the respective analytical solutions in the studied cases can be compensated by the numerical analysis. For that purpose, the finite element approach has been involved, and the ANSYS software [61] was used. The respective numerical analysis can be performed on various levels of approximation, in the general case the structure can be regarded as the 3D object, which demands the application of 3D solid finite elements. The main drawback of such an approach is the size of the numerical task, which should be solved to get reliable results. The reason for this is the need for the application of very dense mesh in the vicinity of notches to precisely predict and illustrate the stress concentration effect. The compromise in this is the application of the 2D shell or plane models, which are fully justified due to the small, in comparison with the remaining dimensions, thickness of the analyzed plate. In general, the ANSYS software offers different types of finite elements, which can be used for the analysis of orthotropic, layered laminates or their combinations with isotropic materials. Among them, there are the two finite strain shell elements, two layered structural solid elements and 3D layered structural shell elements, which allows for the inclusion of laminates with even 250 independently defined layers. The essential difference between the above-mentioned elements relies on the approximation level of the searched displacements. These can be modeled by polynomials with linear or quadratic shape functions. Here, in this paper, the finite strain shell elements with linear approximation (SHELL181) were used for the calculations. The main reason for this was because the thickness of the analyzed structure was several times smaller than the remaining dimension of the plate (W/t > 10). The detailed numerical calculations were performed for the formerly investigated experimentally flat samples with the circular hole or the rectangular hole with rounded corners with a fillet radius R2. The obtained results were compared between the plates without overlays and with overlays. The three types of overlays are of common use. These were: overlays in the shape of the circular patch, overlays of the rectangular shape and overlays in the form of two thin stripes located by the side of the cut hole directed along the tension direction. The first two shapes seem more justified for structures with local voids, like holes, while the last one is more suitable for long slabs with cracks or elements with external undercuts. The finite element meshes used in the calculations exhibited the common feature, which was as follows-the use of the relatively dense meshes with small elements around the notches and in their neighborhood. This allows obtaining a precise approximation of the strain and stress values, in the notches and their nearby zones and the application of rather moderate meshes in the areas, which are not the main interest. The quality of the applied approximation is controlled by the h-convergence study in notches, observation of regularity of the strain or stress contours, and the control of the agreement of the results in symmetrically located points, with respect to the vertical middle cross-section of the samples. This is valid, of course, for cases with no composite overlays or for samples with overlays with fiber orientation θ = 0 • or 90 • with respect to the tension direction. The exemplary mesh for the plate with the rectangular hole with rounded corners is presented in Figure 21. The analogous meshes were used for the remaining notches. cuts. The finite element meshes used in the calculations exhibited the common feature, which was as follows-the use of the relatively dense meshes with small elements around the notches and in their neighborhood. This allows obtaining a precise approximation of the strain and stress values, in the notches and their nearby zones and the application of rather moderate meshes in the areas, which are not the main interest. The quality of the applied approximation is controlled by the h-convergence study in notches, observation of regularity of the strain or stress contours, and the control of the agreement of the results in symmetrically located points, with respect to the vertical middle cross-section of the samples. This is valid, of course, for cases with no composite overlays or for samples with overlays with fiber orientation θ = 0° or 90° with respect to the tension direction. The exemplary mesh for the plate with the rectangular hole with rounded corners is presented in Figure 21. The analogous meshes were used for the remaining notches. The detailed numerical calculations were performed for the stacking sequence as follows (see Figure 3b): the layer with the number 9 corresponds to the steel plate, while the layers from 1 to 8, and 10 to 17 are the fiber-reinforced composite layers, with the respective angle orientation. The FE analysis was performed for the flat bar with the length of 200 mm, the width of 45 mm and the thickness of 4 mm (Figure 2). In such a flat element the respective rectangular or circular holes were made. The sizes of the openings were chosen in such a way that the nominal stress in the weakened cross-section was equal to 8.33 MPa, which meant that the diameter of the circular hole was set to 15 mm and for the rectangular void, the distance between the opposite straight sides was assumed to be the same. The samples were subjected to the tension force F = 1000 N in the vertical direction at the bottom, while the upper edge was fully clamped. Such a configuration fully followed the experimental tests. In the first step, the samples with no overlays were analyzed to get the basis for the stress concentration comparisons. Figure 22 illustrates the numerically obtained stress concentrations observed in the weakened hole sample cross-section. The detailed numerical calculations were performed for the stacking sequence as follows (see Figure 3b): the layer with the number 9 corresponds to the steel plate, while the layers from 1 to 8, and 10 to 17 are the fiber-reinforced composite layers, with the respective angle orientation. The FE analysis was performed for the flat bar with the length of 200 mm, the width of 45 mm and the thickness of 4 mm (Figure 2). In such a flat element the respective rectangular or circular holes were made. The sizes of the openings were chosen in such a way that the nominal stress in the weakened cross-section was equal to 8.33 MPa, which meant that the diameter of the circular hole was set to 15 mm and for the rectangular void, the distance between the opposite straight sides was assumed to be the same. The samples were subjected to the tension force F = 1000 N in the vertical direction at the bottom, while the upper edge was fully clamped. Such a configuration fully followed the experimental tests. In the first step, the samples with no overlays were analyzed to get the basis for the stress concentration comparisons. Figure 22 illustrates the numerically obtained stress concentrations observed in the weakened hole sample cross-section. The above results are considered the reference solutions taken for comparison in further assessment of the applied overlay effectiveness. Due to the fact that numerous investigations have been made for plates with circular holes, the presented study is concentrated on the experimental and numerical investigations for the samples with rectangular holes.

Circular Patch
In Figures 23-25, the stress distributions in the analyzed samples with composite overlays are presented. This set concerns the results obtained for the circular patch placed around the rectangular hole. In Figure 23, the distribution of tension stress is presented for the radius of the overlay set to Rov = 20 mm and for the fiber orientation θ=0˚. Such a combination provided the minimum stress concentration for the studied sample and overlay configurations. In Figure 24 the comparison of the tension stresses is shown for different fiber orientation angles and for Rov = 20 mm. These stresses are shown in the area just below the circular patch. In Figure 25 the respective tension stresses in the top layer of the overlay are shown.
Due to the fact that the points with the maximum value of the tension stress occupy different positions around the circumference of the void (these points are always located in the notch area), the plots are constructed in such a manner that point C always corresponds to the point with the maximum stress value, but point D is the point located on the outer edge of the sample with the vertical coordinate equal to the 'y' coordinate of point C (Figure 26). In Figure 27a the distribution of tension stress along the line CD is shown. Here the results are given for the θ = 0° fiber orientation and for different sizes of the circular patch. Such fiber orientation provides the maximal reduction of the stress concentration. In Figure 27b the results of tension stress along the CD line are presented for different fiber orientations and for the set value of the outer radius of the circular patch (Rov = 17.5 mm). It is visible that the orientation of fibers along the tension direction in overlays provides the maximum stress reduction. The above results are considered the reference solutions taken for comparison in further assessment of the applied overlay effectiveness. Due to the fact that numerous investigations have been made for plates with circular holes, the presented study is concentrated on the experimental and numerical investigations for the samples with rectangular holes.

Circular Patch
In Figures 23-25, the stress distributions in the analyzed samples with composite overlays are presented. This set concerns the results obtained for the circular patch placed around the rectangular hole. In Figure 23, the distribution of tension stress is presented for the radius of the overlay set to R ov = 20 mm and for the fiber orientation θ = 0 • . Such a combination provided the minimum stress concentration for the studied sample and overlay configurations. In Figure 24 the comparison of the tension stresses is shown for different fiber orientation angles and for R ov = 20 mm. These stresses are shown in the area just below the circular patch. In Figure 25 the respective tension stresses in the top layer of the overlay are shown.
Due to the fact that the points with the maximum value of the tension stress occupy different positions around the circumference of the void (these points are always located in the notch area), the plots are constructed in such a manner that point C always corresponds to the point with the maximum stress value, but point D is the point located on the outer edge of the sample with the vertical coordinate equal to the 'y' coordinate of point C ( Figure 26). In Figure 27a the distribution of tension stress along the line CD is shown. Here the results are given for the θ = 0 • fiber orientation and for different sizes of the circular patch. Such fiber orientation provides the maximal reduction of the stress concentration. In Figure 27b the results of tension stress along the CD line are presented for different fiber orientations and for the set value of the outer radius of the circular patch (R ov = 17.5 mm). It is visible that the orientation of fibers along the tension direction in overlays provides the maximum stress reduction.     Figure 28 is not to scale drawing of the mesh of finite elements for the square patch, while in Figures 29-31, the stress distributions in analyzed samples with composite overlays are presented. This set concerns the results obtained for the square patch placed around the rectangular hole. The distribution of tensile stress in the steel part is presented for different fiber orientations are presented in Figure 29. As previously, the θ = 0° fiber orientation provides the maximum stress reduction. In Figure 30 the respective distributions of vertical stresses on the outer composite layer are presented. Figure 31 shows the distribution of the vertical stresses in the most strenuous weakened cross-section of the steel panel.  Figure 28 is not to scale drawing of the mesh of finite elements for the square patch, while in Figures 29-31, the stress distributions in analyzed samples with composite overlays are presented. This set concerns the results obtained for the square patch placed around the rectangular hole. The distribution of tensile stress in the steel part is presented for different fiber orientations are presented in Figure 29. As previously, the θ = 0° fiber orientation provides the maximum stress reduction. In Figure 30 the respective distributions of vertical stresses on the outer composite layer are presented. Figure 31 shows the distribution of the vertical stresses in the most strenuous weakened cross-section of the steel panel.  Figure 28 is not to scale drawing of the mesh of finite elements for the square patch, while in Figures 29-31, the stress distributions in analyzed samples with composite overlays are presented. This set concerns the results obtained for the square patch placed around the rectangular hole. The distribution of tensile stress in the steel part is presented for different fiber orientations are presented in Figure 29. As previously, the θ = 0 • fiber orientation provides the maximum stress reduction. In Figure 30 the respective distributions of vertical stresses on the outer composite layer are presented. Figure 31 shows the distribution of the vertical stresses in the most strenuous weakened cross-section of the steel panel.   A summary of the obtained numerical results is given in Table 5. The analysis of them says that the bigger size of the overlay gives the lower stress concentration in the notch area in order that the maximum stress reduction is obtained for the square patch, which covers the full width of the sample. Additionally, the most convenient stacking sequence is the one with fibers placed parallel to the tension load direction. This conclusion is obvious for the analyzed set of loadings, but fails in the general case, where the system of external loads becomes more complex. In such a situation the best orientation of fibers should be set in the optimization process. The tension stresses in the vertical direction obtained for the outer composite layer are rather low. This observation joined with a rather moderate reduction of the maximum stress in a steel notched part, means that the applied overlay are too stiff and their tension stiffness should be reduced. This observation opens an area for the search of optimal-in a view of stress reduction-choice of the overlays, including material of fibers, number of applied layers, fiber orientation, and stacking sequence, which in sum define the overlay stiffness.

Square Patch
In this set of analyzed problems, it is worth mentioning that the maximum values for the tension stress in the notch appearing in the steel part are observed in the different rounded corners, but the maximum values for the remaining three notches remain almost the same within the calculation error (the software shows the maximum sign only for one point, even if the values in the three remaining points are the same).

Discussion
The main aim of this study was to attempt to increase the strength of structures with cut-outs. In the presented preliminary studies, this has been done through the application of the composite overlays on both sides of the metal sample with the cut-out. Such reinforcement by the composite overlays is investigated in the paper in the example of the sample with the square hole with a fillet radius R2. The influence of the composite reinforcement on the sample strength is presented in Figure 32. All the figures present nominal stresses in weakened cross-sections in the metal part. It can be seen that the strength of the sample with reinforcements is higher than the strength of the other tested samples 2-4 with circular and square holes. It should also be noted that the experimentally tested sample no. 5 with composite reinforcements has the cut-out with the largest stress concentration factor K t = 2.479 (Table 4) in comparison with other cut-outs (circular hole-K t = 2.314 and square with filler radius R4-K t = 2.084). Despite, the highest K t , the composite reinforcements significantly increased the strength of the structure in comparison with samples with smaller stress concentrations ( Figure 32). of the composite overlays on both sides of the metal sample with the cut-out. Such reinforcement by the composite overlays is investigated in the paper in the example of the sample with the square hole with a fillet radius R2. The influence of the composite reinforcement on the sample strength is presented in Figure 32. All the figures present nominal stresses in weakened cross-sections in the metal part. It can be seen that the strength of the sample with reinforcements is higher than the strength of the other tested samples 2-4 with circular and square holes. It should also be noted that the experimentally tested sample no. 5 with composite reinforcements has the cut-out with the largest stress concentration factor Kt = 2.479 (Table 4) in comparison with other cut-outs (circular hole-Kt = 2.314 and square with filler radius R4-Kt = 2.084). Despite, the highest Kt, the composite reinforcements significantly increased the strength of the structure in comparison with samples with smaller stress concentrations ( Figure 32). Comparing the obtained results to the Yield limit of sample no. 3, the maximal stress is increased from 365 MPa (sample no. 3 without reinforcement) to 454 MPa (sample with composite reinforcements). This difference of Δσ = 89 MPa means that the load-carrying capacity of the sample with reinforcements is increased by almost 25% before the damage of the overlays.
Taking into account the maximal tensile force at the moment of the failure of the adhesive joint, the minimal average shear stress in the adhesive joint is equal to: This value is significantly lower than the minimum shear strength of the adhesive joint, but it should be noted that in the adhesive joint there are large stress concentrations. Such stress concentrations are around the hole and at the ends of the adhesive joints. The value of the stress concentrations at the end of the bounded joint often achieves high values [40,41,45,51]. This value depends on the stiffness of the adherents, the adhesive shear modulus and thickness, and the geometry of the joint (overlap width and length). The distribution of the shear stresses is calculated using FEM. It allows determination of the stress concentrations around the hole. However, the determination of the stress concentrations at the ends of the adhesive is impossible with the use of the applied numerical model. This problem requires further analyses with the application of the 3D FEM analyses.
The results are presented for the maximal tensile loading (εTOT = 0.57%) just before separation of the first overlay (p. 4 in Figure 14). In Figure 33a the distribution of the Comparing the obtained results to the Yield limit of sample no. 3, the maximal stress is increased from 365 MPa (sample no. 3 without reinforcement) to 454 MPa (sample with composite reinforcements). This difference of ∆σ = 89 MPa means that the load-carrying capacity of the sample with reinforcements is increased by almost 25% before the damage of the overlays.
Taking into account the maximal tensile force at the moment of the failure of the adhesive joint, the minimal average shear stress in the adhesive joint is equal to: This value is significantly lower than the minimum shear strength of the adhesive joint, but it should be noted that in the adhesive joint there are large stress concentrations. Such stress concentrations are around the hole and at the ends of the adhesive joints. The value of the stress concentrations at the end of the bounded joint often achieves high values [40,41,45,51]. This value depends on the stiffness of the adherents, the adhesive shear modulus and thickness, and the geometry of the joint (overlap width and length). The distribution of the shear stresses is calculated using FEM. It allows determination of the stress concentrations around the hole. However, the determination of the stress concentrations at the ends of the adhesive is impossible with the use of the applied numerical model. This problem requires further analyses with the application of the 3D FEM analyses.
The results are presented for the maximal tensile loading (ε TOT = 0.57%) just before separation of the first overlay (p. 4 in Figure 14). In Figure 33a the distribution of the vertical strain on the external composite layer is presented. Good agreement of DIC ( Figure 17d) and FEM (Figure 33a) analyses is achieved. The shear stress in the adhesive is presented in Figure 33b. High concentrations of stresses can be observed in the surroundings of the notches of the rectangular hole in both composite overlays (Figure 33a) and in adhesive epoxy (Figure 33b). The value of shear stresses around the notch significantly exceeds the admissible value of 27 MPa. This justifies further optimization of the adhesive joint consisting in determining the optimal design of the adhesive joint (overlay shape and size, orientation of the layers, adhesive type and thickness, and shape of the hole). vertical strain on the external composite layer is presented. Good agreement of DIC ( Figure 17d) and FEM (Figure 33a) analyses is achieved. The shear stress in the adhesive is presented in Figure 33b. High concentrations of stresses can be observed in the surroundings of the notches of the rectangular hole in both composite overlays ( Figure 33a) and in adhesive epoxy (Figure 33b). The value of shear stresses around the notch significantly exceeds the admissible value of 27 MPa. This justifies further optimization of the adhesive joint consisting in determining the optimal design of the adhesive joint (overlay shape and size, orientation of the layers, adhesive type and thickness, and shape of the hole). The two typical shapes of the composite overlays were studied and the influence of size, shape and different fiber orientations were numerically investigated. The finite element mesh was chosen in such a way that the highest density of the elements was introduced in the area neighboring notches. The generated element meshes over the whole sample mimics the symmetry of the structure. However, due to the application of the composite layers and their fiber orientations, the symmetry conditions were not used for all analyzed cases. The obtained result show the reduction of stress concentrations in the hole corners after application of the overlays. The most convenient fiber orientation for both studied corner radii is that one, in which the fibers were set along the direction of the applied external tension load (θ = 0°). In the remaining fiber orientations, the influence of the overlays on the stress reduction was smaller. It appeared that the bigger the size of the overlay, the higher the reduction of the maximum stress is observed. In this context, the square overlay appeared to be the most effective, giving 20% of the maximum stress reduction for fiber orientation θ = 0°. However one must be aware that in the practical application, the size of the overlay may be limited by different factors, such as the general size of the object, etc. The summary of the effectiveness of the studied overlays can be expressed in the percentage reduction of the stress concentration factor in the steel plate part (see Table 6). If the composite overlays are to be applied, then the proper choice of overlays data becomes more complex and must include the influence of the number of layers, the The two typical shapes of the composite overlays were studied and the influence of size, shape and different fiber orientations were numerically investigated. The finite element mesh was chosen in such a way that the highest density of the elements was introduced in the area neighboring notches. The generated element meshes over the whole sample mimics the symmetry of the structure. However, due to the application of the composite layers and their fiber orientations, the symmetry conditions were not used for all analyzed cases. The obtained result show the reduction of stress concentrations in the hole corners after application of the overlays. The most convenient fiber orientation for both studied corner radii is that one, in which the fibers were set along the direction of the applied external tension load (θ = 0 • ). In the remaining fiber orientations, the influence of the overlays on the stress reduction was smaller. It appeared that the bigger the size of the overlay, the higher the reduction of the maximum stress is observed. In this context, the square overlay appeared to be the most effective, giving 20% of the maximum stress reduction for fiber orientation θ = 0 • . However one must be aware that in the practical application, the size of the overlay may be limited by different factors, such as the general size of the object, etc. The summary of the effectiveness of the studied overlays can be expressed in the percentage reduction of the stress concentration factor in the steel plate part (see Table 6). If the composite overlays are to be applied, then the proper choice of overlays data becomes more complex and must include the influence of the number of layers, the stacking sequence of them and the size of the overlay. The numerical analyses performed in the paper relied on the use of layered shell elements. In such an analysis, it is difficult to model and study the influence of the applied bonding layer (glue layer). In the further analysis, it is planned to perform full 3D analysis with the inclusion of the mechanical properties of the applied bonding layer. Here an additional problem would concern the choice of the kind of adhesive, its thickness, and the mechanical properties. It is also planned to perform the investigations for other radii notches (for rectangular hole) and for different shapes of notches used in practice.

Conclusions
The performed experimental and numerical study let us draw the following conclusions: • the good agreement between the DIC results with the FEM results has confirmed the usefulness of the DIC tool and GOM software in the non-destructive analysis of structural elements; • the effectiveness of the composite overlay application has been confirmed by the reduction of the stress concentration factor in the notch area. In the investigated example, the application of the square composite overlay increased the structure strength even by 25%; • the stress concentration reduction is the function of the fiber orientation and is the highest in the case of fibers set in parallel to the applied loading; • the bigger size of the overlay patch the bigger reduction of the stress concentration is observed.