Next Article in Journal
Experimental Investigation of Thrust Force, Delamination and Surface Roughness in Drilling Hybrid Structural Composites
Next Article in Special Issue
The Effect of Prestressing and Temperature on Tensile Strength of Basalt Fiber-Reinforced Plywood
Previous Article in Journal
Predicting Compression Pressure of Knitted Fabric Using a Modified Laplace’s Law
Previous Article in Special Issue
Determination of Mechanical and Tribological Properties of Silicone-Based Composites Filled with Manganese Waste
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation of the Near-Tip Stress Field of a Notch Terminating at a Bi-Material Interface

1
Faculty of Mechanical Engineering, Bialystok University of Technology, 45C Wiejska Str., 15-351 Bialystok, Poland
2
Faculty of Engineering, Mansoura University, Mansoura 35516, Egypt
3
Faculty of Science, Mansoura University, Mansoura 35516, Egypt
*
Author to whom correspondence should be addressed.
Materials 2021, 14(16), 4466; https://doi.org/10.3390/ma14164466
Submission received: 29 June 2021 / Revised: 29 July 2021 / Accepted: 5 August 2021 / Published: 9 August 2021

Abstract

:
The article deals with the problem of a sharp corner, the tip of which is located on the bi-material interface. The paper presents a qualitative and quantitative description of singular stress fields occurring in the tip area of such a stress concentrator. The qualitative description was obtained by solving the problem of the plane theory of elasticity with appropriately defined boundary conditions. To obtain a quantitative description, it was necessary to determine the values of generalised stress intensity factors (GSIFs). The GSIFs were determined using the developed analytical-numerical method. The calculations were made for various load variants (uniaxial/biaxial tension load, shear load) and notch positions (single/double edge-notched plate, centre-notched plate). Additionally, the impact of notch geometry (height and opening angle) and relative stiffness (Young’s moduli ratio of both components of bi-material) on GSIFs was investigated. It has been noticed that with a decrease in the relative stiffness and an increase in the notch angle or its height, the normalised GSIFs values increased. The obtained results were compared with the data available in the literature and their satisfactory agreement with those presented by other scientists was found.

1. Introduction

Ensuring the high durability of the structure with minimum costs is a priority of today’s world economy. The durability of the structure largely depends on the strength of the materials used for its components.
Currently, various types of composites are widely used, which, compared to traditional construction materials, are characterised by greater strength and at the same time lower specific weight. Composites are not without some drawbacks. They often contain various material defects (voids, inclusions) causing the formation of local stress fields with large gradients. This results in the initiation of new cracks or the propagation of the existing ones. In the case where the stress concentrator is located inside only one material phase, the failure criteria commonly used for isotropic materials can be used to predict fracture [1,2,3,4]. In the opposite situation, e.g., when the defect is located or started at the surface of the connection of different elastic materials, it is necessary to develop new, or modify already existing, material fracture criteria. Suitable fracture hypotheses can be obtained using the qualitative and quantitative description of singular stress fields occurring in close proximity to the tip of the defect under consideration.
Many scientists have dealt with the analytical description of local mechanical fields generated by defects located at the interface (e.g., crack [5,6], structural notch [7,8], and rigid inclusion [9,10,11]). In the literature, much less attention has been paid to the problem of a notch whose tip is located at the interface. Such a defect may appear as a result of the development of voids located in one of the phases of the composite [12,13,14]. It is worth noting that a special case of such a sharp corner (notch with a zero-opening angle) is the crack initiated at the interface [15,16,17,18,19,20,21,22,23,24]. Stress singularities occurring in the tip area of such crack, which is perpendicular to the interface were analysed in [15]. GSIFs values, determined for a short crack located in a bi-material with infinite and finite overall dimensions, are presented in [16] and [17,18,19,20,21,22,23], respectively. The problem of local stress concentration, occurring in the vicinity of the crack tip, not perpendicular to the interface, was discussed in [24,25]. As for the corners with a non-zero opening angle, the authors of the papers [12,26] have dealt with this subject. In [26], eigenequation was derived and the influence of material constants on the roots of this equation was investigated. Moreover, the authors, assuming that the notched bi-material is subjected to tensile loads, determined generalised stress intensity factors. An approximate analytical model, based on the theory of multilayer beams, enabling the calculation of GSIF (for mode I loading) is also presented in [12].
There is no insufficient information in the literature on the complete description of the singular stress fields generated by a notch whose tip is at the interface of two elastic bodies. Therefore, the main goal of the work was to obtain such analytical and numerical solutions. The obtained analytical descriptions (eigenequation, formulas for individual components of the stress tensor) along with the proposed method of GSIFs determination are presented in the first part of the presented work. The values of the generalised stress intensity factors, determined for various load variants and geometric and material properties of the composite, are presented in the second part of the article.

2. Materials and Methods

2.1. Analytical Solutions

The analytical description of the stress fields was obtained by solving (with the accuracy of multiplicative constants called GSIFs) a plane problem of two connected elastic half-spaces, the interface of which is weakened by a sharp corner (Figure 1).
Consider a sharp corner (Figure 1) with a polar coordinate system (with coordinates r and φ) which is located at the notch tip. In such a reference system, there are two independent components of the displacement vector, in the radial (ur) and angular direction (uφ).
The related strain components have the following form:
ε r = u r r ,   ε φ = 1 r u φ φ + u r r ,   ε r φ = 1 r u r φ + u φ r u φ r .
The generalised Hooke’s law can be written as:
σ r = Λ ( ε r + ε φ ) + 2 μ ε r ,   σ φ = Λ ( ε r + ε φ ) + 2 μ ε φ ,   τ r φ = μ ε r φ  
where Λ, μ are Lame’s constants and are respectively: Λ = E ν ( 1 + ν ) ( 1 2 ν ) , μ = E 2 ( 1 + ν ) .
Navier equations are described by the Formula (3):
σ r r + 1 r τ r φ φ + σ r σ φ r = 0 ,   τ r φ r + 1 r σ φ φ + 2 τ r φ r = 0 .  
By using Formulas (1)–(3), equilibrium equations (Lame’s equation) can be obtained:
( Λ + 2 μ ) r ( u r r + 1 r u φ φ + u r r ) μ 1 r φ ( u φ r 1 r u r φ + u φ r ) = 0 ( Λ + 2 μ ) 1 r φ ( u r r + 1 r u φ φ + u r r ) + μ r ( u φ r 1 r u r φ + u φ r ) = 0 } .
The asymptotic form of the displacement function can be written as [27,28]:
u r ( r , φ ) = r λ f ( φ ) ,   u φ ( r , φ ) = r λ g ( φ ) ,  
where f(φ) and g(φ) are angular coordinate functions obtained on the basis of the Airy stress function.
Substituting (5) to Formula (4), a system of differential equations was obtained, from which the general form of the asymptotic solution is determined (6) [11,28]:
u r i = r λ j ( A i   cos ( ( 1 + λ j ) φ ) + B i   sin ( ( 1 + λ j ) φ ) + C i   cos ( ( 1 λ j ) φ ) + D i   sin ( ( 1 λ j ) φ ) ) u φ i = r λ j ( A i   sin ( ( 1 + λ j ) φ ) + B i   cos ( ( 1 + λ j ) φ ) C i κ i + λ j κ i λ j   sin ( ( 1 λ j ) φ ) + D i κ i + λ j κ i λ j   cos ( ( 1 λ j ) φ ) ) σ r i = r λ j 1 μ i ( A i 2 λ j cos ( ( 1 + λ j ) φ ) + B i 2 λ j sin ( ( 1 + λ j ) φ ) + C i ( 3 λ j ) 2 λ j κ i λ j cos ( ( 1 λ j ) φ ) + D i ( 3 λ j ) 2 λ j κ i λ j sin ( ( 1 λ j ) φ ) ) σ φ i = r λ j 1 μ i ( A i 2 λ j cos ( ( 1 + λ j ) φ ) B i 2 λ j sin ( ( 1 + λ j ) φ ) + C i ( 1 + λ j ) 2 λ j κ i λ j cos ( ( 1 λ j ) φ ) + D i ( 1 + λ j ) 2 λ j κ i λ j sin ( ( 1 λ j ) φ ) ) τ r φ i = r λ j 1 μ i ( A i 2 λ j sin ( ( 1 + λ j ) φ ) + B i 2 λ j cos ( ( 1 + λ j ) φ ) + C i ( 1 λ j ) 2 λ j κ i λ j sin ( ( 1 λ j ) φ ) D i ( 1 λ j ) 2 λ j κ i λ j cos ( ( 1 λ j ) φ ) ) } ,  
where: μ i = E i 2 ( 1 + ν i ) —shear modulus κ i = ( 3 ν i ) / ( 1 + ν i ) —a plane stress, κ i = ( 3 4 ν i ) —a plane strain, νi—Poisson’s ratio, i = 1, 2, j = I for symmetric problems (Mode I), j = II for skew-symmetric problems (Mode II).
Particular solutions for the analysed plane problem of notched bi-material were derived by determining the constants Ai, Bi, Ci, Di, and eigenvalue λj. The unknowns sought were obtained on the basis of the following boundary conditions:
  • along the interface, for φ = 0 [29];
    u r 1 = u r 2 ; u φ 1 = u φ 2 ; σ φ 1 = σ φ 2 ; τ r φ 1 = τ r φ 2 ,
  • of the upper surface of the V-notch, for φ = γ;
    σ φ 1 = τ r φ 1 = 0
  • for φ = π / 2
    • symmetry conditions (Mode I)
      τ r φ 2 = u φ 2 = 0
    • skew-symmetry conditions (Mode II)
      σ φ 2 = u r 2 = 0
Eigenequations (7-symmetric problem, 8-skew-symmetric problem) from the zero condition of the boundary condition matrix determinant were determined:
β λ I sin [ 2 γ ] + ( β 1 + 2 α ( 1 α + β + ( α β ) λ I 2 ) + 2 α ( β α ) λ I 2 cos [ 2 γ ] ) sin [ π λ I ] + + α ( 1 + α β ) sin [ ( π 2 γ ) λ I ] + ( 1 + α ) ( α β ) sin [ ( π + 2 γ ) λ I ] = 0 ,  
β λ I I sin [ 2 γ ] + ( 1 β + 2 α ( 1 + α α λ I I 2 + β ( λ I I 2 1 ) ) + 2 α ( α β ) λ I I 2 cos [ 2 γ ] ) sin [ π λ I I ] + + α ( β α 1 ) sin [ ( π 2 γ ) λ I I ] ( 1 + α ) ( α β ) sin [ ( π + 2 γ ) λ I I ] = 0 ,  
where: α = μ 1 / μ 2 1 ( 1 + κ 1 ) , β = μ 1 ( 1 + κ 2 ) μ 2 ( 1 + κ 1 ) .
The roots of the above equations correspond to the eigenvalues λj.
Assuming in Formulas (7) and (8) that α = 0 , β = 1 , γ = γ π / 2 , the resulting eigenequations are identical to those for the notch problem in isotropic material [27]:
λ I sin [ 2 γ ] + sin [ 2 γ λ I ] = 0 ,   λ I I sin [ 2 γ ] sin [ 2 γ λ I ] = 0 .
The roots of Equations (7) and (8) cannot be found analytically. They were determined numerically using the Berents method applied to the proprietary calculation program written in the Wolfram Language. The eigenvalues λj obtained in this way, determined for notches with arbitrarily assumed opening angles, are shown graphically in Figure 2.
It was found that the parameters λj, regardless of the notch geometry and the mechanical properties of the bi-material, always assume real values. Moreover, the strength of the stress singularity increases with a decrease in the notch angle and an increase in the relative stiffness.
To obtain an analytical description of individual components of the stress tensor, it was necessary to define GSIFs in advance. Since eigenvalues λj always assume real values, the generalised stress intensity factors Kj were defined, similar to the work [27], as follows:
K I = lim r 0 2 π r 1 λ I σ φ 2 ( r , π / 2 ) K I I = lim r 0 2 π r 1 λ I I τ r φ 2 ( r , π / 2 ) } .
Due to the rather complicated form of the solutions obtained, only the stress formulas in the material in which the potential crack will propagate are presented below (for the notch shown in Figure 1 it is the material marked with the number 2).
σ φ 2 = ( K I r λ I 1 ( ( 1 + λ I ) cos [ 1 2 ( λ I 1 ) ( π + 2 φ ) ] F I 2 + α μ 1 ( λ I κ 2 ) sin [ π λ I 2 + φ + λ I φ ] F I 1 ) 2 π α μ 1 ( λ I κ 2 ) F I 0 + K I I r λ I I 1 ( ( 1 + λ I I ) cos [ π λ I I 2 + ( λ I I 1 ) φ ] F I I 2 + α ( λ I I κ 2 ) cos [ π λ I I 2 + φ + λ I I φ ] F I I 1 ) 2 π α μ 1 ( λ I I κ 2 ) F I I 0 ) τ r φ 2 = ( K I r λ I 1 ( ( λ I 1 ) cos [ π λ I 2 + ( λ I 1 ) φ ] F I 2 + α μ 1 ( λ I κ 2 ) cos [ π λ I 2 + φ + λ I φ ] F I 1 ) 2 π α μ 1 ( λ I κ 2 ) F I 0 + + K I I r λ I I 1 ( ( λ I I 1 ) cos [ 1 2 ( λ I I 1 ) ( π + 2 φ ) ] F I I 2 + α ( λ I I κ 2 ) sin [ π λ I I 2 + φ + λ I I φ ] F I I 1 ) 2 π α μ 1 ( λ I I κ 2 ) F I I 0 ) σ r 2 = ( K I r λ I 1 ( ( λ I 3 ) cos [ 1 2 ( λ I 1 ) ( π + 2 φ ) ] F I 2 + α μ 1 ( λ I κ 2 ) sin [ π λ I 2 + φ + λ I φ ] F I 1 ) 2 2 π α μ 1 ( λ I κ 2 ) F I 0 + K I I r λ I I 1 ( ( λ I I 3 ) cos [ π λ I I 2 + ( λ I I 1 ) φ ] F I I 2 + α ( λ I I κ 2 ) cos [ π λ I I 2 + φ + λ I I φ ] F I I 1 ) 2 π α μ 1 ( λ I I κ 2 ) F I I 0 ) } ,  
where:
F I 1 = ( α ( λ I 1 ) λ I cos [ γ + π λ I 2 γ λ I ] + ( 1 + α β + α λ I 2 ) cos [ γ π λ I 2 + γ λ I ] + + ( β 2 α 1 ) λ I cos [ γ + π λ I 2 + γ λ I ] + ( α β ) ( λ I 1 ) cos [ γ 1 2 ( π + 2 γ ) λ I ] ) ,
F I 2 = ( α β + α λ I ) μ 1 + β μ 2 ) ( cos [ γ + π λ I 2 + γ λ I ] 2 α ( λ I sin [ γ ] sin [ 1 2 ( π 2 γ ) λ I ] + sin [ γ λ I ] sin [ γ + π λ I 2 ] ) ) ,
F I 0 = ( 2 cos [ γ ] ( λ I cos [ π λ I 2 ] cos [ γ λ I ] + ( 1 β + λ I + 2 α ( 1 + λ I ) ) sin [ π λ I 2 ] sin [ γ λ I ] ) + + sin [ γ ] ( ( 1 β + 4 α λ I 2 ) sin [ 1 2 ( π 2 γ ) λ I ] + ( 1 + β + 2 λ I + 4 α λ I 2 β λ I ) sin [ 1 2 ( π + 2 γ ) λ I ] ) ) ,
F I I 1 = μ 1 ( α λ I I ( 1 + λ I I ) cos [ γ + π λ I I 2 γ λ I I ] ( 1 + α β + α λ I I 2 ) cos [ γ π λ I I 2 + γ λ I I ] + ( β 1 2 α ) λ I I cos [ γ + π λ I I 2 + γ λ I I ] + ( α β ) ( 1 + λ I I ) cos [ γ 1 2 ( π + 2 γ ) λ I I ] ) ,
F I I 2 = ( ( α β + α λ ) μ 1 + β μ 2 ) ( ( 1 + α ) cos [ γ + π λ I I 2 + γ λ I I ] + + α λ I I cos [ γ π λ I I 2 + γ λ I I ] α ( 1 + λ I I ) cos [ γ + π λ I I 2 γ λ I I ] ) ,
F I I 0 = ( 2 cos [ γ ] ( λ I I cos [ π λ I I 2 ] cos [ γ λ I I ] ( β 1 + 2 α ( λ I I 1 ) + λ I I ) sin [ π λ I I 2 ] sin [ γ λ I I ] ) + + sin [ γ ] ( ( 1 β + 4 α λ I I 2 ) sin [ 1 2 ( π 2 γ ) λ I I ] + ( 1 + β + 2 ( β 1 2 α ) λ I I ) sin [ 1 2 ( π + 2 γ ) λ I I ] ) ) .
For the quantitative description of stresses, it is necessary to determine the values of the Kj coefficients (GSIFs). A method for finding generalised stress intensity factors is discussed in the next section.

2.2. The Method for Determining Generalised Stress Intensity Factors Kj

For the considered problem of a notch with a tip located on the bi-material interface, there are no exact solutions enabling the determination of the value of generalised stress intensity factors. The multiplied constants Kj (GSIFs) used in Formula (11) can be found using approximate methods. In this paper, the analytical and numerical method for determining GSIFs presented in the work [30] was used. It is an asymptotic method, based on the comparison of analytically and numerically obtained stress distributions in the vicinity of the defect’s tip. So, for its application, it was necessary to derive, based on the obtained analytical solutions, approximating functions and to determine the appropriate stresses using numerical methods (finite element method (FEM)). The methodology of obtaining the latter is discussed in the next section. Below, the approach of the applied method of GSIFs determination is discussed and the formulas for approximating functions are derived.
The following functions were used to extrapolate the hoop σφ2 and tangential τ2 stresses numerically determined at two neighbouring points located at a distance of rn and rn+1 from the notch tip:
σ φ 2 ( r n , π / 2 ) = K I ( r ) 2 π r 1 λ I ( 1 + c I r n ) , σ φ 2 ( r n + 1 , π / 2 ) = K I ( r ) 2 π r 1 λ I ( 1 + c I r n + 1 ) τ r φ 2 ( r n , π / 2 ) = K I I ( r ) 2 π r 1 λ I I ( 1 + c I I r n ) , τ r φ 2 ( r n + 1 , π / 2 ) = K I I ( r ) 2 π r 1 λ I I ( 1 + c I I r n + 1 ) }
where CI, CII are constants that can be eliminated from the equations.
Applying the (10) and (11) to the extrapolating functions (12), the formulas for GSIFs (approximating functions) were obtained in the following form:
K I ( r ) = 2 π ( r n r n + 1 ) 1 λ I ( r n λ I σ φ 2 ( r n + 1 , π / 2 ) r n + 1 λ I σ φ 2 ( r n , π / 2 ) ) r n r n + 1 K I I ( r ) = 2 π ( r n r n + 1 ) 1 λ I I ( r n λ I I τ r φ 2 ( r n + 1 , π / 2 ) r n + 1 λ I I τ r φ 2 ( r n , π / 2 ) ) r n r n + 1 } .
The K j ( r ) factors were calculated on the basis of the above Formula (13), assuming for σφ2 and τ2, respectively, the values of hoop and tangential stresses obtained with the use of FEM at n + 1 nodes.
The accuracy of the determined GSIFs values depends on the selection of the area from which the stresses obtained using FEM are implemented into the approximating functions. This area was determined using the criterion of selecting nodes proposed in [30].
The K j ( r ) coefficients determined in all selected nodes should theoretically be identical. However, due to potential errors in numerical calculations, the found values of the generalised stress intensity factors may differ slightly. To minimise such an error, the obtained results were averaged according to the following formula:
K I = n = 1 n + 1 K I ( r ) n + 1 , K I I = n = 1 n + 1 K I I ( r ) n + 1 .
To verify the method used, the GSIFs were calculated for the problem of a crack perpendicular to the interface and started on it. The obtained results were compared with the exact solution [16] and the approximate one [18]. A satisfactory agreement was obtained in both cases—the difference was about 1.2%. A similar convergence of results was obtained in the notch problem, the tip of which is located on the border of the connection of two materials. The relative difference between the compared GSIFs values, found with the use of the developed method and in the work [26], was about 1.8%.

2.3. FEM Modelling

Currently, various types of numerical simulations are often performed before the implementation and production of the product. Numerical calculations can be performed using various methods. One of them is FEM. This method can be used, for example, to analyse issues related to friction [31,32], flow [33,34] or to predict the operation of piezoelectric transducers [35,36,37].
In the presented work, FEM was used to determine generalised stress intensity factors. Numerical calculations were performed in the ANSYS environment, for which, using the built-in programming language Ansys Parametric Design Language (APDL), a proprietary module enabling direct determination of GSIFs was implemented. This module also takes into account the previously mentioned criterion of selecting nodes. The following types of specimens were modelled in numerical simulations:
  • a rectangular plate with a single edge sharp corner under uniaxial tension (Figure 3a);
  • a rectangular plate with a double edge sharp corner under uniaxial/biaxial tension (Figure 3b);
  • a rectangular plate with a central sharp corner under uniaxial/biaxial tension (Figure 4a);
  • a rectangular plate with a central sharp corner under pure shear loading (Figure 4b).
As the tested samples have one or two planes of symmetry, only their halves or quarters were modelled (shaded area in the figures above). Symmetry (for tension samples—Figure 3 and Figure 4a) and anti-symmetry (for shear samples—Figure 4b) boundary conditions were assumed in the symmetry planes.
The samples were discretized using quadrangular, eight-node finite elements (Figure 5). The mesh of division into finite elements was densified in the tip region. Furthermore, the tip of the notch was surrounded by a special triangular finite element with a shape function, which could simulate the singularities of displacements of the r0.5 type. In this way, a better representation of the singular stress fields occurring in the vicinity of the corner tip was obtained.
The applied load σy1 was constant and equal to 1 Pa. However, the value of the load σy2 was determined on the basis of Formula (15), resulting from the condition of continuity of strain −εy1 = εy2:
σ y 2 = E 2 ( σ y 1 σ x ν 1 ) E 1 + σ x ν 2 .  
For tensile specimens, numerical calculations were performed for variable proportions of longitudinal and transverse loads −σx/σy1. Moreover, the simulations were prepared for various relative stiffness of the individual components of the bi-material– Γ = μ 1 / μ 2   , where μ i = E i 2 ( 1 + ν i ) . In all simulations, it was assumed that the Young’s modulus of material 1 was constant and be equal to: E1 = 1 × 109 Pa. However, the Young’s modulus of material 2 –E2- was variable and depended on the parameter Γ. The Poisson coefficients also depended on this parameter, which were respectively: ν1 = ν2 = 0.3 (for Γ = 1), ν1 = 0.3, ν2 = 0.35 (for Γ > 1), and ν1 = 0.35, ν2 = 0.3 (for Γ < 1).
Numerical tests were performed for specimens with various notch angles ψ and proportions of characteristic dimensions a/w. Moreover, it was arbitrarily assumed that the notch height a is constant and amounts to 1 m. The ratio of the height and width of the specimens was also kept constant—h/w = 2. The plane stress condition was assumed in all simulations.

3. Results and Discussion

The quantitative description of singular stress fields is obtained by finding generalised stress intensity factors. For their determination, using the previously described method, the numerical data (hoop and tangential stresses at the interface of the bi-material) and eigenvalues λj are necessary. The latter, for the samples described in Section 2.2, were calculated on the basis of the dependence (7) and (8) and are presented in Table 1 and Table 2.
GSIFs units—[Pa m1−λj]—depend on the geometrical and material parameters of the specimens. This prevents direct comparison of the results obtained. Such inconvenience can be eliminated by normalising GSIFs [20]. In the presented work, the determined values of the Kj coefficients were normalised using the following formula:
F j = K j / ( σ y 1 π a 1 λ j ) ,   ( j = I , I I )
The results obtained for each type of specimens are presented below.

3.1. Rectangular Plate with a Single Edge Sharp Corner under Uniaxial Tension

Table 3 shows the normalised FI factors calculated for a tensile bi–material with a single edge notch. The calculations were made for various relative stiffness and notch apex angles.
By analysing the results presented in Table 3, it can be seen that the normalised FI values increase with
  • an increase in the notch angle 2ψ;
  • a decrease in the relative stiffness Γ (this tendency is consistent with the distribution of the normalised stress intensity factors determined for the crack initiated at the interface [16,18]).
It was also found that as the height of the notch increased, regardless of the material parameters and its opening angle, the FI coefficients increased (Figure 6).

3.2. Rectangular Plate with a Double Edge Sharp Corner under Uniaxial/Biaxial Tension

Similar tests, as for the case described in Section 3.1, were performed for the problem of a double edge-notched plate. The obtained results are presented in Table 4 and Figure 7 and Figure 8.
The stress intensity coefficients were calculated for two load variants: uniaxial (σx = 0) and biaxial (σx/σy1 ≠ 0) tension. The influence of the tested parameters—relative stiffness Γ, vertex angle 2ψ, and relative notch height a/w—on the values of normalised stress intensity factors was identical as for the single edge-notched plate problem.
As for the influence of the ratio of loads perpendicular and parallel to the interface (σx/σy1) on the values of the FI coefficients, on the basis of the analyses performed (Table 4, Figure 8) it can be concluded that with increasing load σx:
  • the FI value significantly decreases for Γ ≤ 1;
  • increases slightly for the case where Γ > 1.
Of course, for the case of a crack in a homogeneous material (Γ = 1, 2ψ = 0°), the application of an additional load σx does not affect the stress intensity factor.

3.3. Rectangular Plate with a Central Sharp Corner under Uniaxial/Biaxial Tension

The same tests, as for the case described in Section 3.2, were performed for the problem with the centre-notched plate. The obtained results are presented in Table 5 and Figure 9 and Figure 10.
The performed analyses showed that the factors Fj always increase with a decrease in the relative stiffness Γ and an increase in the height of the notch and its apex angle.
For samples subjected to the biaxial tension load, similarly to the problem of the double-sided notch, an increase in the load σx causes either a decrease in the values of the FI (for Γ ≤ 1) coefficients or their increase (for Γ > 1).

3.4. Rectangular Plate with a Central Sharp Corner under Pure Shear Loading

The test results for the notched samples subjected to tangential load (Mode II) are shown in Table 6 and Figure 11.
As for elements subjected to shear loads, the values of the normalised factors FII (Table 6) change in the same way as in the case of applying tensile loads (increase with the decrease of the relative stiffness and the increase of the notch tip angle).
In Figure 11, for an element with a central notch subjected to pure shear (for arbitrarily selected geometric and material parameters), the stress distributions obtained from the analytical solution were compared with the results obtained through FEM. Good agreement of both solutions was obtained in the apex region of about 10% of the notch height. Similar comparisons were also made for the specimens discussed in Section 3.1, Section 3.2 and Section 3.3 (not included in the paper). In each case, a similar compliance of the analytical description with the FEM solution was obtained.

4. Conclusions

The paper presents analytical and numerical solutions for the plane problem of a sharp corner, assuming that its tip touches the line separating two different elastic materials.
Two variants of loading were considered—tensile (transverse and/or longitudinal to the interface) and shear in the plane of symmetry of the notch. For both load cases, eigenequations were determined, which depended on the mechanical properties of both components of the composite and the notch tip angle. It was found that for each loading variant, there is always one singular term described by the real eigenvalues λj. Moreover, it was observed that the stress singularity strength decreased with increasing notch apex angle and relative stiffness Γ.
Furthermore, formulas were derived to determine the individual components of the stress tensor occurring in the notch tip area with the use of generalised stress intensity factors Kj. The Kj coefficients were also calculated for three variants of the notch location: a single edge sharp notch, a double edge sharp, and a central sharp notch. The calculations were made for various notch apex angles as a function of material constants.
It was found that the normalised values of generalised stress intensity factors increase with a decrease in the relative stiffness Γ and an increase in the notch height and its apex angle. Moreover, it was found that in the case of specimens subjected to biaxial tension, an increase in the load perpendicular to the interface (σx) causes either a decrease in the values of the normalised stress intensity factors (Γ ≤ 1) or their increase (Γ ≤ 1).
The derived analytical solution describing the individual stress components was compared with the results obtained using the FEM. Both solutions were found to be very compatible in the apex region of about 10% of the notch height.
The research results presented in this article can be used by other researchers in many ways, e.g., as comparative data. The analytical and numerical description of singular stress fields can also be used to develop a fracture criterion of structural elements with this type of material defect. The development of such a criterion will be the aim of the future works of the author.

Author Contributions

G.M., supervised all the analyses, developed the analytical solutions, and wrote the paper; D.S. and A.B., carried out the FEM analyses; M.M.A., A.E. and M.S., visualization, validation, and formal analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This publication was financed through the program of the Ministry of Science and Higher Education of Poland named “Regional Initiative of Excellence” in 2019–2022 project number 011/RID/2018/19.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. At the time the project was carried out, there was no obligation to make the data publicly available.

Acknowledgments

This research was partially financed through subsidy of the Ministry of Science and Higher Education of Poland for the discipline of mechanical engineering at the Faculty of Mechanical Engineering Bialystok University of Technology WZ/WM-IIM/4/2020.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Griffits, A.A. The phenomena of rupture and flow in solids. Philos. Trans. R. Soc. London Ser. A Contain. Pap. A Math. Phys. Character 1921, 221, 163–198. [Google Scholar]
  2. Sih, G.C. Strain-energy-density factor applied to mixed mode crack problems. Int. J. Fract. 1974, 10, 305–321. [Google Scholar] [CrossRef]
  3. McClintock, F.A. Ductile fracture instability in shear. J. Appl. Mech. 1958, 25, 582–588. [Google Scholar] [CrossRef]
  4. Yosibash, Z.; Priel, E.; Leguillon, D. A failure criterion for brittle elastic materials under mixed-mode loading. Int. J. Fract. 2006, 141, 291–312. [Google Scholar] [CrossRef]
  5. Sih, G.C.; Chen, E.P.; Sih, G.C.; Chen, E.P. Cracks in materials possessing homogeneous anisotropy. In Cracks in Composite Materials; Springer: Dordrecht, The Netherlands, 1981; pp. 1–101. [Google Scholar]
  6. Sun, C.T.; Jih, C.J. On strain energy release rates for interfacial cracks in bi–material media. Eng. Fract. Mech. 1987, 28, 13–20. [Google Scholar] [CrossRef]
  7. Krishnan, A.; Roy Xu, L. An experimental study on the crack initiation from notches connected to interfaces of bonded bi-materials. Eng. Fract. Mech. 2013, 111, 65–76. [Google Scholar] [CrossRef]
  8. Mieczkowski, G. Stress fields and fracture prediction for an adhesively bonded bimaterial structure with a sharp notch located on the interface. Mech. Compos. Mater. 2017, 53, 305–320. [Google Scholar] [CrossRef]
  9. Ballarini, R. A rigid line inclusion at a bimaterial interface. Eng. Fract. Mech. 1990, 37, 1–5. [Google Scholar] [CrossRef]
  10. Mieczkowski, G. Description of stress fields and displacements at the tip of a rigid, flat inclusion located at interface using modified stress intensity factors. Mechanika 2015, 21, 91–98. [Google Scholar] [CrossRef] [Green Version]
  11. Mieczkowski, G. Stress fields at the tip of a sharp inclusion on the interface of a bimaterial. Mech. Compos. Mater. 2016, 52, 601–610. [Google Scholar] [CrossRef]
  12. Carpinteri, A.; Paggi, M.; Pugno, N. Numerical evaluation of generalised stress-intensity factors in multi-layered composites. Int. J. Solids Struct. 2006, 43, 627–641. [Google Scholar] [CrossRef]
  13. Sisodia, S.M.; Bull, D.J.; George, A.R.; Gamstedt, E.K.; Mavrogordato, M.N.; Fullwood, D.T.; Spearing, S.M. The effects of voids in quasi-static indentation of resin-infused reinforced polymers. J. Compos. Mater. 2019, 53, 4399–4410. [Google Scholar] [CrossRef]
  14. Mehdikhani, M.; Gorbatikh, L.; Verpoest, I.; Lomov, S.V. Voids in fiber-reinforced polymer composites: A review on their formation, characteristics, and effects on mechanical performance. J. Compos. Mater. 2019, 53, 1579–1669. [Google Scholar] [CrossRef]
  15. Zak, A.R.; Williams, M.L. Crack point stress singularities at a bi–material interface. J. Appl. Mech. Trans. ASME 1960, 30, 142–143. [Google Scholar] [CrossRef]
  16. Cook, T.S.; Erdogan, F. Stresses in bonded materials with a crack perpendicular to the interface. Int. J. Eng. Sci. 1972, 10, 677–697. [Google Scholar] [CrossRef]
  17. Wang, W.C.; Chen, J.T. Theoretical and experimental re-examination of a crack perpendicular to and terminating at the bimaterial interface. J. Strain Anal. Eng. Des. 1993, 28, 53–61. [Google Scholar] [CrossRef]
  18. Lin, K.Y.; Mar, J.W. Finite element analysis of stress intensity factors for cracks at a bi–material interface. Int. J. Fract. 1976, 12, 521–531. [Google Scholar] [CrossRef]
  19. Meguid, S.A.; Tan, M.; Zhu, Z.H. Analysis of cracks perpendicular to bimaterial interfaces using a novel finite element. Int. J. Fract. 1985, 73, 1–23. [Google Scholar] [CrossRef]
  20. Dai-Heng, C. A crack normal to and terminating at a bimaterial interface. Eng. Fract. Mech. 1994, 49, 517–532. [Google Scholar] [CrossRef]
  21. Keikhaie, M.; Keikhaie, N.; Keikhaie, R.; Kaykha, M.M. Stress Intensity Factors in Two Bonded Elastic Layers Containing Crack Perpendicular on the Interface with Different Elastic Properties. J. Mod. Phys. 2015, 6, 640–647. [Google Scholar] [CrossRef] [Green Version]
  22. Bogy, D.B. On the plane elastostatic problem of a loaded crack terminating at a material interface. J. Appl. Mech. Trans. ASME 1971, 38, 911–918. [Google Scholar] [CrossRef]
  23. Chang, J.; Xu, J.Q. The singular stress field and stress intensity factors of a crack terminating at a bimaterial interface. Int. J. Mech. Sci. 2007, 49, 888–897. [Google Scholar] [CrossRef]
  24. Náhlík, L.; Knésl, Z.; Klusák, J. Crack initiation criteria for singular stress concentrations Part III: An Application to a Crack Touching a Bimaterial Interface. Eng. Mech. 2008, 15, 99–114. [Google Scholar]
  25. Selvarathinam, A.S.; Weitsman, Y.J. Fracture in angle-ply ceramic matrix composites. Int. J. Fract. 2000, 102, 71–84. [Google Scholar] [CrossRef]
  26. Li, Y.; Song, M. Method to calculate stress intensity factor of V-notch in bi-materials. Acta Mech. Solida Sin. 2008, 21, 337–346. [Google Scholar] [CrossRef]
  27. Williams, M.L. Stress Singularities Resulting from Various Boundary Conditions in Angular Corners of Plates in Extension. J. Appl. Mech. 1952, 19, 526–528. [Google Scholar] [CrossRef]
  28. Parton, V.Z.; Perlin, P.I. Mathematical Methods of the Theory of Elasticity; Mir Publishers: Moscow, Russia, 1984. [Google Scholar]
  29. Hein, V.L.; Erdogan, F. Stress singularities in a two-material wedge. Int. J. Fract. Mech. 1971, 7, 317–330. [Google Scholar] [CrossRef]
  30. Mieczkowski, G. Determination of stress intensity factors for elements with sharp corner located on the interface of a bi–material structure or homogeneous material. Acta. Mech. 2021, 232, 709–724. [Google Scholar] [CrossRef]
  31. Borawski, A.; Szpica, D.; Mieczkowski, G. Verification tests of frictional heat modelling results. Mechanika 2020, 26, 260–264. [Google Scholar] [CrossRef]
  32. Łukaszewicz, A. Nonlinear Numerical Model of Friction Heating during Rotary Friction Welding. J. Frict. Wear 2018, 39, 476–482. [Google Scholar] [CrossRef]
  33. Borawski, A.; Szpica, D.; Mieczkowski, G.; Borawska, E.; Awad, M.M.; Shalaby, R.M.; Sallah, M. Theoretical Analysis of the Motorcycle Front Brake Heating Process during High Initial Speed Emergency Braking. J. Appl. Comput. Mech. 2020, 6, 1431–1437. [Google Scholar]
  34. Yevtushenko, A.A.; Kuciej, M.; Grzes, P.; Wasilewski, P. Temperature in the railway disc brake at a repetitive short-term mode of braking. Int. Commun. Heat Mass Transf. 2017, 84, 102–109. [Google Scholar] [CrossRef]
  35. Mieczkowski, G.; Borawski, A.; Szpica, D. Static electromechanical characteristic of a three-layer circular piezoelectric transducer. Sensors 2020, 20, 222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Liu, Z.; Chen, J.; Zou, X. Modeling the piezoelectric cantilever resonator with different width layers. Sensors 2021, 21, 87. [Google Scholar] [CrossRef] [PubMed]
  37. Asadi Dereshgi, H.; Dal, H.; Yildiz, M.Z. Piezoelectric micropumps: State of the art review. Microsyst. Technol. 2021, 1–29. [Google Scholar] [CrossRef]
  38. Treifi, M.; Oyadiji, S.O. Bi–material V-notch stress intensity factors by the fractal-like finite element method. Eng. Fract. Mech. 2013, 105, 221–237. [Google Scholar] [CrossRef]
  39. Dunn, M.L.; Suwito, W.; Cunningham, S. Stress intensities at notch singularities. Eng. Fract. Mech. 1997, 57, 417–430. [Google Scholar] [CrossRef]
Figure 1. A notch terminating at a bi–material interface.
Figure 1. A notch terminating at a bi–material interface.
Materials 14 04466 g001
Figure 2. A solution of eigenequations (ν1 = ν2 = 0.3, plane stress condition), (a) for Mode I (7), (b) for Mode II (8).
Figure 2. A solution of eigenequations (ν1 = ν2 = 0.3, plane stress condition), (a) for Mode I (7), (b) for Mode II (8).
Materials 14 04466 g002
Figure 3. Geometry and method of fixing and loading specimens: (a) with a single edge-notched plate, and (b) with a double edge-notched plate.
Figure 3. Geometry and method of fixing and loading specimens: (a) with a single edge-notched plate, and (b) with a double edge-notched plate.
Materials 14 04466 g003
Figure 4. Geometry method of fixing and loading of specimens with a centre-notched plate subjected to (a) uniaxial/biaxial tension and (b) pure shear loading.
Figure 4. Geometry method of fixing and loading of specimens with a centre-notched plate subjected to (a) uniaxial/biaxial tension and (b) pure shear loading.
Materials 14 04466 g004
Figure 5. A typical finite element mesh used for modelling the specimens.
Figure 5. A typical finite element mesh used for modelling the specimens.
Materials 14 04466 g005
Figure 6. Variation of the normalised stress intensity factors FI with the relative stiffness Γ and notch height a/w, (a) notched bi-material with apex angle 2ψ = 30°, (b) notched bi-material with apex angle 2ψ = 60°.
Figure 6. Variation of the normalised stress intensity factors FI with the relative stiffness Γ and notch height a/w, (a) notched bi-material with apex angle 2ψ = 30°, (b) notched bi-material with apex angle 2ψ = 60°.
Materials 14 04466 g006
Figure 7. Variation of the normalised stress intensity factors FI with the relative stiffness Γ and notch height a/w, (a) notched bi-material with 2ψ = 0° apex angle, (b) notched bi-material with 2ψ = 90° apex angle.
Figure 7. Variation of the normalised stress intensity factors FI with the relative stiffness Γ and notch height a/w, (a) notched bi-material with 2ψ = 0° apex angle, (b) notched bi-material with 2ψ = 90° apex angle.
Materials 14 04466 g007
Figure 8. Variation of the normalised stress intensity factors FI with the relative stiffness Γ and the notch angle 2ψ, (a) σx/σy1 = 1, (b) σx/σy1 = 0.5, a/w = 0.4.
Figure 8. Variation of the normalised stress intensity factors FI with the relative stiffness Γ and the notch angle 2ψ, (a) σx/σy1 = 1, (b) σx/σy1 = 0.5, a/w = 0.4.
Materials 14 04466 g008
Figure 9. Variation of the normalised stress intensity factors FI with the relative stiffness Γ and notch height a/w, (a) notched bi-material with apex angle 2ψ = 30°, (b) notched bi-material with apex angle 2ψ = 90°.
Figure 9. Variation of the normalised stress intensity factors FI with the relative stiffness Γ and notch height a/w, (a) notched bi-material with apex angle 2ψ = 30°, (b) notched bi-material with apex angle 2ψ = 90°.
Materials 14 04466 g009
Figure 10. Variation of the normalised stress intensity factors FI with the relative stiffness Γ and the notch angle 2ψ, (a) for σx/σy1 = 1, (b) for σx/σy1 = 0.5, a/w = 0.4.
Figure 10. Variation of the normalised stress intensity factors FI with the relative stiffness Γ and the notch angle 2ψ, (a) for σx/σy1 = 1, (b) for σx/σy1 = 0.5, a/w = 0.4.
Materials 14 04466 g010
Figure 11. Comparison of selected stress distributions (pure shear loading) obtained by the FEM (points) and according to the analytical Formula (11) using the factors from Table 6, (continuous lines), (a) hoop stress, (b) tangential stress, (φ = 0°, 2ψ = 30°, Γ = 0.5).
Figure 11. Comparison of selected stress distributions (pure shear loading) obtained by the FEM (points) and according to the analytical Formula (11) using the factors from Table 6, (continuous lines), (a) hoop stress, (b) tangential stress, (φ = 0°, 2ψ = 30°, Γ = 0.5).
Materials 14 04466 g011aMaterials 14 04466 g011b
Table 1. The eigenvalues λI (Mode I).
Table 1. The eigenvalues λI (Mode I).
2ψ
[°]
γ
[°]
λI
Γ = 0.1Γ = 0.2Γ = 0.5Γ = 1Γ = 2Γ = 5Γ = 10
0900.681450.640750.56383/0.56383 *0.50.42944/0.42944 *0.325790.25150
30750.729360.671110.573860.501450.430420.335140.26647
60600.795360.720260.599750.512220.431660.336270.27426
90450.866120.786940.649840.544480.446390.336060.27385
ν1 = 0.3, ν2 = 0.35 for Γ > 1, ν1 = 0.35, ν2 = 0.3 for Γ < 1; *—from reference [18].
Table 2. The eigenvalues λII (Mode II).
Table 2. The eigenvalues λII (Mode II).
2ψ
[°]
γ
[°]
λII
Γ = 0.1Γ = 0.2Γ = 0.5Γ = 1Γ = 2Γ = 5Γ = 10
0900.681450.640750.563830.50.429440.325790.25150
30750.778210.737730.660740.598190.528280.418050.33128
60600.927760.881730.796160.730900.662010.552260.45604
9045110.991050.908530.832060.730520.64428
ν1 = 0.3, ν2 = 0.35 for Γ > 1, ν1 = 0.35, ν2 = 0.3 for Γ < 1.
Table 3. Values of normalised stress intensity factors FI calculated for a rectangular plate with a single edge sharp corner under uniaxial tension, a/w = 0.2.
Table 3. Values of normalised stress intensity factors FI calculated for a rectangular plate with a single edge sharp corner under uniaxial tension, a/w = 0.2.
ΓFI
2ψ = 0°2ψ = 30°2ψ = 60°2ψ = 90°
0.14.3176.0708.76511.990
0.23.0763.7224.7135.991
0.51.9082.0192.2362.609
11.3611.3751.431/2.220 **1.579/2.471 **
1.367 * 2.230 *2.478 *
20.9830.9830.9891.042
50.6360.6590.6590.659
100.4520.4900.5020.504
*—From reference [38]; **—calculated for a/w = 0.4.
Table 4. Values of normalised stress intensity factors FI calculated for the rectangular plate with a double edge sharp corner under uniaxial/biaxial tension, a/w = 0.4.
Table 4. Values of normalised stress intensity factors FI calculated for the rectangular plate with a double edge sharp corner under uniaxial/biaxial tension, a/w = 0.4.
ΓFI
2ψ = 0°2ψ = 30°2ψ = 60°2ψ = 90°
σx = 0σxy1 = 2σx = 0σx/σy1 = 2σx = 0σx/σy1 = 2σx = 0σx/σy1 = 2
0.14.3323.4885.8473.688.4694.08611.4984.600
0.22.9502.4303.4893.004.4122.4245.5802.578
0.51.6851.5251.7671.401.9581.3182.2811.343
11.1321.1321.1401.031.1890.9501.3240.969
20.7720.8060.7770.7780.7810.7830.8210.832
50.4800.4930.4970.510.4960.5090.5170.518
100.3440.3460.3680.380.3760.3950.3810.406
Table 5. Values of normalised stress intensity factors FI calculated for the rectangular plate with a central sharp corner under uniaxial/biaxial tension, a/w = 0.4.
Table 5. Values of normalised stress intensity factors FI calculated for the rectangular plate with a central sharp corner under uniaxial/biaxial tension, a/w = 0.4.
ΓFI
2ψ = 0°2ψ = 30°2ψ = 60°2ψ = 90°
σx = 0σx/σy1 = 2σx = 0σx/σy1 = 2σx = 0σx/σy1 = 2σx = 0σx/σy1 = 2
0.14.1463.3955.8183.7668.7403.79312.362.415
0.22.7422.4103.4972.4604.6002.3006.121.270
0.51.6201.5001.7671.4642.0671.3042.560.710
11.109
1.004 *
1.1091.146
1.027 *
1.0761.254
1.112 *
0.9531.47
1.263 *
0.596
0.996 **1.028 **1.115 **1.267 **
20.7710.8030.7820.7970.8180.8210.910.512
50.4820.5050.4980.5260.5020.5080.530.406
100.3450.3530.3660.3810.3790.3850.390.337
**—From reference [39], *—calculated for a/w = 0.1.
Table 6. Values of normalised stress intensity factors FII calculated for a rectangular plate with a central sharp corner under pure shear loading, a/w = 0.4.
Table 6. Values of normalised stress intensity factors FII calculated for a rectangular plate with a central sharp corner under pure shear loading, a/w = 0.4.
ΓFII
2ψ = 0°2ψ = 30°2ψ = 60°2ψ = 90°
0.113.73916.14619.152-
0.26.5017.7969.387-
0.52.2792.8453.5304.669
11.0341.3451.7312.124
20.4510.6130.8361.074
50.1390.1970.2930.425
100.0620.0880.1350.218
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mieczkowski, G.; Szpica, D.; Borawski, A.; Awad, M.M.; Elgarayhi, A.; Sallah, M. Investigation of the Near-Tip Stress Field of a Notch Terminating at a Bi-Material Interface. Materials 2021, 14, 4466. https://doi.org/10.3390/ma14164466

AMA Style

Mieczkowski G, Szpica D, Borawski A, Awad MM, Elgarayhi A, Sallah M. Investigation of the Near-Tip Stress Field of a Notch Terminating at a Bi-Material Interface. Materials. 2021; 14(16):4466. https://doi.org/10.3390/ma14164466

Chicago/Turabian Style

Mieczkowski, Grzegorz, Dariusz Szpica, Andrzej Borawski, Mohamed M. Awad, Ahmed Elgarayhi, and Mohammed Sallah. 2021. "Investigation of the Near-Tip Stress Field of a Notch Terminating at a Bi-Material Interface" Materials 14, no. 16: 4466. https://doi.org/10.3390/ma14164466

APA Style

Mieczkowski, G., Szpica, D., Borawski, A., Awad, M. M., Elgarayhi, A., & Sallah, M. (2021). Investigation of the Near-Tip Stress Field of a Notch Terminating at a Bi-Material Interface. Materials, 14(16), 4466. https://doi.org/10.3390/ma14164466

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop