Microsoft Word - numero_55_art_18_2951 M. M. Konieczny et alii, Frattura ed Integrità Strutturale, 55 (2021) 241-257; DOI: 10.3221/IGF-ESIS.55.18 241 Research of stress distribution in the cross-section of a bimetallic perforated plate perpendicularly loaded with concentrated force Mateusz Marcin Konieczny, Henryk Achtelik, Grzegorz Gasiak Opole University of Technology, Poland mateuszmarcinkonieczny@wp.pl, kmpkm@po.edu.pl, g.gasiak@po.edu.pl ABSTRACT. The paper presents the stress distribution along the plate thickness in a bimetallic steel – titanium circular perforated plate produced in the technological process of explosion welding. The steel layer is the layer that transfers the load in the plate, while the titanium layer is used to improve the properties of the plate, e.g. corrosion resistance, thermal transmittance, etc. in the plate. Two cases of fastening were considered, i.e. a freely supported and fixed plate. Such plates are used in various engineering structures, e.g. simply supported plates can be used in loose material screens, while plates are fixed in heat exchangers. The load was assumed as a concentrated force applied perpendicularly to the plate surface. The results obtained numerically using the finite element method were compared with the results calculated according to the analytical equations. It has been shown that the difference in the results of equivalent von Mises stress calculations does not exceed 13%. The research results presented in the paper can be used by engineers to design bimetallic perforated plates perpendicularly loaded to their surface. KEYWORDS. Bimetallic perforated plate; Concentrated force; Stress analysis; Analytical calculations; Numerical calculations Citation: Konieczny, M. M., Achtelik, H., Gasiak, G., Research of stress distribution in the cross-section of a bimetallic perforated plate perpendicularly loaded with concentrated force, Frattura ed Integrità Strutturale, 55 (2021) 241-257. Received: 13.11.2020 Accepted: 17.12.2020 Published: 01.01.2021 Copyright: © 2019 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. INTRODUCTION n engineering structures, structural elements consisting of many layers [1, 2] and of many materials [3] are often used plated elements called metallic composites [4, 5]. Metallic composites (plating) are obtained by so-called explosive welding. Explosive welding, or bonding, makes it possible to join together many groups of construction materials by means of the energy released during the detonation of the explosive [6]. The unquestionable advantage of the explosive technology is the possibility of joining together materials that currently cannot be permanently joined by conventional methods such as: soldering, gluing or welding. However, the key benefit of using metal composites is the ability to significantly reduce the cost of making a given structure by reducing the amount of expensive materials with unique properties. An example of this type of material is bimetal steel - titanium, which, when used in engineering structures, maintains the assumed strength and anti-corrosion parameters, ensuring work safety, while reducing the amount of expensive titanium used. I https://youtu.be/K24shMRyUyw M. M. Konieczny et alii, Frattura ed Integrità Strutturale, 55 (2021) 241-257; DOI: 10.3221/IGF-ESIS.55.18 242 The main applications of steel – titanium plated materials include high-pressure tanks and highly specialized elements of process equipment, e.g. heat exchanger tube sheets (Fig. 1). Elements made of bimetal steel – titanium are most often used due to the improvement of physical properties, thermal properties, chemical properties (e.g. corrosion resistance) as well as mechanical parameters of a given structure in which such elements are present. In heat exchangers, an important issue is the analysis of the complex effect of the perforated plate on the connected shell of the chemical apparatus and the associated stress distribution in the shell, in the tubes and in the perforated plate itself [7, 8, 9]. Contemporary engineering structures, compared to their predecessors, are characterized by a reduced safety factor. Therefore, the precise calculation of temporary strength, stiffness and stability, in particular of thin-walled structural elements, is of increasing importance. Literature studies indicate the widening perspective of the use of plated materials [10, 11, 12, 13]. The paper proposes an analytical and numerical approach to determining the value of stresses along the plate thickness in a bimetallic circular axially symmetrical perforated plate, freely supported at the edge and fixed at the edge and centrally loaded with a concentrated force P perpendicular to the plate surface. The research conducted at work was carried out on the basis of the finite element method (FEM) and the analytical method according to mathematical formulas. A novelty in this work is the development and application of the analytical method to determine the stress state in the transition zone between the titanium layer and the steel layer in a bimetallic perforated plate loaded with e.g. a concentrated force. Figure 1: Model of bimetallic perforated plate. Figure 2: Model of bimetallic perforated plate consisting of steel thickness perforated plate H and titanium perforated plate thickness a containing holes with diameters d1,…,d5 and sections T1 and T2. M. M. Konieczny et alii, Frattura ed Integrità Strutturale, 55 (2021) 241-257; DOI: 10.3221/IGF-ESIS.55.18 243 MATERIAL AND GEOMETRY he calculation was based on a bimetallic circular axisymmetric perforated plate with dimensions: diameter D = 300 mm and comprising two layers, i.e. the base layer B in the form of structural steel with a thickness of H = 10 mm and the applied layer (plating) in the form of titanium with the thickness of a = 2.5 mm. 100 holes with different radii located on the plate. These holes were arranged in five circles with 20 holes in each circle. On the first outer circle, the plate had holes d1 = 20.5 mm in diameter and on the fifth, inner circle holes d5 = 9.5 mm in diameter (Fig. 2). Steel plate grade S355J2 was adopted as the base material, and titanium sheet was adapter as the applied material with the following mechanical and material parameters (Tab. 1) [14]. The materials used in the work, titanium and steel included in the bimetallic perforated plate, were modeled as elastic materials. S355J2 steel and titanium Applied layer A - titanium Re [MPa] Rm [MPa] E [GPa] v [-] G [GPa] A5 [%] 189-215 308-324 100 0.37 36.5 43-56 C Fe H N O Ti 0.10 0.20 0.015 0.03 0.18 99.5 Base layer B - S355J2 steel Re [MPa] Rm [MPa] E [GPa] v [-] G [GPa] A5 [%] 382-395 598-605 220 0,30 84,6 24-34 C Si Mn P S Cu Fe 0.22 0.55 1.60 0.025 0.025 0.45 all Table 1: Strength properties and chemical composition of S355J2 steel and titanium, where: Re – yield strength [MPa], Rm – tensile strength [MPa], E – Young’s modulus [GPa], v – Poisson’ s ratio [-], A5 – tensile elongation [%]. ANALYTICAL METHOD nlike the classical form of plate bending theory [15], in which the analysis is concerned with the state of equilibrium of internal forces acting on the separated full plate element, in this case the state of equilibrium of the plate element containing the halves of the holes is considered (Fig. 2) [16] and described in the cylindrical coordinate axis system r, θ, z. The following variables can be found in Fig. 3: 1m , 2m – circumferential torque intensity in the steel and titanium parts, respectively; 1r m , 2r m – radial torque intensity in the steel and titanium parts, respectively;  1 ,  2 – circumferential stress in the steel and titanium parts;  1r ,  2r – radial stress in the steel and titanium parts; E1 , E2, 1v , 2v – Young's modules and Poisson's ratios for steel and titanium parts respectively; h1, h2 – values determining the position of the inert layer in the cross-section of the plated plate; h – total thickness of the plate; t(r) – intensity of transverse force acting in cross-sections of the plated plate. The conditions of equilibrium of internal forces presented in Fig. 3 result in the following form of a differential equation [15] describing the state of plate effort: T U M. M. Konieczny et alii, Frattura ed Integrità Strutturale, 55 (2021) 241-257; DOI: 10.3221/IGF-ESIS.55.18 244      r r dm m m t r dr r (1) The intensity of moments mr1, mr2, mθ1, mθ2 in a bimetallic perforated plate are expressed by the respective relations [16]:         1 2 * ** r r r r d m m m C B B dr r (2)             1 2 ** *dm m m C B B dr r (3) where:  – is the angle of inclination tangential to the curved central surface of the plate.  * 1 2B B B  ** 1 1 2 2B v B v B       3 31 1 1 22 13 1 E B h h v ;         3 32 2 2 22 23 1 E B h a h v ; – bending strength of the steel and titanium layer, respectively. Figure 3: Wedge element cut from a bimetallic perforated plate in a state of static equilibrium. The rC , C coefficients determine the degree of weakening of the cross-section due to the existence of discontinuities in the form of perforations (Fig. 2) in the radial and circumferential direction. They were adopted in the respective form [16]:     Δ Δ r r d C r ;    Δ Δ r d C r (4) where: Δr , Δ – finite increments of radial coordinate r and circumferential one θ; d - diameter of the perforation. After introducing compounds (2), (3) into Eqn. (1), the following form was obtained: the equilibrium of internal forces acting on the element of the bimetallic circular perforated plate (Fig. 3): M. M. Konieczny et alii, Frattura ed Integrità Strutturale, 55 (2021) 241-257; DOI: 10.3221/IGF-ESIS.55.18 245                   2 ** 2 * 2 * 1 1 1 r r r t rC Cd B d C r dr Cdr B r C B (5) where the function of the transverse force intensity t(r) for the most common methods of plate support was assumed to take the form [16]:      221 3 k t r k r k r r (6) where: k1, k2, k3 – lateral force function coefficients. The solution to Eqn. (5) is given by the function of the plate deflection angle  (r):                            1 2 43 λ λ31 2 1 2* 1 2 1 2 1 2 1 3 λ 3 λ 1 λ 1 λ 4 λ 4 λr k rk r k r r D r D r C B (7) in which 1λ and 2λ are roots of the characteristic equation:          ** 2 * λ 1 λ 0 r r C CB C CB (8) In function (7), which is the general integral of the differential Eqn. (5), the constant coefficients k1, k2, k3 depend on the load method, while D1 and D2 are integration constants whose numerical values are determined on the basis of the boundary conditions given in this study. After determining the function φ(r) from compounds (2) and (3), the intensities of radial and circumferential moments in the base material (steel) mr1, mθ1 and in the plating layer (titanium) mr2, mθ2 are determined. The deflection of the bimetallic perforated plate with the plating layer is determined by the relationship:       3w r r dr D (9) where D3 – integration constant. The stress in the cross-sections of the bimetalic perforated plate (Fig. 3) is derived on the basis of the compounds: 1) Radial stress  r: a) in the base material of the plate (steel)  1r :      1 1 3 3 1 2 3 r r r m z C h h , dla   1 2h z h (10) b) in the applied layer of the plate (titanium)  2r :         2 2 3 3 1 2 3 r r r m z C h a h , dla    2 2h z h a (11) 2) Circumferential stress   : a) in the base material of the plate (steel)  1 :         1 1 3 3 1 2 3m z C h h , dla   1 2h z h (12) M. M. Konieczny et alii, Frattura ed Integrità Strutturale, 55 (2021) 241-257; DOI: 10.3221/IGF-ESIS.55.18 246 b) in the applied layer of the plate (titanium)  2 :            2 2 3 3 2 2 3m z C h a h , dla    2 2h z h a (13) The equivalent von Mises stress was determined from the relationship: a) for the steel layer:        1 1 1 1 1 2 2 red r r (14) b) for the titanium layer:        2 2 2 2 2 2 2 red r r (15) In the formulas for determining the intensity of bending moments in the radial direction (2) and circumferential direction (3) as well as in the formulae applied for determining the stresses in the radial direction (10) and (11) and circumferential directions (12) and (13) h1 and h2 determine the position of the layer inert sections of the plate (Fig. 3):      2 2 1 2 2 1 22 E H E a h E H E a (16)  1 2h H h (17) where: h1, h2 – values determining the state of the inactive layer in the cross section of the steel plate; E1 , E2 – Young's modules for the steel and titanium parts; H, a – thickness of steel and titanium plate, respectively. Inactive layer of the plate – it is a layer in which the deformations are equal to zero [15, 16]. Bimetallic perforated plate freely supported on the perimeter and loaded with a concentrated force P acting perpendicular to the surface In the first case, the analysis the state of stress in a bimetallic perforated plate loaded with a concentrated force applied perpendicularly to the surface with the value P = 10 kN was presented below. The geometry of the plate is presented in Fig. 2. The concentrated force P was applied centrally to the pressure stamp with a diameter of b1 = 14 mm. A free support of the plate on the outer perimeter was assumed. The method of support and loading the plate is shown in Fig. 4. The boundary conditions are assumed to take the following form (Fig. 4): a) for the radius r – ½ b = 6 mm – plate deflection angle φ = 0; b) for the radius r = ½ D = 150 mm – radial moment strength  0rm . In this case, the coefficients k1, k2, k3 inside the function of the lateral force intensity t(r) (6) are equal to, respectively: k1 = 0; k2 = 2 P ; k3 = 0 The roots of the Eqn. (8) 1λ and 2λ are:   1λ 1.127  ;    2λ 1.065  M. M. Konieczny et alii, Frattura ed Integrità Strutturale, 55 (2021) 241-257; DOI: 10.3221/IGF-ESIS.55.18 247 Integration constants D1 and D2 in the Eqn. (7) are equal to:    71 21.21545 * 10  D k mm ;    52 25.1032 * 10  D k mm ; Cross-section weakening coefficients (4) C and rC in the analyzed plate are respectively:    0.590 C ;   0.496 rC ; Figure 4. Bimetallic perforated plate freely supported on the perimeter and loaded with a concentrated force P acting perpendicular to the surface. For calculations it was assumed: P = 10 kN – load in the form of concentrated force distributed around the circumference of the hole with a diameter; b = 12 mm – hole diameter; b1 – 14 mm – pressure stamp diameter; D = 300 mm – outer diameter of the plate; H = 10 mm – base layer thickness (steel plate); a = 2.5 mm – applied layer thickness (titanium plate); h1 = 5.64 mm (Fig. 2); h2 = 4.36 mm (Fig. 2); The mechanical properties of the steel plate and titanium plate are taken from Tab. 1. On the basis of the Eqns. (10) – (15), values of radial stresses  ar and circumferential  a and equivalent von Mises stresses  ared were determined for the adopted geometry of the bimetallic perforated plate given in Fig. 4. For example, in Fig. 5 the values of the circumferential stress  a , radial stress  ar and equivalent von Mises stresses  a red are plotted along the thickness h of the bimetallic perforated plate in sections T1 and T2, i.e. on the radius r = R5 = 60 mm and r = 130 mm (Fig 2), loaded from the side of the titanium layer with a concentrated force applied perpendicularly to the center of the plate P = 10 kN. Bimetallic perforated plate fixed supported on the perimeter and loaded with a concentrated force P acting perpendicular to the surface In the second case, the analysis the state of stress in a bimetallic perforated plate loaded with a concentrated force applied perpendicularly to the surface with the value P = 10 kN was presented below. The geometry of the plate is presented in Fig. 2. The concentrated force P was applied centrally to the pressure stamp with a diameter of b1 = 14 mm. A fixed of the plate on the outer perimeter was assumed. The method of support and loading the plate is shown in Fig. 6. The boundary conditions are assumed to take the following form (Fig. 6): a) for the radius r – ½ b = 6 mm – plate deflection angle φ = 0; b) for the radius r = ½ D = 150 mm – plate deflection angle φ = 0. M. M. Konieczny et alii, Frattura ed Integrità Strutturale, 55 (2021) 241-257; DOI: 10.3221/IGF-ESIS.55.18 248 Figure 5. Stress distribution, respectively:  a ,  ar , and  a red in the cross-section of a bimetallic perforated plate freely supported and loaded with a concentrated force with the value P = 10 kN: a) in the section T1 on the radius r = 5R = 60 mm; b) in the section T2 on radius r = 130 mm (Fig. 2). Figure 6. Bimetallic perforated plate fixed on the perimeter and loaded with a concentrated force P acting perpendicular to the surface In this case, the coefficients k1, k2, k3 inside the function of the lateral force intensity t(r) (6) are equal to, respectively: k1 = 0; k2 = 2 P ; k3 = 0 The roots of the Eqn. (8) 1λ and 2λ are:   1λ 1.127  ;    2λ 1.065  Integration constants D1 and D2 in the Eqn. (7) are equal to:    71 21.322 * 10  D k mm ;    52 23.219 * 10  D k mm M. M. Konieczny et alii, Frattura ed Integrità Strutturale, 55 (2021) 241-257; DOI: 10.3221/IGF-ESIS.55.18 249 Cross-section weakening coefficients (4) C and rC in the analyzed plate are respectively:    0.590 C ;   0.496 rC ; Data as in the above case were adopted for the calculations. On the other hand, the mechanical properties of the steel plate and the titanium plate were taken from Tab. 1. On the basis of the Eqns. (10) – (15), values of radial stresses  ar and circumferential  a and equivalent von Mises stresses  ared were determined for the adopted geometry of the bimetallic perforated plate given in Fig. 6. For example, in Fig. 7 the values of the circumferential stress  a , radial stress  ar and equivalent von Mises stresses  a red are plotted along the thickness h of the bimetallic perforated plate in sections T1 and T2, i.e. on the radius r = R5 = 60 mm and r = 130 mm (Fig 2), loaded from the side of the titanium layer with a concentrated force applied perpendicularly to the center of the plate P = 10 kN. Figure 7. Stress distribution, respectively:  a ,  ar , and  a red in the cross-section of a bimetallic perforated plate fixed and loaded with a concentrated force with the value P = 10 kN: a) in the section T1 on the radius r = 5R = 60 mm; b) in the section T2 on radius r = 130 mm (Fig. 2). NUMERICAL METHOD he subject of numerical calculations was the stress distribution in a bimetallic circular perforated plate, in the first case freely supported (Fig. 4) and in the second case fixed (Fig. 6) and loaded with a concentrated force applied perpendicular to the center of the plate. The free support of the plate was realized by designing a ring support on which a bimetallic perforated plate was placed. Minimal friction contact is applied between the plate and the support. On the other hand, the plate was fixed by removing all degrees of freedom from the finite elements located on the outer contour of the perforated plate. T M. M. Konieczny et alii, Frattura ed Integrità Strutturale, 55 (2021) 241-257; DOI: 10.3221/IGF-ESIS.55.18 250 The three-dimensional 3D computational model of the bimetallic circular axially-symmetrical perforated plate was created and discretized by the possibility of modeling and creating a finite element mesh in ANSYS [17]. The finite elements used for discretization were SOLID186 finite elements. The calculation model contained the total number of 3379707 finite elements, while the total number of nodes was 7027927. The entire plate had 25 layers of finite elements (20 layers of base plate and 5 layers of applied layer) (Fig. 8). Figure 8. Finite elements of a bimetallic circular perforated plate in ANSYS program. The base layer of the B – steel plate is combined with the applied layer of the A – titanium plate by using the "connections" function, providing the option to bond two or more elements [17]. It was assumed that the combination of the titanium layer and the steel layer of the plates is excellent. Modeling the connection of these two layers, contact was used in the form of inseparable bonding – as both layers were bonded permanently. Bimetallic perforated plate freely supported on the perimeter and loaded with a concentrated force P acting perpendicular to the surface For example, the distribution of equivalent von Mises stress in the cross section of the plate is presented in Fig. 9. The figure illustrates the distribution of equivalent von Mises stress in a bimetallic perforated plate loaded with a force of P = 10 kN in the applied layer A – titanium, in the applied transition zone A’- titanium, in the based transition zone B’ – steel and based layer B - steel determined for the radius r = R5 = 60 mm. The following numerical calculation results (Fig. 9) demonstrate the places where the equivalent von Mises stress values occur, i.e. point AP’ – the place of maximum equivalent M. M. Konieczny et alii, Frattura ed Integrità Strutturale, 55 (2021) 241-257; DOI: 10.3221/IGF-ESIS.55.18 251 von Mises stress near the hole in the titanium transition zone, BP’ – the value of the maximum equivalent von Mises stress near the hole in the transition zone of structural steel, AP– the place of maximum equivalent von Mises stress near the hole in the titanium layer, BP – the value of the maximum equivalent von Mises stress near the hole in the structural steel layer. The maximum equivalent von Mises stress around the hole and its coordinates in the vicinity of the hole were respectively: red maxσ 97.997 MPa in A layer, where AP [x = 6.0 mm, y = 68.51 mm, z = 0.0 mm], red maxσ 168.83 MPa in B layer, where BP [x = 6.0 mm, y = 68.81 mm, z = 12.5 mm], while red maxσ 66.417 MPa in A’ zone, where AP’ [x = 6.0 mm, y = 68.51 mm, z = 2.49 mm] and red maxσ 130.27 MPa in B’ zone, where BP’ [x = 6.0 mm, y = 68.51 mm, z = 2.51 mm]. Figure 9. Distribution of equivalent von Mises stress  nred in the cross section of the plate determined on the radius r = R5 = 60 mm. Whereas in Fig. 7 the values of the circumferential stress  a , radial stress  ar and equivalent von Mises stresses  a red are plotted along the thickness h of the bimetallic perforated plate in sections T1 and T2, i.e. on the radius r = R5 = 60 mm and r = 130 mm (Fig 2), loaded from the side of the titanium layer with a concentrated force applied perpendicularly to the center of the plate P = 10 kN. Bimetallic perforated plate fixed supported on the perimeter and loaded with a concentrated force P acting perpendicular to the surface For example, the distribution of equivalent von Mises stress in the cross section of the plate is presented in Fig. 11. The figure illustrates the distribution of equivalent von Mises stress in a bimetallic perforated plate loaded with a force of P = 10 kN in the applied layer A – titanium, in the applied transition zone A’- titanium, in the based transition zone B’ – steel and based layer B - steel determined for the radius r = R5 = 60 mm. The following numerical calculation results (Fig. 11) demonstrate the places where the equivalent von Mises stress values occur, i.e. point AP’ – the place of maximum equivalent von Mises stress near the hole in the titanium transition zone, BP’ – the value of the maximum equivalent von Mises stress near the hole in the transition zone of structural steel, AP– the place of maximum equivalent von Mises stress near the hole in the titanium layer, BP – the value of the maximum equivalent von Mises stress near the hole in the structural steel layer. The maximum equivalent von Mises stress around the hole and its coordinates in the vicinity of the hole were respectively: red maxσ 51.155 MPa in A layer, where AP [x = 8.88 mm, y = -54.81 mm, z = 0.0 mm], red maxσ 91.816 MPa in B layer, where BP [x = 8.88 mm, y = -54.81 mm, z = 12.5 mm], while red maxσ 32.258 MPa in A’ zone, where AP’ [x = 8.88 mm, y = -54,81 mm, z = 2.49 mm] and red maxσ 71.527 MPa in B’ zone, where BP’ [x = 8.88 mm, y = -54.81 mm, z = 2.51 mm]. M. M. Konieczny et alii, Frattura ed Integrità Strutturale, 55 (2021) 241-257; DOI: 10.3221/IGF-ESIS.55.18 252 Figure 10. Stress distribution, respectively:  n ,  nr , and  n red in the cross-section of a bimetallic perforated plate freely supported and loaded with a concentrated force with the value P = 10 kN: a) in the section T1 on the radius r = 5R = 60 mm; b) in the section T2 on radius r = 130 mm (Fig. 2). Figure 11. Distribution of equivalent von Mises stress  nred in the cross section of the plate determined on the radius r = R5 = 60 mm. M. M. Konieczny et alii, Frattura ed Integrità Strutturale, 55 (2021) 241-257; DOI: 10.3221/IGF-ESIS.55.18 253 Whereas in Fig. 12 the values of the circumferential stress  a , radial stress  ar and equivalent von Mises stresses  a red are plotted along the thickness h of the bimetallic perforated plate in sections T1 and T2, i.e. on the radius r = R5 = 60 mm and r = 130 mm (Fig 2), loaded from the side of the titanium layer with a concentrated force applied perpendicularly to the center of the plate P = 10 kN. Figure 12. Stress distribution, respectively:  n ,  nr , and  n red in the cross-section of a bimetallic perforated plate fixed and loaded with a concentrated force with the value P = 10 kN: a) in the section T1 on the radius r = 5R = 60 mm; b) in the section T2 on radius r = 130 mm (Fig. 2). RESULTS OF THE RESEARCH he problem of the stress analysis on the plate thickness, taken up in the paper, enables the appropriate selection of the base layer of the steel plate as well as the layer of the applied plate – titanium. The base layer is the layer that carries the load in the plate, in this case it is the steel layer of the plate. However, in this case, the titanium layer of the plate is used due to the improvement of the plate properties, e.g. corrosion resistance, thermal permeability, etc. The thickness of the applied layers, depending on the type of use of the plated plate, may be from 1.5–15% of the thickness of the base plate. Taking into account the thickness of the layer of the applied plate during design also affects the strength of the bimetallic perforated plate. For example, Fig. 13 shows the values equivalent von Mises stresses obtained numerically by the finite element method  nred and analytically  ared in a bimetallic circular perforated plate freely supported on the edge and loaded centrally with a concentrated force P = 10 kN perpendicular to the plate surface and the percentage difference between the results obtained for these cases red . For example, Fig. 14 presents the values of equivalent von Mises stresses obtained numerically using the finite element method  nred and analytically  a red in a bimetallic circular perforated plate fixed at the edge and loaded centrally with a concentrated force P = 10 kN acting perpendicular to the plate surface and the percentage difference between the results obtained for these cases red . Apart from numerical calculations the state of stress using the finite element method, an attempt was made to develop analytical equations allowing to determine the stress values in a bimetallic perforated plate subjected to a normal load in the form of a concentrated force. The results the state of stress calculations obtained according to the analytical method were T M. M. Konieczny et alii, Frattura ed Integrità Strutturale, 55 (2021) 241-257; DOI: 10.3221/IGF-ESIS.55.18 254 compared with the results obtained by the finite element method according to the ANSYS program. The results of these comparisons are presented in Fig. 13 – in the case of freely supported plate and in Fig. 14 – in the case of fixed edges, respectively. Fig. 13 shows that the maximum difference in the stress values is equivalent von Mises stress   7.49%ared . On the other hand, Fig. 14 shows that the maximum difference in the stress values is equivalent von Mises stress   12.34%ared . The calculations carried out in the work with the analytical and numerical method using the ANSYS program show that the maximum difference in the results obtained with these methods in the case of equivalent von Mises stress is  ared = 12.34%. (around 13%), i.e.:    σ σ *100% σ a n a red red red a red where:  nred ,  n red –equivalent vo Mises stress determined by analytical and numerical methods, respectively. Figure 13: Comparison of equivalent von Mises stress  red determined analytically and numerically in the cross-section of the bimetallic perforated plate freely supported and loaded with a concentrated force P = 10 kN determined on the radius r = R5 = 60 mm. Figure 14: Comparison of equivalent von Mises stress  red determined analytically and numerically in the cross-section of the bimetallic perforated plate fixed and loaded with a concentrated force P = 10 kN determined on the radius r = R5 = 60 mm. M. M. Konieczny et alii, Frattura ed Integrità Strutturale, 55 (2021) 241-257; DOI: 10.3221/IGF-ESIS.55.18 255 A significant difference in the results obtained with the analytical method and the finite element method results from the fact that the bimetallic plate has a perforation, which causes a significant heterogeneity of the material, which was fully taken into account in the analytical solution. Therefore, it is believed that the error of about 13% is the discrepancy between the results obtained by the analytical method and the results obtained with the ANSYS program. The paper shows that the ANSYS program allows to determine the state of stress in a sufficiently accurate state in a bimetallic perforated plate at any point, i.e. on the bridge between the holes, around the holes, along the thickness of the titanium layer and the steel layer. It also enables the determination of stress concentration zones, and thus the location of dangerous places in the designed plate. Knowing the coordinates of stress concentration points, it is possible to predict the zones of crack initiation and propagation in bimetallic perforated plates. Despite significant differences in the stress values determined by the analytical method, it is believed that the analytical solution of the state of stress in a bimetallic perforated plate may be helpful for engineers to estimate the effort of the designed multilayer structures of the perforated plate type. It should be mentioned that the use of analytical solutions when designing bimetallic perforated plates allows to obtain the stress values in the structure relatively quickly. The paper shows that the method of supporting a bimetallic perforated plate has a significant impact on its state of stress. Tabs. 2 and 3 show that for the same load, the stresses in the plate freely supported along the perimeter are higher in relation to the plate fixed at the edge. The graphs in Figs. 5, 7, 10 and 12 show that both the analytical method and the ANSYS program, engineers can present the curves of the circumferential stress  , radial stress r , equivalent von Mises stress  red distributions along the plate thickness, both in the applied layer – titanium, the base layer – steel and in the transition zone between the overlay and the edge layer. The existence of a significant difference in the stress values in the transition zone between the applied layer A’ - titanium and the base layer B’ – steel. Figs. 9 and 11 show as an example the location of points along the plate thickness, where we plate with the maximum values of equivalent von Mises stresses obtained by the ANSYS program. The analytical method does not provide such possibilities. CONCLUSIONS he following conclusions can be drawn from the above analysis: 1) The attempt to use the analytical method to determine the state of stress in a bimetallic perforated plate resulted in the discrepancy of the calculation results in the equivalent von Mises stress values in comparison with the finite element method, up to approximately 13%; 2) It is proposed that the analytical method be recommended to design engineers as a tool for the estimation of the state of stress of multilayer perforated plates; 3) It was shown that in the transition zone of the plate between the applied layer – titanium and the base layer – steel, there is a step change in stress values. REFERENCES [1] Cheng Wu, X.J. (1999). A higher order theory for plane stress conditions of laminates consisting of isotropic layers, Journal of Applied Mechanics, 66, pp. 95-100. [2] Jafari, M., Mounssavian, H. and Chaleshtari, M.H. (2018). Optimum Design of Perforated Orthotropic and Laminated Composite Plates under in–plane Loading by Genetic Algorithm, Structural and Multidisciplinary Optimization, 57, pp. 341-357, DOI: 10.1007/s00158-017-1758-5. [3] Shaikh, E.N., Panchal, K.C. and Patel, D.B. (2015). Stress analysis of an infinite plate with different shaped cutouts in composite plate, International Journal of Science Technology and Management, 4 (1), pp. 307-317. [4] Kosturek, R., Wachowski, M., Śnieżek, L. and Gloc, M. (2018). The Influence of the Post Weld Heat Treatment on the Microstructure of Inconel 625/Carbon Steel Bimetal Joint Obtained by Explosive Welding, www.preprints.org., pp. 1- 15. [5] Konieczny, M., Achtelik, H. and Gasiak, G. (2020). Stress distribution in plated perforated plate loaded centrally with concentrated force, Zmęczenie materiału w eksploatacji maszyn roboczych, PO, 2, pp. 47-62 (in Polish). [6] Walczak, W. (1989). Metal explosion welding. WTN, Warsaw (in Polish). T M. M. Konieczny et alii, Frattura ed Integrità Strutturale, 55 (2021) 241-257; DOI: 10.3221/IGF-ESIS.55.18 256 [7] Konieczny, M., Achtelik, H. and Gasiak, G. (2020). Finite Element Analysis (FEA) and experimental stress analysis in circular perforated plates loaded with concentrated force, Frattura ed Integrità Strutturale, 51, pp. 164-173, DOI: 10.3221/IGF-ESIS.51.13. [8] Konieczny, M., Achtelik, H. and Gasiak, G. (2019). Research of maximum stresses zones in circular perforated plates made of S235JR steel loaded with concentrated force, Materials Engineering, 2 (228), pp. 39-46, DOI: 10.15199/28.2019.2.4. [9] Konieczny, M., Achtelik, H. and Gasiak, G. (2020). Analysis of Reduced Stress Distribution in Circular Perforated Plates, Fundamentals of Machine Design: Selected problems, Military University of Technology, pp. 133-145. [10] Kurek, A., Niesłony, A. and Szulc, Z. (2013). Design of process equipment made from explosively cladded materials, including the imposed material thickness in calculations, Przegląd Mechaniczny, 12, pp. 22-27. [11] Karolczuk, A., Kowalski, M., Bański, R. and Żok, F. (2013). Fatigue phenomena in explosively welded steel – titanium clad components subjected to push-pull loading, International Journal of Fatigue, 48, pp. 101-108, DOI: 10.1016/j.ijfatigue.2012.10.007 . [12] Karolczuk, A. and Kowalski, M. (2013). Residual stress determination based on the hole drilling method in explosively welded bimetallic composite, Third International Conference on Material Modelling, Warsaw, p. 140. [13] Konieczny, M., Achtelik, H. and Gasiak, G. (2021). Location of stress concentration zones in a two-layer axially symmetrical perforated plate with force applied normally to its surface, Engineering Structures, 226 , pp. 1-14, DOI: 10.1016/j.engstruct.2020.111355. [14] Data from Explomet Company research (Z.T.W. EXPLOMET, Gałka, Szulc, Sp. j. ul. Oświęcimska 100H, 45-641 Opole, Poland). [15] Timoshenko, S. and Woinowsky – Krieger, S. (1959). Theory of plates and shells, 2nd ed. New York, McGraw – Hill Book Company, pp. 62-85. [16] Achtelik, H., Gasiak, G. and Grzelak, J. (2008). Strength tests of axially symmetric perforated plates for chemical reactors: Part 1 - The simulation of stress state, International Journal of Pressure Vessels and Piping, 85, pp. 248-256, DOI: 10.1016/j.ijpvp.2007.08.009. [17] User’s Guide ANSYS 2020 R1, Ansys, Inc., USA. NOMENCLATURE H, a – thickness of steel and titanium plate, respectively; h – thickness of plate; D – diameter of the circular axially-symmetrical plated perforated plate; b – diameter of the hole in the central part of the plate; b1 – pressure stamp diameter; Cr, Cθ – coefficients of weakening of the plate cross-section due to the existence of discontinuities in the form of holes in the radial and circumferential directions, respectively; B*, B** - plate stiffness coefficients;    r – tilt angle of the tangent to the curved central surface of the plate; k1, k2, k3 – coefficients of the lateral force function; mθ1, mθ2 – circumferential torque intensity in the steel and titanium parts, respectively; mr1, mr2 – radial moment intensity in the steel and titanium parts, respectively; 1λ , 2λ – roots of the characteristic equation for the differential equation of the perforated plate balance; D1, D2, D3 – integration constants; , r – radial and peripheral coordinates of the plate;  1 ,  2 – circumferential stress in the steel and titanium parts, respectively;  1r ,  2r – radial stress in the steel and titanium parts, respectively; E1 , E2, 1v , 2v – Young's modules and Poisson's coefficients for steel and titanium parts; h1, h2 – values determining the position of the inactive layer in the cross-section of the plated plate; 1B , 2B – bending stiffness of the steel and titanium layer, respectively; t(r) – intensity of the transverse force acting in the cross-sections of the plated plate; w = w(r) – deflection of the perforated plate; M. M. Konieczny et alii, Frattura ed Integrità Strutturale, 55 (2021) 241-257; DOI: 10.3221/IGF-ESIS.55.18 257  1red ,  2red – equivalent von Mises stress in steel and titanium parts, respectively; Rn – radius of the circles on which the holes on the plate are arranged; dn – diameter of the holes; P – load in the form of concentrated force; r – radius of the circles on which the measuring points are located; T1, T2 – sections; B – base (steel) layer of the plate; B’ – basic transition zone (from the steel side) of the plate; A’ – transition zone applied (from the titanium side) to the plate; A – applied (titanium) layer of the plate.