• Structural Engineering and Mechanics
• /
• v.70 no.1
• /
• pp.43-54
• /
• 2019

Due to various numerical problems, crack analysis of reinforced concrete members using the finite element method is confronting with substantial difficulties, rendering the prediction of crack patterns and crack widths a formidable task. The root cause is that the conventional analysis methods are not capable of tracking the crack sequence and accounting for the stress relief and re-distribution during cracking. To address this deficiency, the crack queuing algorithm has been proposed. Basically, at each load increment, iterations are carried out and within each iteration step, only the most critical concrete element is allowed to crack and the stress re-distribution is captured in subsequent iteration by re-formulating the cracked concrete element and re-analysing the whole concrete structure. To demonstrate the effectiveness of the crack queuing algorithm, crack analysis of concrete members tested in the literature is performed with and without the crack queuing algorithm incorporated.

• Proceedings of the Korea Concrete Institute Conference
• /
• 1999.04a
• /
• pp.233-238
• /
• 1999

Localized failure analysis of concrete structures can be carried out effectively by modeling fracture process zone of concrete during crack initiation and propagation. But, the analysis techniques are still insufficient for crack modeling because of difficulties in numerical analysis procedure which describe progressive crack. In this paper, a finite element with embedded displacement discontinuity is introduced to remove the difficulties of remeshing for crack propagation in discrete crack model during progressive failure analysis of concrete structures. The performance of this so-called embedded crack approach for concrete failure analysis is verified by several analysis examples. The analysis results show that the embedded crack approach retains mesh size objectivity and can simulate localized failure under mixed mode loading. It can be concluded that the embedded crack approach cab be an effective alternate to the smeared and discrete crack approaches.

• Journal of Mechanical Science and Technology
• /
• v.18 no.11
• /
• pp.1989-1995
• /
• 2004

Finite element analysis(FEA) is the most popular numerical method to simulate plasticity-induced fatigue crack closure and can predict fatigue crack closure behavior. Finite element analysis under plane stress state using 4-node isoparametric elements is performed to investigate the detailed closure behavior of fatigue cracks and the numerical results are compared with experimental results. The mesh of constant size elements on the crack surface can not correctly predict the opening level for fatigue crack as shown in the previous works. The crack opening behavior for the size mesh with a linear change shows almost flat stress level after a crack tip has passed by the monotonic plastic zone. The prediction of crack opening level presents a good agreement with published experimental data regardless of stress ratios, which are using the mesh of the elements that are in proportion to the reversed plastic zone size considering the opening stress intensity factors. Numerical interpolation results of finite element analysis can precisely predict the crack opening level. This method shows a good agreement with the experimental data regardless of the stress ratios and kinds of materials.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.11 s.254
• /
• pp.1455-1460
• /
• 2006

The relationship between fatigue crack growth behavior and cyclic crack tip opening displacement is studied. An elastic-plastic finite element analysis (FEA) is performed to examine the growth behavior of fatigue crack, where the contact elements are used in the mesh of the crack tip area. We investigate the relationship between the reversed plastic zone size and the changes of the cyclic crack tip opening displacement along the crack growth. We investigate the effect of the element size when predict fatigue crack opening behavior using the cyclic crack tip opening displacement obtained from FEA. The cyclic crack tip opening displacement is related to fatigue crack opening behavior.

• Journal of Mechanical Science and Technology
• /
• v.14 no.4
• /
• pp.401-407
• /
• 2000

An elastic-plastic finite element analysis is performed to investigate detailed closure behavior of fatigue cracks and the numerical results are compared with experimental results. The finite element analysis performed under plane stress using 4-node isoparametric elements can predict fatigue crack closure behavior. The mesh of constant element size along crack surface can not predict the opening level of fatigue crack. The crack opening level for the constant mesh size increases linearly from initial crack growth. The crack opening level for variable mesh size, is almost flat after crack tip has passed the monotonic plastic zone. The prediction of crack opening level using the variable mesh size proportioning the reversed plastic zone size with the opening stress intensity factors presents a good agreement with the experimental data regardless of stress ratios.

• Journal of Welding and Joining
• /
• v.21 no.6
• /
• pp.77-83
• /
• 2003

The finite element analysis was performed for the cracks existing in residual stress fields in order to investigate the effects of configuration of residual stress distribution to the fatigue crack opening behaviour. And the variation of stress distributions adjacent to the crack caused by uploading was examined. The finite element model with contact elements for the crack plane and plane stress elements for the base material and the analytical method based on the superposition principle to estimate crack opening behaviour and the stress distribution adjacent to the crack subjected to uploading were used. The results of the analysis showed that crack opening behaviors and variations of stress distribution caused by uploading were changed depending on the configuration of residual stress distribution. When the crack existed in the region of compressive residual stress and the configuration of compressive residual stress distribution were inclined, a partial crack opening just behind of a crack tip occurred during uploading. Based on the above results, it was clarified that the crack opening behaviour in the residual stress field could be predicted accurately by the finite element analysis using these analytical method and model.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.202-207
• /
• 2004

An elastic-plastic finite element analysis of fatigue crack closure is performed for plane strain conditions. The stabilization behavior of crack opening level and the effect of mesh size on the crack opening stress are investigated. In order to obtain a stabilized crack opening level for plane strain conditions, the crack must be advanced through approximately four times the initial monotonic plastic zone. The crack opening load tends to increase with the decrease of mesh size. The mesh size nearly equal to the theoretical plane strain cyclic plastic zone size may provide reasonable numerical results comparable with experimental crack opening data. The crack opening behavior is influenced by the crack growth increment and discontinuous opening behavior is observed. A procedure to predict the most appropriate mesh size for different stress ratio is suggested. Crack opening loads predicted by the FE analysis based on the procedure suggested resulted in good agreement with experimental ones within the error of 5 %. Effect of the distance behind the crack tip on the crack opening load determined by the ASTM compliance offset method based on the load-displacement relation and by the rotational offset method based on the load-differential displacement relation is investigated. Optimal gage location and method to determine the crack opening load is suggested.

• Journal of Mechanical Science and Technology
• /
• v.16 no.8
• /
• pp.1117-1126
• /
• 2002

A conducting crack in an electrostrictive ceramic under combined electric and mechanical loading is investigated. Analysis based on linear dielectric model predicts that the surfaces of the crack are not open completely but they are contact near the crack tip. The complete solution for the crack with a contact zone in a linear electrostrictive ceramic under combined electric and mechanical loading is obtained by using the complex variable formula. The asymptotic problems for a semi-infinite crack with a partial opening zone as well as for a fully open semi-infinite crack in a nonlinear electrostrictive ceramic are analyzed in order to investigate the effect of the electrical nonlinearity on the stress intensity factor under small scale nonlinear conditions. Particular attention is devoted to a finite crack in the nonlinear electrostrictive ceramic subjected to combined electric and mechanical loading. The stress intensity factor for the finite crack under small scale nonlinear conditions is obtained from the asymptotic analysis.

• Journal of Ship and Ocean Technology
• /
• v.11 no.2
• /
• pp.26-33
• /
• 2007

For the assessment of the retardation of fatigue crack propagation behavior due to overload, new FE analysis algorithms considering compressive residual stress redistribution near crack tip was proposed in this paper. The size of plastic zone near crack tip was obtained by elasto-plastic analysis and it was compared with Irwin's equation. The amount of residual stress redistribution was assessed by subsequent elasto-plastic analysis, and the difference of residual stress distributions between constant amplitude load and overload was obtained. In the analysis of fatigue crack propagation, the applied SIF range was evaluated by ASTM E647, and the effect of residual stresses on crack propagation was considered using the effective SIF concept. The test results of crack propagations were compared with the predicted data obtained by the analysis.

• Journal of the Korean Society for Precision Engineering
• /
• v.15 no.7
• /
• pp.91-97
• /
• 1998

The analysis model is the infinite body consisted of power law creep material containing a rigid inclusion with line crack shape subjected to the arbitrarily directional stress on an infinite boundary. The crack analysis is performed using the complex pseudo-stress function. The strain rate intensity factor is determined in the closed form as new fracture mechanics parmeter which represents the magnitudes of stress and strain rate near the tip in power law creep material.