IPASJ International Journal of Mechanical Engineering (IIJME) Web Site: http://www.ipasj.org/IIJME/IIJME.htm Email: [email protected] ISSN 2321-6441 A Publisher for Research Motivation........ Volume 3, Issue 4, April 2015 Numerical modeling and analytical validation of stress and stress intensity factor for SENT tensile specimen of P265GH steel material M.LAHLOU1, M.RACHIK2, M.BECHTAOUI3, A.HACHIM4, M. EL GHORBA5 1 Laboratory of Control and Mechanical Characterization of Materials and Structures, National Higher School of Electricity and Mechanics, BP 8118 Oasis, Hassan II University, Casablanca, Morocco 2 Roberval Laboratory. UTC, Royallieu Research Center. Person street from Roberval, BP20520,60205 Compiègne Cedex, France 3 ISEM/Higher Institute of Maritims Studies, Laboratory of Mechanics, Km 7 Road El Jadida Casablanca, Morocco 4 ISEM/Higher Institute of Maritims Studies, Laboratory of Mechanics, Km 7 Road El Jadida Casablanca, Morocco 5 Laboratory of Control and Mechanical Characterization of Materials and Structures, National Higher School of Electricity and Mechanics, BP 8118 Oasis, Hassan II University, Casablanca, Morocco ABSTRACT In a material subjected to a mechanical loading, the stress distribution is not uniform in the vicinity of a defect, causing a stress concentration in this area which result a sudden fracture. Therefore, the effect of the notch is generally reflected by an increase in stress at the bottom of the notch. The objective of this paper is to establish a numerical finite element modeling for a SENT tensile specimen (Single Edge Notch Tension) using CASTEM2013 computer code. The studied material is P265GH steel commonly used in sheet form in boilers and pressure vessels. The results show that the stresses exhibit a parabolic trend with a maximum value at the bottom of the notch at first, and then stress stabilization is noticed at the nominal stress value. The length of the critical notch decreases with increasing stress. Keywords: Notch, Finite element model, Stress, Stress intensity factor, Pressure vessels. 1. INTRODUCTION Through the development of information technology, new tool of production have been used in many industrial sectors that become currently inescapable: numerical finite element modeling; this method is a numerical way to resolve the problems of mechanical systems which allows determining an approximate solution on a spatial domain [1]. Nevertheless, in metallic structures, cracks are mostly initiated at geometric discontinuities of notches or defects. The geometric parameters and discontinuities govern cracks initiation or propagation and therefore affect the resistance of structures during their use [2]. In industry, for economic or security reasons it is seek to know the degree of defects harmfulness and residual life time of structures; This requires the development of models based on fracture mechanics. The behavior simulation with FEM has been presented by many authors with the aim to improve the knowledge of predict trends. A. HACHIM [3,4] presented a finite element based approach to simulate a Double Edge Notch Tension specimen of S355 Steel; he studied the behavior of the material in the presence of defects. Y.HIROSHI [5] presented the critical stress intensity factor on SENT specimen. A. EL Hakimi [6] studied the correction function i0, by applying a constant pressure along the lips of crack. According to the results of the calculations, the integral J in the case of a defect at the base of the transition is always higher than that of the similar case in a straight tube [7]. In the present article, the P265GH steel used in pressure vessels, and its current specifications, are introduced first. The tensile properties for selected steel are then evaluated. A numerical finite element modeling based on the tensile properties has been developed to determine the life prediction. Finally, a parametric study of tenacity was performed. 2. EXPERIMENTATION To extract the mechanical characteristics of the P265GH steel used in our program, tensile tests of standard specimens (Figure 1a and 1b) were conducted in different directions of rolling (longitudinal and transversal). The test curves showing the stress versus elongation are given in Figure 2: Volume 3, Issue 4, April 2015 Page 42 IPASJ International Journal of Mechanical Engineering (IIJME) Web Site: http://www.ipasj.org/IIJME/IIJME.htm Email: [email protected] ISSN 2321-6441 A Publisher for Research Motivation........ Volume 3, Issue 4, April 2015 Figure 1-a : Image of the standard test specimen Figure 1-b : Dimensions of the standard test specimen Figure 2: Stress strain curve test By comparing the mechanical characteristics of the specimen in both rolling directions, it is found that there is a negligible difference between the two curves. The mechanical characteristics of P265GH steel, at the ambient temperature, are reported in the table1: Table1. Mechanical properties of the material Young's modulus E (MPa elastic limit: e (Pa) 2.105 320 Breaking stress: g (Pa) 470 Elongation % Poisson's ratio ν 35 0,3 We notice that the elongation is about 35% which is higher than 14% required by the CODAP[8]. Therefore, this P265GH steel used is well adapted for pressurized structures. Volume 3, Issue 4, April 2015 Page 43 IPASJ International Journal of Mechanical Engineering (IIJME) Web Site: http://www.ipasj.org/IIJME/IIJME.htm Email: [email protected] ISSN 2321-6441 A Publisher for Research Motivation........ Volume 3, Issue 4, April 2015 3. NUMERICAL MODELING The calculation code Cast3m 2013 is used to construct a finite element model to analyze SENT specimen behavior subjected to tensile stress. In what follows, we describe FE modeling. 3.1Geometry The geometry and dimensions of the studied specimen are shown in Figure 3. Since the study is restricted to the I, The specimen was subjected to tensile solicitations. mode Figure 3. Geometry and dimensions of the studied specimen 3.2Mesh and boundary conditions By taking into consideration the symmetry of the problem, only half of the test specimen is discretized. Because the numerical results are intended for analysis of fracture mechanics, special attention is paid to mesh principally in crack and its vicinity (Mesh Refinement using Barsoum elements)[9]. Details of the mesh are illustrated in Figure 4a and 4b. x Figure 4-a : Mesh Figure 4-b :Mesh in the vicinity of a notch 3.3Loading The simulated load is a tensile solicitation along the longitudinal axis of the specimen. To avoid bending or twisting parasite and to ensure that the tensile stress is perfectly aligned; it is applied on the specimen via a rigid triangle indicated by the arrow 4a. The selected loads are calibrated in such a way that the applied nominal stresses are respectively 148 MPa, 284MPa and 356MPa. 4. RESULTS AND DISCUSSION 4.1Evolution of the stress in the ligament The curves in Figure5 illustrate the evolution of the stress in the ligament of the specimen, ie, along the x axis of Figure 4a (origin at the crack bottom) for the three levels of applied stress: Δσ = 148 MPa, 284MPa, 356MPa Volume 3, Issue 4, April 2015 Page 44 IPASJ International Journal of Mechanical Engineering (IIJME) A Publisher for Research Motivation........ Volume 3, Issue 4, April 2015 Web Site: http://www.ipasj.org/IIJME/IIJME.htm Email: [email protected] ISSN 2321-6441 Figure 5: Evolution of the Von Mises stress in the ligament of the specimen for the three levels of applied stress (Δσ = 356MPa, 284MPa, 148MPa). For Δσ = 148 MPa, there are three zones: The first zone [0, 0.5mm] (in the vicinity of the notch bottom) corresponds to the interval of stress values [470, 653MPa] .The maximum stress is 653 MPa which exceeds the strength of the material. The stress then decreases until reaching the breaking stress σg = 470MPa. The second zone [0.5, 2mm] matching the stress values that are included in the interval [320, 470MPa]. The stress is greater than the elastic limit of the material; this zone is the seat of plastic deformations. The third zone [2, 25mm] that corresponds to the interval of stress values [148, 320MPa]; We are observing a slight diminution of stress to stabilize at the applied stress value which is below the elastic limit. Thus, this zone is elastically deformed. Volume 3, Issue 4, April 2015 Page 45 IPASJ International Journal of Mechanical Engineering (IIJME) A Publisher for Research Motivation........ Volume 3, Issue 4, April 2015 Web Site: http://www.ipasj.org/IIJME/IIJME.htm Email: [email protected] ISSN 2321-6441 For Δσ = 284 MPa, there are three zones: The first zone [0, 1.7mm] (in the vicinity of the notch bottom) corresponds to the interval of stress values [470, 1241MPa] .The maximum stress is 1241 MPa which exceeds the strength of the material. The stress then decreases until reaching the breaking stress σg = 470MPa. The second zone [1.7, 4.5mm] matching the stress values that are included in the interval [320, 470MPa]. The stress is greater than the elastic limit of the material; this zone is the seat of plastic deformations. The third zone [4.5, 25mm] that corresponds to the interval of stress values [284, 320MPa]; We are observing a slight diminution of stress to stabilize at the applied stress value which is below the elastic limit. Thus, this zone is elastically deformed. For Δσ = 356 MPa, there are two zones: The first zone [0, 3.9mm] (in the vicinity of the notch bottom) corresponds to the interval of stress values [470, 1567MPa] .The maximum stress is 1567 MPa which exceeds the strength of the material. The stress then decreases until reaching the breaking stress σg = 470MPa. The second zone [3.9, 25mm] that corresponds to the interval of stress values [356, 470MPa]; We are observing a slight diminution of stress to stabilize at the applied stress value which is above the elastic limit. Thus, this zone is plastically deformed. 4.2 Evolution of the stress intensity factor The curves in Figure 6 represent the evolution of analytical and numerical stress intensity factor in the ligament of the specimen, for the three levels of applied stress: Δσ = 148 MPa, 284MPa, 356MPa. Volume 3, Issue 4, April 2015 Page 46 IPASJ International Journal of Mechanical Engineering (IIJME) A Publisher for Research Motivation........ Volume 3, Issue 4, April 2015 Web Site: http://www.ipasj.org/IIJME/IIJME.htm Email: [email protected] ISSN 2321-6441 Figure 6 : Evolution of stress intensity factor in the ligament of the specimen, for the three levels of applied stress: Δσ = 148 MPa, 284MPa, 356MPa. The analysis of the curves in Figure 6 shows that there is a significant increase in the stress intensity factor in function of the crack propagation and the applied stress. The variation of the numerical stress intensity factor is comparable to that of the analytical. The values of the critical notch length are reported in table 2. Table2. The values of the critical notch length depending on the stress levels the critical notch length Stress (MPa) ac (mm) 148 12 284 5 356 2 5. CONCLUSION In the sectors of unsafe structures such as pressure vessels and in the presence of defect, it is essential to detect precisely the degree of defect harmfulness. Numerical finite element modeling method is an extremely efficient tool to address this issue. A numerical model using Cast3m 2013 is performed on a tensile specimen (SENT) to study the evolution of the stress and stress intensity factor along the ligament of the specimen for three levels of applied stress (Δσ = 356MPa, 284MPa, 148MPa). Found results show that the stresses follow a parabolic tendency until reaching a maximum value at the root of the defect, and then stress stabilization is observed at the nominal stress value. The values of the critical notch length decrease with the increase of the applied stress. The finite element model adopted for this work is commonly used and can be extended to real applications. References [1] Martin Kempeneers, éléments finis statiquement admissibles et estimation d’erreur par analyse duale, Mémoire de thèse doctorat Université de Liège Faculté des Sciences Appliquées, 31 mars 2006, [2] Damien Fournier, Analyse et Développement de Méthodes de Ra nement hp en Espace pour l’Equation de Transport des Neutrons, Mémoire de thèse doctorat Université de Provence Marseille, 10 octobre 2011, [3] A. Hachim, Etude numérique et validation expérimentale des mécanismes d’endommagement et de fissuration de l’acier S355 : application aux équipements sous pression, Mémoire de thèse doctorat, Faculté des Sciences-Aïn Chock,Casablanca, 27/04/2013 [4] A. Hachim, Numerical Evaluation and Experimental Validation of Stress Concentration and Crack Propagation a Bi-notched Specimen of Steel S355, International Journal of Research and Reviews in Mechatronic Design and Simulation (IJRRMDS) Vol. 2, No. 3, September 2012, ISSN: 2046-6234 Volume 3, Issue 4, April 2015 Page 47 IPASJ International Journal of Mechanical Engineering (IIJME) A Publisher for Research Motivation........ Volume 3, Issue 4, April 2015 Web Site: http://www.ipasj.org/IIJME/IIJME.htm Email: [email protected] ISSN 2321-6441 [5] Hiroshi Yoshihara, Mode I Critical Stress Intensity Factor of MediumDensity Fiberboard Obtained by SingleEdgeNotched Bending Test ,Original scientifi c paper • Izvorni znanstveni rad. Accepted – prihvaceno: 6. 2. 2014. [6] A. El Hakimi, Numerical and experimental study of the harmfulness of defects in cylindrical and spherical shells under pressure, Ph.D. Thesis University of Technology of Compiègne, June 30, 2006. [7] A. Saffih, S. Hariri, Numerical study of elliptical cracks in cylindrers with a thickness transition, International Journal of Pressure Vessels and Piping, Vol 83, No. 1, pp. 35-41, 2006. [8] CODAP: Code de Construction des Appareils à Pression non soumis à la flame 2005 [9] Barsoum,Furthur application of quadratic isoparametric elements to linear fracture mechanics of plate bending and general shells. Int.J.Num.Meth,Engng , 11,167-169. (1976). Volume 3, Issue 4, April 2015 Page 48
© Copyright 2024