An effective crack tip region finite element sub-model for fracture mechanics analysis

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.


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.


Author(s):  
V Sura ◽  
S Mahadevan

Shattered rim cracking, propagation of a subsurface crack parallel to the tread surface, is one of the dominant railroad wheel failure types observed in North America. This crack initiation and propagation life depends on several factors, such as wheel rim thickness, wheel load, residual stresses in the rim, and the size and location of material defects in the rim. This article investigates the effect of the above-mentioned parameters on shattered rim cracking, using finite element analysis and fracture mechanics. This cracking is modelled using a three-dimensional, multiresolution, elastic–plastic finite element model of a railroad wheel. Material defects are modelled as mathematically sharp cracks. Rolling contact loading is simulated by applying the wheel load on the tread surface over a Hertzian contact area. The equivalent stress intensity factor ranges at the subsurface crack tips are estimated using uni-modal stress intensity factors obtained from the finite element analysis and a mixed-mode crack growth model. The residual stress and wheel wear effects are also included in modelling shattered rim cracking. The analysis results show that the sensitive depth below the tread surface for shattered rim cracking ranges from 19.05 to 22.23 mm, which is in good agreement with field observations. The relationship of the equivalent stress intensity factor (Δ K eq) at the crack tip to the load magnitude is observed to be approximately linear. The analysis results show that the equivalent stress intensity factor (Δ K eq) at the crack tip depends significantly on the residual stress state in the wheel. Consideration of as-manufactured residual stresses decreases the Δ K eq at the crack tip by about 40 per cent compared to that of no residual stress state, whereas consideration of service-induced residual stresses increases the Δ K eq at the crack tip by about 50 per cent compared to that of as-manufactured residual stress state. In summary, the methodology developed in this article can help to predict whether a shattered rim crack will propagate for a given set of parameters, such as load magnitude, rim thickness, crack size, crack location, and residual stress state.


Author(s):  
Viktor Gaidaichuk ◽  
Kostiantyn Kotenko

The problem of dynamic deformation of a three-layer cylindrical shell under non-stationary loads in the case of rigid clamping of the shell ends is considered. The article presents the results of assessing the stress-strain state of a three-layer cylindrical shell, taking into account its structural feature, the ratio of the sheathing thickness and the physical and mechanical characteristics of a one-piece polymer filler. Calculations were performed by software complex Nastran. The values of displacements and stresses were calculated by the algorithm of direct transient dynamic process. The step duration of the time interval was 0.0000025 s, and the total number of steps was 200. The choice of the type of three-dimensional finite element was due to obtaining more detailed and accurate calculation results. The finite element model included 19000 three-dimensional finite elements and numbered 20800 nodes. The influence of geometrical parameters of shell layers with different physical and mechanical properties of one-piece filler on the stress-strain state under axisymmetric internal impulse load is investigated. Numerical results on the dynamics of the three-layer structure, obtained by the finite element method, allow to characterize the stress-strain state of the three-layer elastic structure of the cylindrical type at any time in the studied time interval. Optimization of the shell design is recommended. Changing the ratio of the thickness of the internal and external shells of the shell significantly affects the stress-strain state of the shell and its performance. Increasing the thickness of the internal layer of the shell significantly contributes to the increase of the latter. Comparison of the given results with materials of other similar researches and positions, testify to objectivity of the made approach.


Author(s):  
Curtis Sifford ◽  
Ali Shirani

Abstract This paper presents the application of the rules from ASME Section VIII, Division 3 of the ASME Boiler and Pressure Vessel Code for a fracture mechanics evaluation to determine the damage tolerance and fatigue life of a flowline clamp connector. The guidelines from API 579-1 / ASME FFS-1 Fitness-For-Service for the stress analysis of a crack-like flaw have been considered for this assessment. The crack tip is modeled using a refined mesh around the crack tip that is referred to as a focused mesh approach in API 579-1 / ASME FFS-1. The driving force method is used as an alternative to the failure assessment diagram method to account for the influence of crack tip plasticity. The J integral is determined using elastic-plastic finite element analysis and converted to an equivalent stress intensity factor to be compared to the fracture toughness of the material. The fatigue life is calculated using the Paris Law equation and the stress intensity factor calculated from the finite element analysis. The allowable number of design cycles is determined using the safety factors required from Division 3 of the ASME Pressure Vessel Code.


Author(s):  
Curtis Sifford ◽  
Ali Shirani

This paper presents the application of the rules from ASME Section VIII, Division 3 of the ASME Boiler and Pressure Vessel Code for a fracture mechanics evaluation to determine the damage tolerance and fatigue life of a flowline clamp connector. The guidelines from API 579-1 / ASME FFS-1 Fitness-For-Service for the stress analysis of a crack-like flaw have been considered for this assessment. The crack tip is modeled using a refined mesh around the crack tip that is referred to as a focused mesh approach in API 579-1 / ASME FFS-1. The driving force method is used as an alternative to the failure assessment diagram method to account for the influence of crack tip plasticity. The J integral is determined using elastic-plastic finite element analysis and converted to an equivalent stress intensity factor to be compared to the fracture toughness of the material. The fatigue life is calculated using the Paris Law equation and the stress intensity factor calculated from the finite element analysis. The allowable number of design cycles is determined using the safety factors required from Division 3 of the ASME Pressure Vessel Code.


Author(s):  
V.V. LEONTYEV

The method for analyzing of stress-strain state characteristics of unloaded riveted joints performed with OST 1 11781-74 rivets has been developed using Coupled Euler-Lagrange finite element approach implemented in the CAD / CAE system Abaqus. A comparative analysis of the stress-strain state characteristics of the examined riveted joint’s finite element models using the Lagrangian and the Coupled Lagrangian-Eulerian finite element approaches has been conducted. A three- dimensional finite element model based on the CLE method has been proposed for further study of fatigue strength and durability of the loaded riveted joints.


Sign in / Sign up

Export Citation Format

Share Document