Next Article in Journal
Synthesis of Sodium Cobalt Fluoride/Reduced Graphene Oxide (NaCoF3/rGO) Nanocomposites and Investigation of Their Electrochemical Properties as Cathodes for Li-Ion Batteries
Next Article in Special Issue
Influence of the Undercut Anchor Head Angle on the Propagation of the Failure Zone of the Rock Medium
Previous Article in Journal
The Bond-Slip Relationship at FRP-to-Brick Interfaces under Dynamic Loading
Previous Article in Special Issue
Load Controlled Fatigue Behaviour of Microplasma Arc Welded Thin Titanium Grade 5 (6Al-4V) Sheets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stress Concentration Factors for Welded Plate T-Joints Subjected to Tensile, Bending and Shearing Loads

by
Krzysztof L. Molski
and
Piotr Tarasiuk
*
Faculty of Mechanical Engineering, Bialystok University of Technology, Wiejska 45C, 15-351 Bialystok, Poland
*
Author to whom correspondence should be addressed.
Materials 2021, 14(3), 546; https://doi.org/10.3390/ma14030546
Submission received: 14 December 2020 / Revised: 20 January 2021 / Accepted: 21 January 2021 / Published: 24 January 2021
(This article belongs to the Special Issue Fatigue and Fracture of Materials)

Abstract

:
The paper deals with the problem of stress concentration at the weld toe of a plate T-joint subjected to axial, bending, and shearing loading modes. Theoretical stress concentration factors were obtained from numerical simulations using the finite element method for several thousand geometrical cases, where five of the most important geometrical parameters of the joint were considered to be independent variables. For each loading mode—axial, bending, and shearing—highly accurate closed form parametric expression has been derived with a maximum percentage error lower than 2% with respect to the numerical values. Validity of each approximating formula covers the range of dimensional proportions of welded plate T-joints used in engineering applications. Two limiting cases are also included in the solutions—when the weld toe radius tends to zero and the main plate thickness becomes infinite.

Graphical Abstract

1. Introduction

Welded T-connections are commonly used engineering structures. Their significant advantage lies in low costs of production due to the use of simple elements, like tubes of circular or rectangular cross sections and plates, usually connected by welding. There are many types of such joints produced in various configurations: as tubular, tubular-plate, and plate couples. The fatigue failure of welded structural components subjected to cyclic loads is the most frequent cause of damage, which begins at critical zones where stress concentration occurs [1]. Various standards and recommendations presented in References [2,3,4,5,6,7] are very helpful in designing welded structures.
Stress concentration factor (SCF) values determined for the known geometrical and loading conditions are of great importance in fatigue design because they enable relating the loading history to the critical zone of the structure. Some examples of fatigue analysis in tubular joints applied to various types of structures and devices are given in References [8,9,10,11,12] where stress concentration plays a key role. Fatigue damage analysis considering crack initiation and propagation processes in welded T-joints was published in References [13,14,15]. Numerous formulas of stress concentration factors for various types of welded T-joints, based on the finite and boundary element analyses, were presented in References [16,17,18,19,20].
Some important phenomena associated with welding processes consist of a random character of key geometrical parameters influencing the SCF value. In such cases, appropriate statistic data obtained directly from real structures and regarding changed values and shapes of weldments are necessary. Examples of such data in the form of histograms can be found in the literature [21,22,23,24,25]. Various types of the weld face geometry—plane, convex, or concave—may be obtained during manufacturing, which generally depend on the welding method, position, and parameters of the process, e.g., welding speed, current, etc. Some examples of real weld shapes and profiles influencing SCF values in cases of as-welded and intentionally improved weld toe zones are also given in References [25,26,27,28,29]. Many techniques have also been developed to improve the endurance limit and reliability of the structure (grinding, milling, re-melting techniques, etc.) [30]. Each process produces specific radii influencing local stress concentration.
The application of the known parametric formulas in fatigue design encounters some difficulties arising from their accuracy, range of validity, and different ways of defining SCF. On the other hand, assessment of fatigue life requires high accuracy of SCF solutions and should cover the whole range of values of all key parameters influencing SCFs.
The present study deals with determination of stress concentration factors for non-load carrying welded plate T-joints subjected to axial, bending and shearing loading modes. The main objective was to develop highly accurate SCF formulas of a wide ranging validity, including two limiting cases: when the toe radius ρ tends to zero and when the main plate thickness t becomes infinite.

2. Methodology

2.1. General Assumptions

Five key geometrical parameters—ρ, a, θ, t, and T—have been selected and considered as independent variables, where t and T represent the main and attachment plate thickness, respectively, ρ is the weld toe radius, θ indicates the weld angle, and a is the weld throat thickness defined here as the shortest distance between the joint root and the weld face. Leg length parameters, h and hp, characterizing the weld size, are also used in the literature, e.g., in References [16,17,18]. Shape and geometrical dimensions of a plate T-joint, as well as three loading modes considered in the present work, are depicted schematically in Figure 1. Small circles drawn symmetrically at the weld toe zone indicate locations of the maximum stresses.
The present analysis considered the following assumptions.
  • Joint material is linear elastic, isotropic, and homogeneous.
  • Elastic properties of the weld material and the base material, represented by elastic constants, are the same.
  • The welded joint is free from residual stresses, structural irregularities, and imperfections, including lack of penetration defects.
  • Both welds are symmetrical with respect to the axis of the attachment plate, which is perpendicular to the main plate.
  • Weld faces are plane and contour of the weldment is smooth with a transition toe radius ρ > 0.
  • External load – axial, bending, and shearing – is applied at the opposite sides of the main plate sufficiently far from the welds.
  • Small deformations occur in the whole body.
  • Five geometrical parameters: ρ, a, θ, t, and T are considered as independent variables varying in the following ranges: 0 < ρ/a ≤ 1.3, 0 < a/t ≤ 1.3, 1 ≤ T/a ≤ 4, and 30° ≤ θ ≤ 60°.
  • Stress concentration factors for axial, bending, and shearing loads are defined as Ktt = σ1maxt, Ktb = σ1max/σb, and Kts = τmaxs, respectively.
Theoretical stress concentration factor SCF is valid for linear-elastic, homogeneous, and isotropic materials as well as bodies of defined geometry, loading and displacement boundary conditions. Any departure from these assumptions, as material anisotropy or geometrical irregularities (flaws, imperfections, etc.) may significantly change the stress distribution in the body, including the maximum stress value at the notch root.
The aim of the work was to obtain SCF approximating formulae covering possible proportions between geometrical parameters of the T-joint used in engineering practice. For this reason, the supposed margin for the normalized parameters was very wide. The weld toe radius ρ 1.3 times larger than the weld throat may be considered useless for practical applications, but such a limiting value can guarantee that all smaller radii fall into the range of application. The second case when a/t→0 may also take place in engineering applications while a relatively thin plate is welded to a very thick main plate.

2.2. General Approach

The procedure of reaching main objectives of the present study consisted of several steps. First, conveniently defined new geometrical parameters had to be specified. Due to the fact that SCF values do not vary when all dimensions of the body proportionally change, the following normalized non-dimensional quantities were introduced.
X = ρ / ( ρ + a ) = ρ / a ρ / a + 1
Y = a / ( a + t ) = a / t a / t + 1
Z = T / a
In this way, four independent variables X, Y, Z, and θ, containing the five key dimensions of the joint, were taken into consideration. The upper weld toe radius, ρ1, shown in Figure 1, was chosen as equal to 0.1a and had no influence on stress concentration factors. The second step consisted of numerical modelling of the joint using the finite element method, where mutual proportions between the basic joint dimensions systematically changed. In this way, several thousands of SCF values, covering the assumed range of validity, were determined. More detailed information on the numerical modelling and SCF calculating procedures are presented in the next section. The third step of the approach consisted of choosing a general form of the mathematical representation of SCFs’ approximating functions. An expression applied in the present study was similar to that presented by the authors in Reference [31], where a cruciform joint of θ = 45°, subjected to tensile and bending loads, was analysed. However, the problem considered in the present work was more general than the previous one because the weld angle θ varied in the range of 30–60°, increasing the number of variables in the parametric formula. In the numerical Finite Element Method (FEM) modelling, the weld angle θ was changed by 2.5° and the whole calculating procedure, carried out for a given angle, had to be repeated 12 times. Hence, a general form of the approximating function is given by Equation (4).
K t = X n P ( X , Y , θ , Z 0 ) κ ( X , Y , Z , θ , Z 0 ) ,
where a singular term Xn accounts for the stress concentration effects when ρ→0. The function P(X,Y,θ,Z0), represented by polynomials, can be derived from numerical SCF solutions normalized with respect to the singular term, Kt/Xn for different X, Y, and θ, while Z = Z0. Equation (4) indicates that the P function cannot directly relate SCF to Z because it is derived for a particular, arbitrarily chosen value of Z = Z0. For this reason, the additional correction function is κ(X,Y,Z,θ,Z0), satisfying Equation (5),
κ ( X , Y , Z = Z 0 , θ ) = 1
was necessary to account for the relative stiffener thickness on SCF. In other words, the function κ can be interpreted as a multiplier, which enables us to calculate SCF for any value of Z with respect to the one corresponding to Z0, for the known parameters X, Y, and θ.
The exponent n, denoted as ns for shearing load, is necessary to describe stress concentration effects when the weld toe radius tends to zero. In References [31,32], it was proven that the correct n values are the same as the stress field singularity exponents corresponding to the characteristic angle of the sharp corner [33], into which the weld toe region is converted while the weld toe radius ρ tends to zero. According to Reference [32], numerical values of the exponent n can be calculated for axial and bending load, using Equation (6),
n = 0.63662 θ 0.09330 θ 2 1 + 0.77635 θ + 0.04075 θ 1.5 0.00499 θ 2 + 0.13365 θ 2.5
which is valid in the range of 0 ≤ θ ≤ π/2 with an accuracy of five significant digits.
For anti-plane deformation produced by the shearing load, the exact value of the exponent ns equals:
n s = θ θ + π
where θ is in radians.
The next step of the procedure consisted of approximating the P(X,Y,Z0,θ) function to the normalized numerical SCF data using the least squares method. After particular constants (coefficients and exponents) of the P function were determined, the approximation accuracy was verified. In this way, one part of the solution, where Z = Z0, was completed. The final step of the present approach consisted of finding unknown coefficients and exponents of the correction function κ(X,Y,Z,θ,Z0) accounting for the influence of Z values on SCFs. After performing the second validation of the full solution including κ, the close form approximating functions for calculating SCFs were determined. The whole procedure was consecutively repeated for each loading mode.

3. Numerical FEM Modelling and Some SCFs’ Results

3.1. Tensile and Bending Load

Systematic numerical modelling of plate T-joints was carried out by using the Finite Element Method (FEM) incorporated in the ANSYS 19 Multiphysics program. The PLANE182 finite elements that were used were defined by four nodes having 2 degrees of freedom at each node. The elastic material constants were assumed to be similar to those for steel, where Young’s modulus E = 210 GPa and Poisson’s ratio ν = 0.3. However, according to the theory of elasticity, stress field in the body and SCF values do not depend on the elastic constants of the material in all cases modeled. Loading and displacement boundary conditions imposed on one half of the body, as well as its shape, are shown in Figure 2. According to the principle of de Saint Venant, the minimum length of the main plate measured from the weld toe was about 4.5t. This value was confirmed by conducting preliminary numerical tests.
Approximately 820,000 finite elements were used for each geometrical case of the modelled T-joint. Special attention was paid to the finite element mesh density at the weld toe zone, where high values of the stress gradients occur. For this reason, the dimensional ratio of neighbouring finite elements was about 1.2 and approximately 40 to 70 finite elements were used along the weld toe circular arc, as described by the radius ρ. One example of a finite element mesh is shown in Figure 3.
An example of the solution obtained for the tensile load for the nominal uniform stress σt = 1 MPa applied over the right end of the body is shown in Figure 4. The first principal stress results offer the most convenient representation because the particular value of the SCF is directly accessible. Therefore, the nominal stress in all numerical models was always equal to one.
In the presented case, the maximum principal stress σ1max equals 1.56459 MPa, which is directly seen on the scale below the picture. This value is also interpreted as the stress concentration factor Ktt.

3.2. Shearing Load

The formulation of two-dimensional problems for the anti-plane state of deformation is different than for in-plane loading modes. In the case of the anti-plane problem, only one displacement coordinate W in each point is sufficient to describe deformation of the body and W(x,y) may be treated as a scalar potential function ψ(x,y). Therefore, a problem in an anti-plane state of deformation may be considered as a boundary value problem governed by Laplace’s equation, represented in Cartesian coordinates by Equation (8).
2 Ψ x 2 + 2 Ψ y 2 = 0 ,
which is valid in the whole body and satisfies boundary conditions along the surrounding contour. The fact that the same relationship also holds for the temperature field Ttemp(x,y) in plane steady-state heat conduction problems leads to the conclusion that thermal analogy may be used to obtain solutions of stress concentration factors for an anti–plane state of deformation. A more detailed description of this approach can be found in Reference [32] where SCFs for butt-welded plates subjected to shearing load were determined using FEM.
In the present analysis, the ANSYS 19 Multiphysics program with a Thermal module and PLANE55 finite element was used. A PLANE55 finite element is defined by four nodes with a single degree of freedom corresponding to the temperature at each node. Meshing of the modelled area was the same as in the previous cases for tensile and bending loads for identical joint geometry. The shape of the modelled body as well as mixed boundary conditions are shown in Figure 5.
Nominal uniform heat flux qnom = 1 W/m2 was applied over the right-hand edge of the body, while zero temperature was applied to the left end along the plane of symmetry of the joint. Since the upper face and the lower face of the joint are free from external shearing loads, they have to be insulated in the thermal model. It is clear that the shear stress components, τxz and τyz, given by Equation (9):
τ x z = G W x ; τ y z = G W y ,
and related to the partial derivatives of the potential function W(x,y) in particular directions, are proportional to the corresponding heat flux components, qx and qy, represented by Equation (10).
q x = k T temp x ; q y = k T temp y .
Numerical SCF values were calculated as a ratio of the maximum magnitude of the temperature gradient |∇Ttemp|max at the weld toe zone, shown by a small circle in Figure 5, to the magnitude of the nominal temperature gradient |∇Ttemp|nom over the right side of the body. It is well known that SCF values in such cases do not depend on the conductivity k of the medium. Therefore, the same result is obtained by comparing corresponding heat flux quantities qmax/qnom. In order to simplify relations between the heat flux and the magnitude of the temperature gradient, it was convenient to take k = 1 W/(Km).
One example of a steady-state heat conduction solution is shown in Figure 6 for the nominal heat flux qnom = 1 W/m2 applied over the right edge of the main plate. Magnitudes of the flux q, obtained at each point of the body, make it possible to find directly Kts value equal 1.54769, which is shown on the scale below the picture.
The maximum temperature corresponding to the unit nominal heat flux applied along the right end of the body was equal to 7.089 K. This temperature value resulted from the numerical solution of the boundary-value problem with mixed boundary conditions and was not directly applied over the right edge as a boundary condition.
In subsequent models of T-joints, particular values of X, Y, and θ were systematically changed, while the relative stiffener thickness T/a was constant and equal to one (Z0 = 1). One example of such a set of SCF results for a T-joint subjected to shearing load for θ = 45° is presented in Table 1.
Similar 12 sets of SCFs’ data were obtained for other weld angles θ in the range of 30–60° with a step of 2.5°. In this way, 3328 SCF numerical solutions were obtained for each loading mode, which means that 9984 cases of different T-joint shapes had to be modelled and solved for all loading modes.

4. SCF Approximation Formulas

4.1. Numerical and Graphical Representation of P Functions

After normalizing each SCF value with respect to a proper singular term X n , three approximating functions Pt, Pb, and Ps, corresponding to each loading mode, were determined in the form of polynomials. Their mathematical representation is shown in brackets of Equations (A1)–(A3) given in Appendix A. Graphical representations of these functions, corresponding to the weld angle θ = 45° for three loading modes, are shown in Figure 7, Figure 8 and Figure 9. Numerical values of the exponents n = −0.3264 and ns = −0.2 were obtained from Equations (6) and (7), respectively, for θ = π/4.

4.2. Validation of the P Functions

Accuracy of the approximating formulas of P functions was verified by comparing calculated SCFs values to their numerical equivalencies obtained using the finite element method. General formula, expressed by Equation (4), is now reduced to Equation (11).
K t = X n P ( X , Y , θ , Z 0 ) ,
where Z = Z0 = 1 and κ = 1 in the present case. Three examples of such comparisons are presented in Table 2, Table 3 and Table 4 for T-joints subjected to tensile, bending, and shearing loads, respectively, for θ = 45°.
Similar comparisons were made for various weld angles θ in the range of 30–60°. The maximum percentage error of the approximation was lower than 2%.

4.3. Determination of the Correction Functions κ

Additional FEM modelling was carried out in order to determine the influence of the relative attachment plate thickness Z = T/a on SCF, represented mathematically by a correction function κ. Several thousand new cases of a T-joint were modelled and solved for each loading mode and various proportions of T/a in the range of 1 ≤ T/a ≤ 4.
The general mathematical representation of the correction functions κ in the form of Equation (12).
κ ( X , Y , Z , θ , Z 0 ) = 1 + ( Z Z 0 ) [ 1 ( B 1 + B 2 Y 2 ) X m ] E x p [ ( B 3 Y ) p B 4 ]
was derived, which is similar to the one performed by the authors in Reference [30] for cruciform welded joints. Particular values of the exponents, m and p, and coefficients B1-B4, were determined for each loading mode using the least squares method. The best fitting of the functions κt, κb and κs, to the numerical FEM SCF results were obtained for m = 1, for tensile and bending load, and for m = 2 for shearing load, and for p equal to 2.4, 2.6, and 2.0 for tensile, bending, and shearing loads, respectively. Unfortunately, the coefficients Bi depended also on the weld angle θ. Therefore, additional approximations for Bi = fi(θ) were necessary. Mathematical representations of the correction functions κt, κb, and κs, suitable for each loading mode, are given in Appendix A. Some examples of the correction functions κt, κb, and κs for tensile, bending, and shearing loads, respectively, are shown in Figure 10, Figure 11 and Figure 12, for arbitrarily chosen parameters: θ = 45° and T/a equal to 2 and 4.

4.4. Validation of SCF Approximation Functions

Accuracy of the SCF approximation functions, given by the Formulas (A1)–(A3), containing correction functions κt, κb, and κs, > 1, was verified for all loading modes and for the weld angle θ varying in the range of 30–60°. Some examples of the validation are presented in Table 5, Table 6 and Table 7 for arbitrarily chosen angles θ equal to 30°, 45°, and 55°.
Similar validation was made for other weld angles θ with a step of 2.5°. The maximum percentage error found for all analysed cases was lower than 2% with respect to the SCFs’ results obtained numerically using the FEM. Due to the fact that the weld angle θ = 45° is commonly used in engineering applications of such T-joints, extended parametric Equations (A1)–(A3) given in Appendix A, were reduced to Equations (A4)–(A6) presented in Appendix B.

5. Discussion

Several thousand SCF solutions obtained using the FEM as well as corresponding parametric Equations (A1)–(A3) of high accuracy and a wide range of validity, made it possible to draw some more general conclusions regarding the influence of particular geometrical parameters on SCF. The general rule is that particular characteristic geometrical details of the joint located closer to the reference point of the maximum stress σ1max affect this value more significantly than others. In order to explain this phenomenon more precisely, three qualitatively different geometrical cases of a T-joint are shown in Figure 13. Additionally, some alternative, parametric SCF formulas available in the literature [13,15,16,17,19] and related to T-joints are also given in Appendix C.
In case A, shown in Figure 13, the main plate is relatively thin compared to other dimensions of the joint, tA << L, while the magnitude of the weld toe radius ρ is of the same order as the plate thickness tA. In such a case, SCF values depend mainly on the tA/ρ ratio. For example, when tA/ρ changes in the range of 0.8–5.0, Ktt varies from 1.16 to 1.63 and does not depend on the weld angle θ in the range of 30–60°. Similar behaviour is observed for tA/ρ in the range of 0.8–4.0, where Ktb varies from 1.16 to 1.66. In case B, when the weld toe radius ρ is much smaller than in the previous case and L >> tB, SCF values should depend on the quantities tB/ρ and θ.
Equations (A9) and (A12) presented in Appendix C and numerical SCF data shown in Table 8 confirm such a conclusion.
However, accuracy of the formulas given in Appendix C depends on the mathematical representation chosen by the authors, assumed range of validity, and the accuracy of data used in approximations. Some comparisons of SCFs obtained from various equations to the FEM results carried out in the present work are shown in Table 9 and Table 10 for θ = 45°.
In case C, shown in Figure 13, mutual relations between L and tC become important, especially when the L is much shorter than tC. The presence of such a quantity, L/tC, can be observed directly in Equations (A10) and (A14). In a general case, such an influencing parameter depends on the thickness of the stiffener, T, leg lengths, weld angle θ, and the weld toe radius ρ. Therefore, it seems reasonable to use all these quantities in order to derive “effective length” of the joint and, as a consequence, reach better SCF approximations than by a simple use of L. It is also worth noting that SCFs used in Equations (A10) and (A14) are defined in a different manner than in other formulas presented in this work. They are related to the maximum longitudinal stress component at the weld toe and not to the maximum principal stress σ1max. Therefore, SCF values obtained from Equations (A10) and (A14) and shown in Table 9 and Table 10, are a few percent underestimated compared with appropriate FEM results.
Particular data given in Table 9 and Table 10 represent only some comparative examples of SCF values with respect to the FEM results. Therefore, these results cannot be generalised to assess accuracy of each formula. Some of these results denoted by (#) mean that they formally lie out of range of validity assumed by the authors. However, in many cases, the accuracy remains satisfactory.
It may be assumed that the range of validity of the formulas (A8) and (A13), proposed by Tsuji [17], is similar to that for Equations (A7) and (A11) reported by Ushirokawa and Nakayama [16].

6. Conclusions

Systematic numerical FEM modelling including more than 22,600 cases of welded plate T-joints made it possible to derive three approximating SCF’s formulas for axial, bending, and shearing loads. Five geometrical parameters: ρ, a, θ, t, and T were considered as independent variables varying in relatively wide ranges 0 < ρ/a ≤ 1.3, 0 < a/t ≤ 1.3, 1 ≤ T/a ≤ 4, and 30° ≤ θ ≤ 60° covering dimensional proportions of welded plate T-joints used in engineering applications. Two limiting cases are also included in the solutions—when the weld toe radius tends to zero and when the main plate thickness becomes infinite. The accuracy of the formulas is better than 98% compared with numerical FEM results. The use of the proper exponents, n and ns, corresponding to the stress field singularities of a sharp corner dependent on the angle θ, was the key condition in obtaining such an accuracy and range of validity of proposed parametric equations. In the cases of shearing loads, a plane FEM model based on thermal analogy was successfully used.
In spite of the fact that SCF alone is not sufficient for calculating fatigue strength of a structural element, all the parametric formulas given in a closed form may be used as a part of computer-aided procedures for the fatigue strength assessment of welded plate T-joints.
Presented solutions may be applied as a tool for a computer aided assessment in fatigue design of welded plate T-joints. Particularly in:
(a)
comparative studies of stress concentration in various geometrical forms of welded T-joints, including the use of fictitious values of the weld toe radii,
(b)
hot spot method applied in fatigue design of such joints,
(c)
weight function method used in fracture mechanics models for crack initiated at the weld toe,
(d)
dealing with the necessity of additional mechanical improvements in the weld toe region.

Author Contributions

Conceptualization, K.L.M. Methodology, K.L.M. Software, P.T. Validation, K.L.M. and P.T. Formal analysis, K.L.M. and P.T. Investigation, K.L.M. and P.T. Resources, K.L.M. and P.T. Data curation, P.T. Writing—original draft preparation, K.L.M. Writing—review and editing, K.L.M. and P.T. Visualization, K.L.M. and P.T. Supervision, K.L.M. Project administration, P.T. Funding acquisition, K.L.M. and P.T. 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

Data is contained within the article.

Acknowledgments

This research was co-founded through a subsidy of the Ministry of Science and Higher Education for the discipline of mechanical engineering at the Faculty of Mechanical Engineering, Bialystok University of Technology, project no. WZ/WM-IIM/4/2020.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

atheoretical weld throat thickness
Gshear modulus
h, hpleg lengths defining the weld size
kthermal conductivity
Ktstress concentration factor (SCF)
Kttstress concentration factor for tensile (axial) load
Ktbstress concentration factor for bending load
Ktsstress concentration factor for shearing load
m, pexponents dependent on the loading mode
nstress field exponent for a sharp corner for axial and bending load
nsstress field exponent for a sharp corner subjected to shearing load
Pregular function represented by polynomials
Pt, Pb, Psfunctions P corresponding to tensile, bending, and shearing load, respectively
qmagnitude of the heat flux
qmaxmagnitude of the maximum heat flux
qnommagnitude of the nominal heat flux at the end of the body
tthickness of the main plate
Tthickness of the attachment plate
Ttemptemperature
Wdisplacement component corresponding to the anti–plane shear
X = ρ/(ρ + a)normalized toe radius parameter
Y = a/(a + t)normalized weld thickness parameter
Z = T/anormalized attachment plate thickness parameter
Z0arbitrarily chosen reference T/a value
|∇Ψ|magnitude of the displacement (temperature) gradient
|∇Ψ|maxmagnitude of the maximum displacement (temperature) gradient
|∇Ψ|nommagnitude of the nominal displacement (temperature) gradient
θweld angle
κcorrection function for the relative attachment plate thickness T/a
κtcorrection function κ for tensile (axial) load
κbcorrection function κ for bending load
κscorrection function κ for shearing load
Ψ(x,y)potential function representing temperature or anti–plane displacement
∂Ψ/nnormal derivative at the bounding contour
ρweld toe radius
σ1first principal stress produced by normal or bending load
σ1maxmaximum value of the first principal stress at the weld toe zone
σtnominal tensile (axial) stress
σbnominal bending stress
τmaxmaximum shear stress at the weld toe due to shearing load
τsnominal shear stress

Appendix A

Approximating SCFs’ Formulas for a Welded Plate T-Joint Subjected to Axial, Bending, and Shearing Loads

Independent variables: ρ, a, θ, t, T;
Normalized quantities: X = ρ/(ρ + a), Y = a/(a + t); Z = T/a;
Range of validity: 0 < ρ/a ≤ 1.3; 0 < a/t ≤ 1.3; 1 ≤ T/a ≤ 4; 30° ≤ θ ≤ 60°.
In all formulas, θ is expressed in radians.
Exponent of the singular term for:
Axial and bending load:
n = 0.63662 θ 0.09330 θ 2 1 + 0.77635 θ + 0.04075 θ 1.5 0.00499 θ 2 + 0.13365 θ 2.5
Shearing load:
n s = θ θ + π
Accuracy: maximum percentage error lower than 2% compared to the FEM results
Tensile (axial) load:
K t t = X n ( A 0 t + A 1 t X + A 2 t X 2 + A 3 t X 3 + A 4 t X 4 ) κ t
where:
A 0 t = A 00 t + A 01 t Y + A 02 t Y 2 + A 03 t Y 3 + A 04 t Y 4
A 00 t = 2.078 0.712 θ 0.076 θ 4
A 01 t = 0.132 + 0.718 θ 0.455 θ 4
A 02 t = 18.982 + 12.585 θ + 0.398 θ 4
A 03 t = 55.711 54.642 θ + 5.304 θ 4
A 04 t = 47.047 + 53.604 θ 7.139 θ 4
A 1 t = A 10 t + A 11 t Y + A 12 t Y 2 + A 13 t Y 3 + A 14 t Y 4
A 10 t = 0.066 0.789 θ + 0.878 θ 4
A 11 t = 0.413 + 0.119 θ 2 + 0.428 θ 4
A 12 t = 6.193 5.495 θ 2 5.077 θ 4
A 13 t = 20.187 + 34.745 θ 2 + 11.092 θ 4
A 14 t = 16.393 27.986 θ 2 13.135 θ 4
A 2 t = A 20 t + A 21 t Y + A 22 t Y 2 + A 23 t Y 3 + A 24 t Y 4
A 20 t = 5.133 21.927 θ + 24.944 θ 2 8.229 θ 4
A 21 t = 2.250 2.429 θ 2 + 0.805 θ 4
A 22 t = 5.156 6.961 θ 2 + 14.020 θ 4
A 23 t = 0.909 + 92.878 θ 2 118.392 θ 4
A 24 t = 16.571 147.711 θ 2 + 151.148 θ 4
A 3 t = A 30 t + A 31 t Y + A 32 t Y 2 + A 33 t Y 3 + A 34 t Y 4
A 30 t = 15.018 + 58.059 θ 60.616 θ 2 + 17.595 θ 4
A 31 t = 7.053 + 5.113 θ 0.340 θ 4
A 32 t = 14.167 + 8.281 θ 2 22.438 θ 4
A 33 t = 19.091 213.131 θ 2 + 226.174 θ 4
A 34 t = 146.976 + 316.815 θ 195.919 θ 4
A 4 t = A 40 t + A 41 t Y + A 42 t Y 2 + A 43 t Y 3 + A 44 t Y 4
A 40 t = 10.494 40.594 θ + 41.995 θ 2 11.917 θ 4
A 41 t = 24.260 73.105 θ + 67.325 θ 2 17.427 θ 4
A 42 t = 1.928 16.706 θ 2 + 18.955 θ 4
A 43 t = 86.411 + 181.383 θ 108.284 θ 4
A 44 t = 117.729 227.646 θ + 117.488 θ 4
κ t = 1 + ( Z 1 ) { 1 [ ( 0.889 + 2.279 θ 0.539 θ 2 ) + ( 12.70 + 10.21 θ 7.17 θ 2 ) Y 2 ] X } E x p { [ ( 12.94 13.94 θ + 6.57 θ 2 ) Y ] 2.4 ( 3.72 4.03 θ + 1.62 θ 2 ) }
Bending load:
K t b = X n ( A 0 b + A 1 b X + A 2 b X 2 + A 3 b X 3 + A 4 b X 4 ) κ b
where:
A 0 b = A 00 b + A 01 b Y + A 02 b Y 2 + A 03 b Y 3 + A 04 b Y 4
A 00 b = 1.833 0.316 θ 2 0.621 θ 3 + 0.394 θ 4
A 01 b = 1.282 + 6.636 θ 10.422 θ 3 + 5.974 θ 4
A 02 b = 16.721 7.442 θ 2 + 54.668 θ 3 33.383 θ 4
A 03 b = 50.505 118.407 θ 2 + 50.936 θ 3 + 12.039 θ 4
A 04 b = 43.771 + 162.845 θ 2 140.901 θ 3 + 30.243 θ 4
A 1 b = A 10 b + A 11 b Y + A 12 b Y 2 + A 13 b Y 3 + A 14 b Y 4
A 10 b = 0.015 0.811 θ 0.974 θ 2 + 1.765 θ 3
A 11 b = 0.585 + 0.319 θ 0.084 θ 4
A 12 b = 7.287 + 53.653 θ 55.081 θ 2 + 0.947 θ 4
A 13 b = 5.158 77.965 θ + 105.085 θ 2
A 14 b = 28.354 41.874 θ 2
A 2 b = A 20 b + A 21 b Y + A 22 b Y 2 + A 23 b Y 3 + A 24 b Y 4
A 20 b = 2.501 11.722 θ + 14.711 θ 2 5.338 θ 4
A 21 b = 20.181 60.484 θ + 51.074 θ 2 14.228 θ 3
A 22 b = 15.157 0.689 θ 2 + 35.741 θ 4
A 23 b = 74.171 + 0.421 θ 2 89.665 θ 4
A 24 b = 108.419 + 93.296 θ 2 + 1.340 θ 4
A 3 b = A 30 b + A 31 b Y + A 32 b Y 2 + A 33 b Y 3 + A 34 b Y 4
A 30 b = 21.534 + 82.796 θ 94.723 θ 2 + 18.151 θ 3 + 14.663 θ 4
A 31 b = 12.022 + 42.247 θ 2 16.989 θ 4
A 32 b = 68.318 111.122 θ 2 28.428 θ 4
A 33 b = 268.940 + 340.766 θ 2 + 18.190 θ 4
A 34 b = 342.766 505.198 θ 2 + 160.946 θ 4
A 4 b = A 40 b + A 41 b Y + A 42 b Y 2 + A 43 b Y 3 + A 44 b Y 4
A 40 b = 30.817 118.209 θ + 137.515 θ 2 34.910 θ 3 14.672 θ 4
A 41 b = 6.060 51.272 θ 3 + 33.481 θ 4
A 42 b = 188.380 + 368.847 θ 453.325 θ 3 + 326.318 θ 4
A 43 b = 534.753 856.175 θ + 926.225 θ 3 645.821 θ 4
A 44 b = 690.666 + 1465.07 θ 1261.73 θ 2 + 396.370 θ 3 + 50.486 θ 4
κ b = 1 + ( Z 1 ) { 1 [ ( 1.00 + 2.23 θ 0.41 θ 2 ) + ( 2.81 + 37.10 θ 21.04 θ 2 ) Y 2 ] X } E x p { [ ( 11.77 13.20 θ + 5.77 θ 2 ) Y ] 2.6 ( 3.84 4.33 θ + 1.68 θ 2 ) }
Shearing load:
K t s = X n s ( A 0 s + A 1 s X + A 2 s X 2 + A 3 s X 3 + A 4 s X 4 ) κ s
where:
A 0 s = 1.4361 0.0912 θ 2 + ( 0.8777 0.0080 θ 2 ) Y 2
A 1 s = 0.1147 0.6461 θ + 0.2553 θ 2 + ( 0.0581 + 0.1094 θ 2 ) Y 2
A 2 s = 0.5070 + 0.4287 θ 2 + ( 0.4582 + 0.2199 θ 2 ) Y 2
A 3 s = 0.7581 0.4544 θ 2 + ( 0.7112 0.1743 θ 2 ) Y 2
A 4 s = 0.6625 + 0.4349 θ 2 + ( 1.1281 0.5013 θ 2 ) Y 2
κ s = 1 + ( Z 1 ) [ 1 ( ( 0.40 + 0.67 θ + 0.70 θ 2 ) + ( 4.17 + 18.54 θ 6.94 θ 2 ) Y 2 ) X 2 ] E x p [ ( ( 6.26 5.74 θ + 2.52 θ 2 ) Y ) 2.0 ( 3.84 3.31 θ + 1.23 θ 2 ) ]

Appendix B

Parametric Formulas for Calculating SCFs, while θ = 45°

Tensile (axial) load:
K t t = X 0.3264 ( A 0 t + A 1 t X + A 2 t X 2 + A 3 t X 3 + A 4 t X 4 ) κ t
where:
A 0 t = 1.490 + 0.523 Y 8.946 Y 2 + 14.813 Y 3 7.663 Y 4
A 1 t = 0.352 0.177 Y + 0.872 Y 2 + 5.466 Y 3 5.868 Y 4
A 2 t = 0.167 + 1.0580 Y 4.115 Y 2 + 13.152 Y 3 17.032 Y 4
A 3 t = 0.115 3.167 Y + 10.737 Y 2 26.319 Y 3 + 27.302 Y 4
A 4 t = 0.018 + 1.742 Y 5.021 Y 2 + 14.844 Y 3 16.359 Y 4
κ t = 1 + ( Z 1 ) [ 1 ( 0.568 + 16.296 Y 2 ) X ] · E x p [ ( 6.044 Y ) 2.4 1.554 ]
Bending load:
K t b = X 0.3264 ( A 0 b + A 1 b X + A 2 b X 2 + A 3 b X 3 + A 4 b X 4 ) κ b
where:
A 0 b = 1.487 + 1.154 Y 7.529 Y 2 + 6.724 Y 3 0.075 Y 4
A 1 b = 0.368 0.366 Y + 1.236 Y 2 1.570 Y 3 + 2.524 Y 4
A 2 b = 0.338 2.711 Y 1.982 Y 2 + 40.313 Y 3 50.359 Y 4
A 3 b = 0.563 + 7.574 Y 11.045 Y 2 51.817 Y 3 + 92.375 Y 4
A 4 b = 0.306 6.040 Y + 5.853 Y 2 + 65.308 Y 3 107.060 Y 4
κ b = 1 + ( Z 1 ) [ 1 ( 0.499 + 13.350 Y 2 ) X ] · E x p [ ( 4.962 Y ) 2.6 1.476 ]
Shearing load:
K t s = X 0.2 ( A 0 s + A 1 s X + A 2 s X 2 + A 3 s X 3 + A 4 s X 4 ) κ s
where:
A 0 s = 1.380 0.883 Y 2
A 1 s = 0.235 + 0.126 Y 2
A 2 s = 0.243 + 0.594 Y 2
A 3 s = 0.478 0.819 Y 2
A 4 s = 0.394 + 0.819 Y 2
κ s = 1 + ( Z 1 ) [ 1 ( 0.558 + 6.110 Y 2 ) X 2 ] · E x p [ ( 3.306 Y ) 2 1.999 ]

Appendix C

Alternative Parametric Formulas for Calculating SCF for Plate T-Joints Subjected to Tensile and Bending Loads

Plate T-joint—axial load
Ushirokawa and Nakayama (1983). Range of validity: 0.025 ≤ ρ/t ≤ 0.35, 20° ≤ θ ≤ 50°, T/t = 1, hp/t = 0.75.
K t t = 1 + [ h / ρ 2.8 ( w / t ) 2 ] 0.65 { 1 E x p [ 0.9 θ w / 2 h ] 1 E x p [ 0.9 ( π / 2 ) w / 2 h ] }
w = ( t + 2 h ) + 0.3 ( T + 2 h p ) .
Tsuji (1990): range of validity not precisely defined.
K t t = 1 + 1.015 [ h / ρ 2.8 ( w / t ) 2 ] 0.446 { 1 E x p [ 0.9 θ w / 2 h ] 1 E x p [ 0.9 ( π / 2 ) w / 2 h ] }
w = ( t + 2 h ) + 0.3 ( T + 2 h p ) .
Monahan (1995), range of validity: 0.02 ≤ ρ/t ≤ 0.066, 30° ≤ θ≤60°, L/t = 2.8.
K t t = 1 + 0.388 ( θ ) 0.37 ( ρ / t ) 0.454
Brennan et al. (2000), range of validity: 0.01 ≤ ρ/t ≤ 0.066, 30° ≤ θ≤60°, 0.3 ≤ L/t ≤ 4.0
K t t = 1.027 + 0.271 ( θ ) 0.216 ( ρ / t ) 0.47 ( L / t ) 0.183
Plate T-joint—bending load
Ushirokawa and Nakayama (1983), range of validity: 0.025 ≤ ρ/t ≤ 0.35, 20° ≤ θ≤50°, T/t = 1, hp/t = 0.75.
K t b = 1 + 1.9 T a n h ( 2 T t + 2 h + 2 ρ t ) { T a n h [ ( 2 h / t ) 1 / 4 1 ρ / t ] } [ 0.13 + 0.65 ( 1 ρ / t ) 4 ( ρ / t ) 1 / 3 ] { 1 E x p [ 0.9 θ w / 2 h ] 1 E x p [ 0.9 ( π / 2 ) w / 2 h ] }
w = ( t + 2 h ) + 0.3 ( T + 2 h p )
Niu and Glinka (1987), range of validity: 0.02 ≤ ρ/t ≤ 0.066, 30° ≤ θ≤60°, T/t = 1, h/t = 1.
K t b = 1 + 0.5121 ( θ ) 0.572 ( ρ / t ) 0.469
Tsuji (1990), range of validity not precisely defined.
K t b = 1 + [ 0.629 + 0.058 L n ( T + 2 h p t ) ] [ ρ t ] 0.431 T a n h ( 6 h t ) { 1 E x p [ 0.9 θ w / 2 h ] 1 E x p [ 0.9 ( π / 2 ) w / 2 h ] }
w = ( t + 2 h ) + 0.3 ( T + 2 h p )
Brennan et al. (2000), range of validity: 0.01 ≤ ρ/t ≤ 0.066, 30° ≤ θ≤60°, 0.3 ≤ L/t ≤ 4.0.
K t t = 1.01 + 0.344 ( θ ) 0.336 ( ρ / t ) 0.468 ( L / t ) 0.233

References

  1. Gurney, T.R. Fatigue of Welded Structures, 2nd ed.; Cambridge University Press: Cambridge, UK, 1979. [Google Scholar]
  2. European Committee for Standardization (CES). EN 1993-1-9:2005 Eurocode 3: Design of Steel Structures—Part. 1–9: Fatigue; CES: Brussels, Belgium, 2005. [Google Scholar]
  3. Hobbacher, A.F. The new IIW recommendations for fatigue assessment of welded joints and components—A comprehensive code recently updated. Int. J. Fatigue 2009, 31, 50–58. [Google Scholar] [CrossRef]
  4. Fricke, W. IIW Recommendations for the Fatigue Assessment of Welded Structures by Notch Stress Analysis; Woodhead Publishing: Cambridge, UK, 2012; pp. 2–41. [Google Scholar]
  5. Kranz, B.; Sonsino, C.M. Verification of FAT Values for the Application of the Notch Stress Concept with the Reference Radii Rref = 1.00 and 0.05 mm. Weld. World 2010, 54, R218–R224. [Google Scholar] [CrossRef]
  6. Fricke, W. IIW guideline for the assessment of weld root fatigue. Weld. World 2013, 57, 753–791. [Google Scholar] [CrossRef]
  7. International Standard Organization (ISO). EN ISO 9692-1:2013 Welding and Allied Processes—Types of Joint Preparation—Part. 1: Manual Metal. Arc Welding, Gas.-Shielded Metal. Arc Welding, Gas. Welding, TIG Welding and Beam Welding of Steels; International Organization for Standardization: Geneva, Switzerland, 2013. [Google Scholar]
  8. Zhao, X.L.; Packer, J.A. Fatigue Design Procedure for Welded Hollow Section Joints, 1st ed.; Abington Publishing: Nashville, TN, USA, 2000. [Google Scholar]
  9. Mashiri, R.; Zhao, X.L.; Grundy, P. Stress concentration factors and fatigue failure of welded T-connections in circular hollow sections under in-plane bending. Int. J. Struct. Stab. Dyn. 2004, 4, 403–422. [Google Scholar] [CrossRef]
  10. Kršćanski, S.; Turkalj, G. FEM stress concentration factors for fillet welded CHS-plate T-joint. Eng. Rev. 2012, 32, 147–155. [Google Scholar]
  11. Saini, D.S.; Karmakar, D.; Ray-Chaudhuri, S. A review of stress concentration factors in tubular and non-tubular joints for design of offshore installations. JOES 2016, 1, 186–202. [Google Scholar] [CrossRef] [Green Version]
  12. Carpinteri, A.; Boaretto, J.; Fortese, G.; Giordani, F.; Iturrioz, I.; Ronchei, C.; Scorza, D.; Vantadori, S. Fatigue life estimation of fillet-welded tubular T-joints subjected to multiaxial loading. Int. J. Fatigue 2017, 101, 263–270. [Google Scholar] [CrossRef]
  13. Niu, X.; Glinka, G. The weld profile effect on stress intensity factors in weldments. Int J. Fract. 1987, 35, 3–20. [Google Scholar]
  14. Bhuyan, G.; Vosikovsky, O. Prediction of fatigue crack initiation lives for welded plate T-joints based on the local stress-strain approach. Int. J. Fatigue 1989, 11, 153–159. [Google Scholar] [CrossRef]
  15. Monahan, C.C. Early Fatigue Crack Growth at Welds; WIT PRESS/Computational Mechanics Publications: Southampton, UK, 1995. [Google Scholar]
  16. Ushirokawa, O.; Nakayama, E. Stress concentration factor at welded joints. Ishikawajima Harima Eng. Rev. 1983, 23, 351–355. [Google Scholar]
  17. Tsuji, I. Estimation of stress concentration factor at weld toe of non-load carrying fillet welded joints. Trans. West. Jpn. Soc. Nav. Archit. 1990, 80, 241–251. [Google Scholar]
  18. Iida, K.; Uemura, T. Stress concentration factor formulas widely used in Japan. Fatigue Fract. Eng. Mater. Struct. 1996, 19, 779–786. [Google Scholar] [CrossRef]
  19. Brennan, F.P.; Peletiesa, P.; Hellierb, A.K. Predicting weld toe stress concentration factors for T and skewed T-joint plate connections. Int. J. Fatigue 2000, 22, 573–584. [Google Scholar] [CrossRef]
  20. Dabiri, M.; Ghafouri, M.; Raftar, H.R.R.; Björk, T. Neural network-based assessment of the stress concentration factor in a T-welded joint. J. Constr. Steel Res. 2017, 128, 567–578. [Google Scholar] [CrossRef]
  21. Bell, R.; Vosikovsky, O.; Bain, S.A. The significance of weld toe undercuts in the fatigue of steel plate T-joints. Int. J. Fatigue 1989, 11, 3–11. [Google Scholar] [CrossRef]
  22. Hou, C.Y. Fatigue analysis of welded joints with the aid of real three-dimensional weld toe geometry. Int. J. Fatigue 2007, 29, 772–785. [Google Scholar] [CrossRef]
  23. Tarasiuk, P.; Molski, K.L.; Szymaniuk, A. Fatigue designing of welded agricultural wheels. Eksploat. Niezawodn. 2013, 15, 123–128. [Google Scholar]
  24. Zhao, W.; Hsu, W.-T. Re-Evaluation of Fatigue Thickness Effect Based on Fatigue Test Database. J. Mar. Sci. Eng. 2020, 8, 895. [Google Scholar] [CrossRef]
  25. Wang, Y.; Luo, Y.; Tsutsumi, S. Parametric Formula for Stress Concentration Factor of Fillet Weld Joints with Spline Bead Profile. Materials 2020, 13, 4639. [Google Scholar] [CrossRef]
  26. Terán, G.; Albiter, A.; Cuamatzi-Meléndez, R. Parametric evaluation of the stress concentration factors in T-butt welded connections. Eng. Struct. 2013, 56, 1484–1495. [Google Scholar] [CrossRef]
  27. Costa, J.D.M.; Jesus, J.S.; Loureiro, A.; Ferreira, J.A.M.; Borrego, L.P. Fatigue life improvement of mig welded aluminium T-joints by friction stir processing. Int. J. Fatigue 2014, 61, 244–254. [Google Scholar] [CrossRef]
  28. Luo, Y.; Ma, R.; Tsutsumi, S. Parametric formulae for elastic stress concentration factor at the weld toe of distorted butt-welded joints. Materials 2020, 13, 169. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Stenberg, T.; Barsoum, Z.; Hedlund, J.; Josefsson, J. Development of a computational fatigue model for evaluation of weld quality. Weld. World 2019, 63, 1771–1785. [Google Scholar] [CrossRef] [Green Version]
  30. Haagensen, P.J.; Maddox, S.J. Recommendations on Methods for Improving the Fatigue Strength of Welded Joints; Woodhead Publishing: Cambridge, UK, 2013; ISBN 9781782420644. [Google Scholar]
  31. Molski, K.L.; Tarasiuk, P.; Glinka, G. Stress concentration at cruciform welded joints under axial and bending loading modes. Weld. World 2020, 64, 1867–1876. [Google Scholar] [CrossRef]
  32. Molski, K.L.; Tarasiuk, P. Stress Concentration Factors for Butt-Welded Plates Subjected to Tensile, Bending and Shearing Loads. Materials 2020, 13, 1798. [Google Scholar] [CrossRef] [Green Version]
  33. Seweryn, A.; Molski, K.L. Elastic stress singularities and corresponding generalized stress intensity factors for angular corners under various boundary conditions. Eng. Fract. Mech. 1996, 55, 529–556. [Google Scholar] [CrossRef]
Figure 1. Geometrical parameters and loading modes of the plate T-joint.
Figure 1. Geometrical parameters and loading modes of the plate T-joint.
Materials 14 00546 g001
Figure 2. Basic geometrical parameters and boundary conditions of the numerical FEM model for tensile and bending loads.
Figure 2. Basic geometrical parameters and boundary conditions of the numerical FEM model for tensile and bending loads.
Materials 14 00546 g002
Figure 3. Example of a finite element mesh for θ = 45°, while X = 0.25, Y = 0.415, and T/t = 1.
Figure 3. Example of a finite element mesh for θ = 45°, while X = 0.25, Y = 0.415, and T/t = 1.
Materials 14 00546 g003
Figure 4. Distribution of the first principal stress σ1 for σt = 1 MPa, while θ = 45°, X = 0.25, Y = 0.415, and T/t = 1.
Figure 4. Distribution of the first principal stress σ1 for σt = 1 MPa, while θ = 45°, X = 0.25, Y = 0.415, and T/t = 1.
Materials 14 00546 g004
Figure 5. Shape and boundary conditions of the numerical FEM model for shearing load.
Figure 5. Shape and boundary conditions of the numerical FEM model for shearing load.
Materials 14 00546 g005
Figure 6. Heat flux q distribution, interpreted as shearing stress τ for τs = 1 MPa, while θ = 45°, X = 0.25, Y = 0.415 and T/t = 1.
Figure 6. Heat flux q distribution, interpreted as shearing stress τ for τs = 1 MPa, while θ = 45°, X = 0.25, Y = 0.415 and T/t = 1.
Materials 14 00546 g006
Figure 7. Graphical representation of the function P t = K t t / X n for θ = 45°, n = −0.3264 and Z0 = 1.
Figure 7. Graphical representation of the function P t = K t t / X n for θ = 45°, n = −0.3264 and Z0 = 1.
Materials 14 00546 g007
Figure 8. Graphical representation of the function P b = K t b / X n for θ = 45°, n = −0.3264 and Z0 = 1.
Figure 8. Graphical representation of the function P b = K t b / X n for θ = 45°, n = −0.3264 and Z0 = 1.
Materials 14 00546 g008
Figure 9. Graphical representation of the function P s = K t s / X n s for θ = 45°, ns = −0.2 and Z0 = 1.
Figure 9. Graphical representation of the function P s = K t s / X n s for θ = 45°, ns = −0.2 and Z0 = 1.
Materials 14 00546 g009
Figure 10. Correction function κt for tensile load, while θ = 45°.
Figure 10. Correction function κt for tensile load, while θ = 45°.
Materials 14 00546 g010
Figure 11. Correction function κb for the bending load while θ = 45°.
Figure 11. Correction function κb for the bending load while θ = 45°.
Materials 14 00546 g011
Figure 12. Correction function κs for the shearing load while θ = 45°.
Figure 12. Correction function κs for the shearing load while θ = 45°.
Materials 14 00546 g012
Figure 13. Characteristic dimensional proportions representing the influence of characteristic geometrical dimensions on SCF.
Figure 13. Characteristic dimensional proportions representing the influence of characteristic geometrical dimensions on SCF.
Materials 14 00546 g013
Table 1. Representation of SCF values for a plate T-joint subjected to shearing load, while θ = 45° and T/a = 1.
Table 1. Representation of SCF values for a plate T-joint subjected to shearing load, while θ = 45° and T/a = 1.
θ = 45°X
Y0.0100.0180.0250.0320.0500.0790.1000.1500.2000.2500.3000.3500.4000.4500.5000.562
0.0103.4743.0932.8862.7502.5002.2662.1541.9661.8371.7381.6571.5881.5281.4751.4271.372
0.0183.4723.0922.8852.7492.5002.2652.1541.9661.8371.7381.6571.5881.5281.4751.4271.372
0.0253.4703.0912.8842.7492.4992.2652.1541.9661.8361.7371.6561.5881.5281.4751.4261.372
0.0323.4693.0902.8842.7482.4992.2652.1541.9661.8361.7371.6561.5881.5281.4741.4261.371
0.0503.4643.0882.8812.7462.4972.2632.1521.9641.8351.7361.6551.5861.5271.4731.4251.370
0.0793.4553.0792.8762.7402.4922.2582.1481.9601.8301.7321.6521.5831.5231.4701.4221.367
0.1003.4533.0762.8712.7352.4872.2542.1431.9561.8271.7291.6481.5801.5201.4671.4191.364
0.1503.4293.0532.8492.7152.4682.2372.1271.9411.8131.7151.6361.5681.5081.4561.4081.353
0.2003.3893.0182.8162.6832.4402.2112.1021.9191.7921.6951.6161.5491.4911.4391.3911.337
0.2503.3352.9702.7712.6412.4012.1762.0691.8881.7641.6691.5911.5251.4671.4161.3691.316
0.3003.2692.9112.7152.5882.3532.1322.0281.8501.7291.6351.5591.4951.4381.3881.3431.291
0.3503.1912.8412.6502.5262.2962.0811.9791.8061.6881.5971.5231.4601.4051.3571.3131.263
0.4003.1042.7632.5782.4572.2342.0241.9251.7571.6421.5541.4831.4221.3691.3231.2811.233
0.4503.0092.6792.5002.3832.1661.9631.8671.7051.5941.5091.4401.3821.3321.2881.2481.204
0.5002.9112.5922.4182.3042.0951.8991.8071.6501.5431.4621.3961.3411.2941.2521.2161.175
0.5622.7822.4772.3122.2032.0031.8161.7281.5791.4781.4021.3411.2901.2471.2101.1771.142
Table 2. Comparison of Ktt values calculated using Equation (A1) (*) to the corresponding FEM results for various X and Y, while θ = 45° and T/a = 1.
Table 2. Comparison of Ktt values calculated using Equation (A1) (*) to the corresponding FEM results for various X and Y, while θ = 45° and T/a = 1.
θ = 45°X
Y0.0500.1000.2000.3000.4000.5000.562
0.0503.916 3.930 *3.091 3.099 *2.412 2.419 *2.070 2.075 *1.845 1.850 *1.679 1.683 *1.593 1.596 *
0.1003.848 3.854 *3.034 3.039 *2.369 2.372 *2.033 2.035 *1.812 1.813 *1.648 1.647 *1.562 1.561 *
0.2003.540 3.530 *2.793 2.785 *2.183 2.177 *1.875 1.870 *1.672 1.668 *1.520 1.517 *1.440 1.439 *
0.3003.099 3.108 *2.451 2.459 *1.927 1.933 *1.666 1.671 *1.496 1.500 *1.371 1.374 *1.307 1.309 *
0.4002.696 2.702 *2.145 2.150 *1.708 1.711 *1.496 1.497 *1.361 1.361 *1.263 1.262 *1.214 1.212 *
0.5002.379 2.374 *1.910 1.906 *1.545 1.542 *1.373 1.371 *1.265 1.264 *1.189 1.190 *1.152 1.153 *
0.5622.209 2.216 *1.784 1.790 *1.460 1.464 *1.309 1.311 *1.216 1.218 *1.152 1.152 *1.121 1.121 *
Table 3. Comparison of Ktb values calculated using Equation (A2) (*) to the corresponding FEM results for various X and Y, while T/a = 1 and θ = 45°.
Table 3. Comparison of Ktb values calculated using Equation (A2) (*) to the corresponding FEM results for various X and Y, while T/a = 1 and θ = 45°.
θ = 45°X
Y0.0500.1000.2000.3000.4000.5000.562
0.0503.997 4.010 *3.149 3.160 *2.455 2.463 *2.103 2.110 *1.872 1.878 *1.700 1.706 *1.611 1.616 *
0.1004.020 4.027 *3.164 3.171 *2.462 2.467 *2.104 2.108 *1.868 1.872 *1.692 1.695 *1.601 1.602 *
0.2003.868 3.856 *3.041 3.033 *2.358 2.351 *2.007 2.001 *1.773 1.768 *1.596 1.593 *1.505 1.501 *
0.3003.494 3.502 *2.746 2.753 *2.127 2.131 *1.806 1.813 *1.595 1.602 *1.440 1.444 *1.357 1.363 *
0.4003.059 3.072 *2.411 2.416 *1.873 1.876 *1.599 1.604 *1.423 1.428 *1.298 1.302 *1.236 1.240 *
0.5002.685 2.671 *2.118 2.106 *1.659 1.649 *1.435 1.426 *1.297 1.290 *1.204 1.198 *1.161 1.157 *
0.5622.469 2.482 *1.952 1.963 *1.540 1.548 *1.346 1.352 *1.231 1.236 *1.157 1.161 *1.123 1.126 *
Table 4. Comparison of Kts values calculated using Equation (A3) (*) to the corresponding FEM results for various X and Y, while T/a = 1 and θ = 45°.
Table 4. Comparison of Kts values calculated using Equation (A3) (*) to the corresponding FEM results for various X and Y, while T/a = 1 and θ = 45°.
θ = 45°X
Y0.0500.1000.2000.3000.4000.5000.562
0.0502.497 2.486 *2.152 2.143 *1.835 1.827 *1.655 1.648 *1.527 1.520 *1.425 1.418 *1.370 1.363 *
0.1002.487 2.474 *2.143 2.133 *1.827 1.818 *1.648 1.640 *1.520 1.513 *1.419 1.412 *1.364 1.357 *
0.2002.440 2.426 *2.102 2.092 *1.792 1.784 *1.616 1.609 *1.491 1.485 *1.391 1.387 *1.337 1.334 *
0.3002.353 2.346 *2.028 2.023 *1.729 1.726 *1.559 1.558 *1.438 1.439 *1.343 1.346 *1.291 1.295 *
0.4002.234 2.235 *1.925 1.927 *1.642 1.645 *1.483 1.486 *1.369 1.375 *1.281 1.288 *1.233 1.241 *
0.5002.095 2.091 *1.807 1.804 *1.543 1.541 *1.396 1.394 *1.294 1.291 *1.216 1.213 *1.175 1.172 *
0.5622.003 1.986 *1.728 1.713 *1.478 1.464 *1.341 1.326 *1.247 1.230 *1.177 1.158 *1.142 1.121 *
Table 5. Comparison of SCF’s values calculated using Equations (A1)–(A3) (*) to the corresponding FEM results for various T/a ratios, while θ = 30°.
Table 5. Comparison of SCF’s values calculated using Equations (A1)–(A3) (*) to the corresponding FEM results for various T/a ratios, while θ = 30°.
θ = 30°KttKtbKts
ρ/at/aT/a = 1T/a = 2T/a = 3T/a = 4T/a = 1T/a = 2T/a = 3T/a = 4T/a = 1T/a = 2T/a = 3T/a = 4
0.05103.469
3.489 *
3.588
3.613 *
3.681
3.708 *
3.754
3.788 *
3.556
3.571 *
3.698
3.721 *
3.813
3.836 *
3.907
3.933 *
2.144
2.156 *
2.203
2.224 *
2.258
2.277 *
2.306
2.321 *
73.378
3.394 *
3.462
3.471 *
3.519
3.529 *
3.557
3.578 *
3.499
3.506 *
3.607
3.612 *
3.687
3.694 *
3.746
3.763 *
2.133
2.146 *
2.186
2.207 *
2.234
2.253 *
2.276
2.292 *
43.110
3.129 *
3.129
3.141 *
3.136
3.151 *
3.138
3.158 *
3.278
3.284 *
3.314
3.309 *
3.331
3.328 *
3.337
3.345 *
2.092
2.113 *
2.130
2.154 *
2.162
2.185 *
2.186
2.212 *
0.25102.334
2.347 *
2.413
2.426 *
2.474
2.487 *
2.522
2.538 *
2.391
2.400 *
2.485
2.498 *
2.561
2.574 *
2.624
2.638 *
1.705
1.714 *
1.752
1.768 *
1.793
1.809 *
1.833
1.844 *
72.273
2.283 *
2.328
2.331 *
2.366
2.368 *
2.390
2.399 *
2.352
2.356 *
2.424
2.425 *
2.477
2.478 *
2.515
2.523 *
1.696
1.706 *
1.738
1.754 *
1.775
1.791 *
1.808
1.822 *
42.094
2.107 *
2.106
2.114 *
2.110
2.120 *
2.111
2.124 *
2.204
2.208 *
2.227
2.224 *
2.238
2.236 *
2.242
2.246 *
1.663
1.680 *
1.694
1.712 *
1.718
1.737 *
1.738
1.758 *
0.5101.976
1.987 *
2.040
2.051 *
2.090
2.101 *
2.129
2.142 *
2.022
2.028 *
2.098
2.109 *
2.161
2.172 *
2.214
2.225 *
1.546
1.552 *
1.588
1.600 *
1.626
1.637 *
1.661
1.668 *
71.925
1.934 *
1.969
1.972 *
1.999
2.001 *
2.019
2.026 *
1.988
1.991 *
2.047
2.047 *
2.090
2.091 *
2.122
2.127 *
1.538
1.545 *
1.575
1.588 *
1.609
1.621 *
1.639
1.648 *
41.775
1.788 *
1.784
1.793 *
1.787
1.797 *
1.788
1.801 *
1.863
1.866 *
1.881
1.878 *
1.890
1.888 *
1.892
1.896 *
1.508
1.522 *
1.535
1.551 *
1.557
1.572 *
1.575
1.591 *
1101.687
1.696 *
1.737
1.748 *
1.776
1.788 *
1.807
1.821 *
1.723
1.728 *
1.783
1.796 *
1.832
1.848 *
1.874
1.892 *
1.405
1.409 *
1.442
1.452 *
1.475
1.484 *
1.506
1.512 *
71.643
1.652 *
1.677
1.682 *
1.700
1.705 *
1.715
1.725 *
1.693
1.695 *
1.738
1.741 *
1.772
1.777 *
1.796
1.807 *
1.398
1.403 *
1.430
1.441 *
1.460
1.469 *
1.486
1.494 *
41.519
1.534 *
1.525
1.537 *
1.527
1.540 *
1.527
1.543 *
1.584
1.588 *
1.597
1.598 *
1.603
1.605 *
1.604
1.611 *
1.371
1.383 *
1.394
1.408 *
1.413
1.427 *
1.429
1.443 *
Table 6. Comparison of SCF’s values calculated using Equations (A1)–(A3) (*) to the corresponding FEM results for various T/a ratios, while θ = 45°.
Table 6. Comparison of SCF’s values calculated using Equations (A1)–(A3) (*) to the corresponding FEM results for various T/a ratios, while θ = 45°.
θ= 45°KttKtbKts
ρ/at/aT/a = 1T/a = 2T/a = 3T/a = 4T/a = 1T/a = 2T/a = 3T/a = 4T/a = 1T/a = 2T/a = 3T/a = 4
0.05103.940
3.938 *
4.200
4.200 *
4.404
4.402 *
4.565
4.572 *
4.095
4.097 *
4.409
4.429 *
4.670
4.684 *
4.889
4.898 *
2.514
2.502 *
2.632
2.630 *
2.735
2.728 *
2.824
2.811 *
73.865
3.853 *
4.078
4.048 *
4.229
4.197 *
4.331
4.323 *
4.081
4.074 *
4.364
4.353 *
4.582
4.568 *
4.749
4.748 *
2.504
2.490 *
2.615
2.608 *
2.709
2.698 *
2.789
2.774 *
43.609
3.588 *
3.701
3.650 *
3.740
3.697 *
3.754
3.736 *
3.941
3.921 *
4.107
4.053 *
4.203
4.155 *
4.256
4.240 *
2.464
2.451 *
2.554
2.540 *
2.623
2.608 *
2.677
2.665 *
0.25102.380
2.384 *
2.523
2.525 *
2.636
2.634 *
2.727
2.725 *
2.456
2.470 *
2.630
2.651 *
2.776
2.790 *
2.899
2.907 *
1.829
1.820 *
1.913
1.911 *
1.986
1.981 *
2.049
2.040 *
72.335
2.333 *
2.451
2.436 *
2.533
2.514 *
2.589
2.581 *
2.450
2.451 *
2.606
2.601 *
2.727
2.715 *
2.820
2.812 *
1.821
1.812 *
1.900
1.895 *
1.967
1.959 *
2.024
2.013 *
42.183
2.177 *
2.230
2.207 *
2.250
2.230*
2.257
2.249 *
2.358
2.351 *
2.411
2.447 *
2.492
2.468 *
2.520
2.511 *
1.792
1.784 *
1.855
1.846 *
1.905
1.894 *
1.943
1.935 *
0.5101.959
1.962 *
2.063
2.066 *
2.147
2.145 *
2.215
2.212 *
2.012
2.024 *
2.140
2.158 *
2.248
2.262 *
2.340
2.349 *
1.603
1.596 *
1.673
1.672 *
1.734
1.730 *
1.789
1.780 *
71.922
1.920 *
2.005
1.993 *
2.063
2.050 *
2.104
2.097 *
1.999
2.004 *
2.110
2.113 *
2.199
2.196 *
2.266
2.267 *
1.596
1.588 *
1.662
1.658 *
1.718
1.711 *
1.766
1.756 *
41.799
1.795 *
1.830
1.814 *
1.842
1.829 *
1.847
1.842 *
1.915
1.915 *
1.976
1.959 *
2.010
1.993 *
2.029
2.022 *
1.571
1.564 *
1.623
1.616 *
1.664
1.655 *
1.696
1.689 *
1101.655
1.656 *
1.726
1.730 *
1.783
1.787 *
1.829
1.835 *
1.690
1.699 *
1.777
1.798 *
1.852
1.874 *
1.916
1.937 *
1.420
1.414 *
1.475
1.475 *
1.524
1.522 *
1.569
1.562 *
71.623
1.620 *
1.676
1.670 *
1.715
1.709 *
1.741
1.741 *
1.673
1.677 *
1.747
1.754 *
1.805
1.813 *
1.851
1.863 *
1.414
1.408 *
1.465
1.463 *
1.510
1.506 *
1.549
1.542 *
41.520
1.517 *
1.536
1.528 *
1.543
1.536 *
1.544
1.543 *
1.596
1.593 *
1.627
1.621 *
1.646
1.642 *
1.655
1.659 *
1.391
1.387 *
1.431
1.428 *
1.462
1.458 *
1.488
1.484 *
Table 7. Comparison of SCF’s values calculated using Equations (A1)–(A3) (*) to the corresponding FEM results for various T/a ratios, while θ = 55°.
Table 7. Comparison of SCF’s values calculated using Equations (A1)–(A3) (*) to the corresponding FEM results for various T/a ratios, while θ = 55°.
θ = 55°KttKtbKts
ρ/at/aT/a = 1T/a = 2T/a = 3T/a = 4T/a = 1T/a = 2T/a = 3T/a = 4T/a = 1T/a = 2T/a = 3T/a = 4
0.05104.005
4.026 *
4.339
4.367 *
4.604
4.628 *
4.816
4.848 *
4.180
4.198 *
4.591
4.641 *
4.939
4.981 *
5.234
5.267 *
2.714
2.716 *
2.875
2.888 *
3.011
3.020 *
3.127
3.132 *
73.941
3.951 *
4.225
4.216 *
4.430
4.419 *
4.574
4.590 *
4.193
4.220 *
4.579
4.593 *
4.887
4.894 *
5.129
5.147 *
2.703
2.703 *
2.857
2.862 *
2.983
2.985 *
3.089
3.089 *
43.707
3.703 *
3.849
3.802 *
3.915
3.879 *
3.938
3.943 *
4.115
4.110 *
4.383
4.336 *
4.549
4.509 *
4.649
4.656 *
2.663
2.659 *
2.791
2.784 *
2.890
2.880 *
2.965
2.961 *
0.25102.347
2.359 *
2.512
2.529 *
2.646
2.659 *
2.753
2.769 *
2.426
2.437 *
2.631
2.658 *
2.807
2.828 *
2.958
2.972 *
1.878
1.879 *
1.984
1.993 *
2.074
2.082 *
2.154
2.156 *
72.309
2.315 *
2.448
2.444 *
2.548
2.543 *
2.618
2.627 *
2.423
2.429 *
2.614
2.621 *
2.767
2.767 *
2.889
2.891 *
1.870
1.870 *
1.971
1.976 *
2.055
2.058 *
2.126
2.127 *
42.175
2.175 *
2.239
2.219 *
2.268
2.254 *
2.278
2.283 *
2.358
2.356 *
2.484
2.459 *
2.562
2.537 *
2.608
2.603 *
1.842
1.840 *
1.926
1.923 *
1.991
1.986 *
2.040
2.040 *
0.5101.931
1.940 *
2.046
2.059 *
2.139
2.149 *
2.215
2.226 *
1.984
1.991 *
2.126
2.146 *
2.250
2.265 *
2.356
2.366 *
1.620
1.620 *
1.705
1.712 *
1.778
1.782 *
1.842
1.842 *
71.899
1.904 *
1.993
1.991 *
2.061
2.058 *
2.108
2.114 *
1.974
1.978 *
2.104
2.108 *
2.208
2.208 *
2.291
2.293 *
1.613
1.612 *
1.694
1.697 *
1.761
1.762 *
1.819
1.817 *
41.789
1.790 *
1.828
1.817 *
1.845
1.838 *
1.851
1.855 *
1.907
1.904 *
1.986
1.968 *
2.034
2.016 *
2.062
2.057 *
1.589
1.588 *
1.655
1.653 *
1.706
1.702 *
1.745
1.745 *
1101.642
1.648 *
1.716
1.726 *
1.777
1.786 *
1.826
1.836 *
1.676
1.683 *
1.768
1.787 *
1.848
1.868 *
1.917
1.935 *
1.424
1.422 *
1.486
1.491 *
1.541
1.544 *
1.590
1.589 *
71.612
1.616 *
1.670
1.670 *
1.711
1.711 *
1.740
1.746 *
1.661
1.664 *
1.740
1.747 *
1.805
1.811 *
1.856
1.865 *
1.418
1.416 *
1.476
1.479 *
1.526
1.527 *
1.570
1.568 *
41.516
1.520 *
1.534
1.532 *
1.542
1.542 *
1.545
1.550 *
1.587
1.586 *
1.627
1.618 *
1.651
1.643 *
1.663
1.664 *
1.396
1.396 *
1.441
1.442 *
1.478
1.478 *
1.507
1.508 *
Table 8. Influence of the weld angle θ and t/ρ ratio on SCF for L > 4t.
Table 8. Influence of the weld angle θ and t/ρ ratio on SCF for L > 4t.
θt/ρ
0.81.01.52.03.04.05.010.020.035.050.075.0100.0
Ktt30°1.161.191.271.331.441.531.611.892.242.572.813.103.33
45°1.161.191.261.331.441.541.631.962.412.873.213.654.01
60°1.161.191.261.331.441.541.631.962.442.933.323.854.28
Ktb30°1.161.201.301.371.501.601.681.992.352.702.953.273.51
45°1.161.201.291.371.521.651.762.182.723.253.654.174.58
60°1.161.201.291.371.531.661.792.272.933.614.134.835.39
Table 9. Accuracy of Equations (A1) and (A7)–(A10) for a T-joint subjected to tension, where: θ = 45°, h = hp, t = T, h/t = 0.75, 0.025 ≤ ρ/t ≤ 0.35, #—out of range.
Table 9. Accuracy of Equations (A1) and (A7)–(A10) for a T-joint subjected to tension, where: θ = 45°, h = hp, t = T, h/t = 0.75, 0.025 ≤ ρ/t ≤ 0.35, #—out of range.
ρ/t0.0250.0500.0750.100.150.250.35
KttFEM3.0022.4172.1391.9671.7571.5441.430
U and N (A7)2.8852.2021.9231.7661.5881.4221.339
Tsuji (A8)2.4262.0471.8741.7691.6411.5111.440
Monahan (A9)2.8942.383# 2.150# 2.009# 1.840# 1.666# 1.571
Brennan et al. (A10)2.7492.270# 2.055# 1.925# 1.769# 1.611# 1.525
(A1)3.0152.4282.1481.9751.7641.5491.433
Table 10. Accuracy of Equations (A2) and (A11)–(A14) for the T-joint subjected to bending, where: θ = 45°, h = hp, t = T, h/t = 0.75, 0.025 ≤ ρ/t ≤ 0.35, #—out of range.
Table 10. Accuracy of Equations (A2) and (A11)–(A14) for the T-joint subjected to bending, where: θ = 45°, h = hp, t = T, h/t = 0.75, 0.025 ≤ ρ/t ≤ 0.35, #—out of range.
ρ/t0.0250.0500.0750.100.150.250.35
KtbFEM3.4342.7402.4042.1951.9351.6631.516
U and N (A11)3.3272.7502.4442.2381.9571.6251.436
Niu and Glinka (A12)# 3.516# 2.818# 2.503# 2.313# 2.086# 1.855# 1.730
Tsuji (A13)3.4712.8332.5392.3592.1421.9161.792
Brennan et al. (A14)3.2172.606# 2.330# 2.164# 1.964# 1.761# 1.652
(A2)3.4302.7412.4082.1991.9401.6691.523
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Molski, K.L.; Tarasiuk, P. Stress Concentration Factors for Welded Plate T-Joints Subjected to Tensile, Bending and Shearing Loads. Materials 2021, 14, 546. https://doi.org/10.3390/ma14030546

AMA Style

Molski KL, Tarasiuk P. Stress Concentration Factors for Welded Plate T-Joints Subjected to Tensile, Bending and Shearing Loads. Materials. 2021; 14(3):546. https://doi.org/10.3390/ma14030546

Chicago/Turabian Style

Molski, Krzysztof L., and Piotr Tarasiuk. 2021. "Stress Concentration Factors for Welded Plate T-Joints Subjected to Tensile, Bending and Shearing Loads" Materials 14, no. 3: 546. https://doi.org/10.3390/ma14030546

APA Style

Molski, K. L., & Tarasiuk, P. (2021). Stress Concentration Factors for Welded Plate T-Joints Subjected to Tensile, Bending and Shearing Loads. Materials, 14(3), 546. https://doi.org/10.3390/ma14030546

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