Crack propagation rate steel

Sign In or Register for Account. Sign in via your Institution. Learn about subscription and purchase options. Product added to cart. View Metrics. Accepted Manuscript Alert. New Issue Alert.

Cited By Web Of Science FEA model based cautious automatic optimal shape control for fuselage assembly J. Five-axis trajectory generation considering synchronisation and nonlinear interpolation errors J. Offshore Mech. Eng February, Basic Eng December, Ind November, Modeling of Fatigue Crack Propagation J. Technol January, Podlasek, Observations on shear fracture propagation behavior, In: Symposium on crack propagation in pipelines, Institute of Gas Engineers, London, Murtagian and G.

Ernst, D. Johnson and H. Rudland, G. Wilkowski and Z. Osegi, An Introduction to mathematical modeling using simulink. United State: Lulu publisher, Prince O. PDF Version View. Agbainor1 and S. Abstract For the stability assessment in piping components, it is important to determine the point of initiation of a crack, and to monitor and diagnose the subsequent crack propagation behavior as well as the failure patterns, to ensure pipeline structural integrity.

Introduction Oil and Gas transmission pipelines usually have a good safety record. This is due to a combination of good design, materials and operating practices. The gamma model was used to determine the next depth of crack by introducing a non-dimensional number and crack fitting or varying the exponent values.

Materials and Methods The proposed use of the material dynamic toughness on the governing equation of the crack motion accounts for the combination of effects of the finite rate of energy transfer by the elastic waves and the static material toughness. Crack area Crack area Figure 2. Numerical loop scheme for solving crack growth problem The model here developed is based in the momentum conservation principle shown in Eq. J app da 2. As a function Replacing Eqs.

Therefore for a steady state crack propagation, when v vLIM Eq. Crack driving force Japp dm. It is assumed that this effective mass is located at the crack tip [9] body inertial forces, specimen geometry changes and stress wave reflections must be considered [5].

Conservation of mass and energy was utilized Where is the system hydraulic discharge coefficient, assumed to be of unitary value for all the experiments, A is the crack opening area, and vfl the fluid flow velocity through the fracture area, which can be computed using Bernoullis theorem, p being the hydrostatic pressure and the fluid density inside the specimen. For a cylindrical specimen of length L Eq. Typical parameterized decompression curve Dynamic material toughness JRD growth is on dependent on k we set the other parameter of k to unity, i.

Crack Area at Specified Conditions tn An 0 0 Figure 4. Plot of crack extension vs. Time From the graph above, there was a gradual increase in the crack growth within the first [secs], and thereafter the crack reaches it terminal where no further crack growth can occur i. Figure 5. A plot showing the comparison be tween Murtagians and Developed model From figure 5 above, it can be seen that unlike the Murtagians model where the constant crack extension value of about [mm] was attained rapidly within the first [secs], the developed model displayed a uniform crack extension with respect to time crack speed before the constant crack extension was attained at [sec].

Figure 6. Figure 8. Table 4 below shows the error that resulted after those numbers of simulations were carried out. Table 4. Optimize values of Monte Carlo wit h error function but increases steadily for critical material toughness 1 10 Figure Monte Carlo estimate of crack ar ea evolution 10 A plot of Y input against time Figure 11 shows the plot of the generated curve for crack growth for Levenberg-Marquardt algorithm with crack growth A as the initial guesses from developed model.

Table 5. Optimum values that satisfies the nonlinear criteria for model parameters X1 A plot of Y Estimate against Tim e The plot above shows the estimate value against time.

At the stage the curves smoothened out. The model developed for crack growth evolution in Simulink Conclusion A software program is developed for the estimation of cracks based on crack tip opening area, CTOA and limiting crack speed, LCS models for steel cracked pipeline as a check of the experimental results of the samples considered.

The algorithm is needed to help site engineers solve relative problems of this nature. References G. Robert, A. Wells, The velocity of brittle fracture, Engineering , This probabilistic model is also extended by considering a damage parameter able to account for mean stress effects. The UniGrow model, originally developed by Noroozi et al. The fatigue crack growth is decomposed into successive crack increments, each one of equal size—the material elementary finite particle size.

The fatigue crack growth rate can then be computed using Eq. Noroozi et al. Peeker and Niemi [ 10 ] propose, alternatively, the use of the Morrow's equation [ 6 ] to compute the failure of the material's representative element:. In particular, Eq. According to Noroozi et al. Application of analytical elastoplastic formulae such as the Neuber's [ 15 ] or Glinka's approaches [ 25 ] to transform the elastic stress field into actual elastoplastic stresses and strains ahead of the crack tip.

Procedures by Moftakhar et al. One way to model crack opening consists in assuming that the compressive residual stress field acting ahead of the crack tip is applied in a symmetrical way, behind the crack tip, directly on crack faces.

This compressive stress distribution acting on crack faces is equivalent to a residual stress intensity factor being used to correct the applied stress intensity factor range leading to a total effective stress intensity factor range that excludes the effects of the compressive stresses. It is proposed by Noroozi et al. Using the total values of the stress intensity factors, the first and second steps steps i and ii are repeated to determine the corrected values for the maximum actual stress and actual strain ranges at the material's representative elements.

Then, Eq. This type of fatigue crack propagation model based on the combination of two parameters crack driving force has recently being proposed by several authors [ 27 , 28 ].

In this chapter, the methodology proposed by Noroozi et al. In the following an elastoplastic constitutive model based on the von Mises yield criterion J2 plasticity and associative flow rule is selected to model compact tension specimen's elastoplastic behaviour and particularly to compute the residual stresses at the fatigue crack region. The von Mises yield criterion points out a yield function, f , defined as follows:. The plastic flow occurs according to the associative flow rule, which means that the plastic strain increment can be computed from the derivatives of the yield function with respect to the stress tensor:.

The individual increments in the plastic strain are summed up using the weighting factors determined in step i to result the total or apparent increment in plastic strain. A perfectly plastic material based on von Mises criterion is assumed.

In the plastic computation algorithm, if the equivalent von Mises stress predicted using an elastic trial exceeds the material yield stress given the yield function, then plastic strains will occur. Then a plastic correction will be required, which means that plastic strains reduce the stress state so that it satisfies the yield criterion return to mapping procedures.

The followed scheme for the integration of the elastoplastic constitutive equations consisted of an elastic trial and a return mapping procedure, as proposed by Simo and Taylor [ 31 ]. ANSYS commercial code was used in this study with the cyclic plasticity model being already part of the software capabilities [ 30 ].

The application of the UniGrow model requires a fatigue damage relation to compute the number of cycles to fail the elementary material blocks. Details of each derivation are found in Castillo et al. This model provides a complete analytical description of the statistical properties of the physical problem being dealt with, including the quantile curves as curves representing the same probability of failure. Any combination of maximum stress and strain amplitude supplying the same value of the SWT parameter should lead to the same fatigue life.

Their physical meanings see Figure 4 are N 0 : threshold value of lifetime and. Note that Eq. The Weibull parameters may be estimated by the maximum likelihood method [ 32 , 33 ].

A procedure proposed by Correia et al. The problem of the CT specimen geometry was addressed by the code developed [ 34 — 39 ]. The required data are the material properties, loading parameters and geometric dimensions of the CT specimen, including the initial and final crack sizes to be simulated.

Figure 3 gives a general overview of the procedure. Procedure to generate probabilistic fatigue crack propagation fields. The residual stresses play a central role in the UniGrow model; therefore, the accuracy of the residual stress evaluation method is vital for attaining consistent fatigue crack growth rate predictions. The PNL1 steel is a pressure vessel steel, the fatigue behaviour of which has been already analysed [ 21 , 35 , 39 , 40 ]. This section presents the cyclic elastoplastic fatigue data and the fatigue crack growth data obtained for the PNL1 steel [ 21 ].

Fatigue crack propagation data obtained for the PNL1 steel. All tests are performed in air, at room temperature, under a sinusoidal waveform at a maximum frequency of 20 Hz. The crack propagation rates are only slightly influenced by the stress ratio. Higher stress ratios provide higher crack growth rates.

The Morrow's equation does not show the vertical asymptote as verified with the Weibull field, being physically more consistent when performing extrapolations for very low number of cycles.