Finite elements analysis of the side grooved I-beam specimen for mode II fatigue crack growth rates determination

Author(s):  
T.M. Lenkovskiy ◽  
V.V. Kulyk ◽  
Z.A. Duriagina ◽  
L.V. Dzyubyk ◽  
V.V. Vira ◽  
...  

Purpose: Carefully investigate the stress-strain state of the side grooved I-beam specimen with edge crack and determine the effect of crack length and crack faces friction on stress intensity factor at transverse shear. Design/methodology/approach: The finite element method was used to estimate the stress-strain state of I-beam specimen at transverse shear. For this purpose, a fullscale, three-dimensional model of the specimen was created, which precisely reproduces its geometry and fatigue crack faces contact. For the correct reproduction of the stress singularity at the crack tip, a special sub-model was used, which has been tested earlier in solving similar problems of fracture mechanics. In order to improve the accuracy of the calculations, for crack plane and cross-section of the specimen on the crack extension modeling, an algorithm for changing the crack length without changing the total number of elements in the model was developed and applied. Young's modulus and Poisson's ratio of structural steels were specified for the model material. The static loading of the model was realized assuming small scale yielding condition. The stress intensity factor was found through the displacement of nodes in the prismatic elements adjacent to the plane and the front of the crack. Findings: Mathematical dependences, which show an increase of stress intensity factor in the I-beam specimen with an increase in the crack length, and its decrease with an increase of crack faces friction factor at transverse shear, were established. The results are compared with the partial cases known from the literature and their good convergence was shown. Research limitations/implications: By analyzing the obtained graphical dependences, it is established that for relative crack lengths less than 0.4 there is a significant influence of the initial notch on the stress-strain state of the specimen, and for the lengths greater than 0.9 an influence of constrained gripping part took place. For this reason, all subsequent calculations were carried out in the range of relative crack length from 0.4 to 0.9, which represents the applicability range of the final calculation formula. Increasing of the crack faces friction factor from 0 to 1 monotonically reduces the stress at the crack tip. For a short crack, this effect is 1.5 times greater than for a long one, which is reflected by the calculation formula. Practical implications: Using the proposed calculation formula, one can calculate the stress intensity factor in the I-beam specimen, and to determine the crack growth resistance characteristics of structural steels at transverse shear. Originality/value: A new, easy-to-use in engineering calculations formula is proposed for stress intensity factor determination in the I-beam specimen at transverse shear. The formula takes into account crack faces friction for various crack lengths.

2017 ◽  
Vol 2 (87) ◽  
pp. 56-65 ◽  
Author(s):  
T.M. Lenkovskiy ◽  
V.V. Kulyk ◽  
Z.A. Duriagina ◽  
R.A. Kovalchuk ◽  
V.G. Topilnytskyy ◽  
...  

Purpose: To create an effective in engineering strength calculation three-dimensional submodel of the near crack tip region in solids for hi-fidelity analysis of their stress-strain state by the finite element method. Design/methodology/approach: To create a volume near the crack tip, regular threedimensional 20-node prismatic isoparametric elements and 15-node special elements with edge length of 12.5 μm with shifted nodes in order to simulate the singularity of stress at the crack tip were used. Using these two types of elements, a cylindrical fragment of diameter of 100 μm was built. In its base is a 16-vertex polygon, and its axis is the crack front line. In the radial direction the size of the elements was smoothly enlarged by creating of 5 circular layers of elements, and in the axial direction 8 layers were created. For convenience of the sub-model usage, the cylindrical fragment was completed by regular elements to a cubic form with edge size 400 μm. For the sub-model approbation, the full-scale three-dimensional models of standard specimens with cracks were built. The stress intensity factor K at normal tension was calculated assuming small scale yielding conditions in a plane between 4th and 5th layers of special elements on the basis of analysis of displacement fields near the crack tip. Findings: An effective three-dimensional sub-model of the near crack tip region is proposed. The sub-model was used to obtain the dependence of the stress intensity factor on the relative crack length at normal tension for four types of standard specimens. The obtained dependences show excellent correlation with known analytical solutions. Research limitations/implications: The concept of finite element meshing at threedimensional modelling of the near crack tip region for high-fidelity stress-strain state analysis was generalized. A sub-model of the near crack tip region was created and used to determine the stress intensity factor at normal tension of four types of standard specimens. It is shown that the proposed methodology is effective for precise analysis of the stressstrain state of solids with cracks within the framework of linear fracture mechanics. Practical implications: By applying the generalized approach and the proposed threedimensional sub-model of the near crack tip region, one can determine the stress-strain state of structure elements and machine parts when analysing their workability by the finite element method. Originality/value: An effective finite-element sub-model for the stress-strain state analysis in the vicinity of the crack tip within the framework of the linear fracture mechanics is proposed.


Author(s):  
D. E. Tulin ◽  

The study of the elastic and elastoplastic stress-strain state of a plate with a semi-elliptical crack when changing its geometry has been carried out. The process of transition of a semi-elliptical crack into a through one was considered. Options of the stress-strain state were determined by means of the FEA method. An approximating dependence of the stress intensity factor on geometry of the plate is suggested. Keywords: strength, semi-elliptical crack, K-calibration, stress intensity factor, elastoplastic deformation, stress state intensity.


2017 ◽  
Vol 8 (5) ◽  
pp. 506-515 ◽  
Author(s):  
Raviraj M.S. ◽  
Sharanaprabhu C.M. ◽  
Mohankumar G.C.

Purpose The purpose of this paper is to present the determination of critical stress intensity factor (KC) both by experimental method and three-dimensional (3D) finite element simulations. Design/methodology/approach CT specimens of different compositions of Al6061-TiC composites (3wt%, 5wt% and 7wt% TiC) with variable crack length to width (a/W=0.3-0.6) ratios are machined from as-cast composite block. After fatigue pre-cracking the specimens to a required crack length, experimental load vs crack mouth opening displacement data are plotted to calculate the KC value. Elastic 3D finite element simulations have been conducted for CT specimens of various compositions and a/W ratios to compute KC. The experimental results indicate that the magnitude of KC depends on a/W ratios, and significantly decreases with increase in a/W ratios of the specimen. Findings From 3D finite element simulation, the KC results at the centre of CT specimens for various Al6061-TiC composites and a/W ratios show satisfactory agreement with experimental results compared to the surface. Originality/value The research work contained in this manuscript was conducted during 2015-2016. It is original work except where due reference is made. The authors confirm that the research in their work is original, and that all the data given in the article are real and authentic. If necessary, the paper can be recalled, and errors corrected.


Sign in / Sign up

Export Citation Format

Share Document