• Nuclear Engineering and Technology
• /
• 제50권7호
• /
• pp.1138-1147
• /
• 2018

Delayed hydride cracking (DHC) is an important failure mechanism for Zircaloy tubes in the demanding environment of nuclear reactors. The threshold stress intensity factor, $K_{IH}$, and critical hydride length, $l_C$, are important parameters to evaluate DHC. Theoretical models of them are developed for Zircaloy tubes undergoing non-homogenous temperature loading, with new stress distributions ahead of the crack tip and thermal stresses involved. A new stress distribution in the plastic zone ahead of the crack tip is proposed according to the fracture mechanics theory of second-order estimate of plastic zone size. The developed models with fewer fitting parameters are validated with the experimental results for $K_{IH}$ and $l_C$. The research results for radial cracking cases indicate that a better agreement for $K_{IH}$ can be achieved; the negative axial thermal stresses can lessen $K_{IH}$ and enlarge the critical hydride length, so its effect should be considered in the safety evaluation and constraint design for fuel rods; the critical hydride length $l_C$ changes slightly in a certain range of stress intensity factors, which interprets the phenomenon that the DHC velocity varies slowly in the steady crack growth stage. Besides, the sensitivity analysis of model parameters demonstrates that an increase in yield strength of zircaloy will result in a decrease in the critical hydride length $l_C$, and $K_{IH}$ will firstly decrease and then have a trend to increase with the yield strength of Zircaloy; higher fracture strength of hydrided zircaloy will lead to very high values of threshold stress intensity factor and critical hydride length at higher temperatures, which might be the main mechanism of crack arrest for some Zircaloy materials.

• 한국공작기계학회:학술대회논문집
• /
• 한국공작기계학회 2001년도 춘계학술대회 논문집(한국공작기계학회)
• /
• pp.72-77
• /
• 2001

The use of hexagonal ceramic inserts for copper extrusion dies offers significant technical and economic advantages over other forms of manufacture. In this paper the data on the loading of the tools is determined from a commercial FEM package as the contact stress distribution on the die-workpiece interface and as temperature distributions in the die. This data can be processed as load input data for a finite element die-stress analysis. Process simulation and stress analysis are thus combined during the design, and a data exchange program has been developed that enables optimal design of the dies taking into account the elastic deflections generated in shrink fitting the die inserts and that caused by the stresses generated in the process. The stress analysis of the dies is used to determine the stress conditions on the ceramic insert by considering contact and interference effects under both mechanical and thermal loads.

• 설비공학논문집
• /
• 제27권10호
• /
• pp.515-519
• /
• 2015

Conjugate heat transfer analysis for an ethylene furnace was carried out based on numerical simulation. Detailed distributions of velocity vectors, chemical species, and temperature inside the furnace are presented and discussed. Von Mises stress and heat flux at the tube surface were also evaluated to elucidate mechanisms regarding failure of the tube. Maximum stress was found at the upstream elbow of the tube, which did not coincide with the location of maximum heat flux. This implies that thermal stress at a steady state would not be a dominant component of the stress. Degradation of the material, as well as the system arrangement, should be considered in order to accurately predict the lifetime of the tube material in the furnace.

• Nuclear Engineering and Technology
• /
• 제56권7호
• /
• pp.2859-2874
• /
• 2024

This paper presents the finite element deformation and failure simulation of a typical Korean high-power reactor vessel under a severe accident characterized by large break loss of coolant (LBLOCA) with in-vessel retention of molten corium through external reactor vessel cooling (IVR-ERVC) conditions. Temperature distributions calculated using Modular Accident Analysis Program Version 5 (MAAP5) as thermal boundary conditions were used, and ABAQUS thermal and structural analyses were performed. After full ablation, the temperature of the inner surface in the thinnest section remained high (920 ℃), but the stress remained relatively low (less than 6 MPa). At the outer surface, the stress was as high as 250 MPa; however, the resulting plastic strain was small owing to the low temperature of 200 ℃. Variations in stress, inelastic strain, and temperature with time in the thinnest section suggest that the plastic and creep strains are saturated owing to stress relaxation, resulting in low cumulative damage. Thus, the lower head of the vessel can maintain its structural integrity under LBLOCA with IVR-ERVC conditions. The sensitivity analysis of internal pressure indicates the occurrence of failure in the thinnest section at an internal pressure >9.6 MPa via local necking followed by failure due to high stresses.

• Structural Engineering and Mechanics
• /
• 제40권3호
• /
• pp.347-371
• /
• 2011

This paper focuses on post-buckling analysis of Timoshenko beams with various boundary conditions subjected to a non-uniform thermal loading by using the total Lagrangian Timoshenko beam element approximation. Six types of support conditions for the beams are considered. The considered highly non-linear problem is solved by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. As far as the authors know, there is no study on the post-buckling analysis of Timoshenko beams under uniform and non-uniform thermal loading considering full geometric non-linearity investigated by using finite element method. The convergence studies are made and the obtained results are compared with the published results. In the study, the relationships between deflections, end rotational angles, end constraint forces, thermal buckling configuration, stress distributions through the thickness of the beams and temperature rising are illustrated in detail in post-buckling case.

• Coupled systems mechanics
• /
• 제13권2호
• /
• pp.115-132
• /
• 2024

This paper introduces an improved shear deformation theory for analyzing the buckling behavior of functionally graded plates subjected to varying temperatures. The transverse shear strain functions employed satisfy the stress-free condition on the plate surfaces without requiring shear correction factors. The material properties and thermal expansion coefficient of the porous functionally graded plate are assumed temperature-dependent and exhibit continuous variation throughout the thickness, following a modified power-law distribution based on the volume fractions of the constituents. Moreover, the study considers the influence of porosity distribution on the buckling of the functionally graded plates. Thermal loads are assumed to have uniform, linear, and nonlinear distributions through the thickness. The obtained results, considering the effect of porosity distribution, are compared with alternative solutions available in the existing literature. Additionally, this study provides comprehensive discussions on the influence of various parameters, emphasizing the importance of accounting for the porosity distribution in the buckling analysis of functionally graded plates.

• 한국항공우주학회지
• /
• 제32권8호
• /
• pp.91-100
• /
• 2004

경사기능재료 판에 대한 열탄성 변형과 응력 해석을 위해 Green 함수 방법이 채택되었다. 3 차원 비정상 온도분포에 대한 해는 적층판 이론에 의해 얻어진다. 열탄성 문제에 대한 기본 방정식은 각각 평면외 (out-plane) 변형과 평면내 (in-plane) 힘에 의해 유도되었다. 굽힘과 평면내 힘에 의한 열탄성 변형과 응력분포는 Galerkin 방법에 근거한 Green 함수를 이용하여 해석하였다. 열탄성 변형과 응력분포 해석을 위한 Galerkin Green 함수의 특성함수들은 사각판의 제차 경계조건을 만족시키는 허용함수들의 급수 형태로 근사화되었다. 단수지시된 사각 판에 대한 수치해석이 수행되었으며, 정사기능재료의 물성치가 판의 비정상 열탄성 거동에 미치는 영향이 검토되었다.

• 열처리공학회지
• /
• 제18권4호
• /
• pp.247-255
• /
• 2005

Although heat treatment is a process of great technological importance in order to obtain desired mechanical properties such as hardness, the process was required a tedious and expensive experimentation to specify the process parameters. Consequently, the availability of reliable and efficient numerical simulation program would enable easy specification of process parameters to achieve desired microstructure and mechanical properties without defects like crack and distortion. In present work, the developed numerical simulation program could predict distributions of microstructure and thermal stress in steels under different cooling conditions. The computer program is based on the finite difference method for temperature analysis and microstructural changes and the finite element method for thermal stress analysis. Multi-phase decomposition model was used for description of diffusional austenite decompositions in low alloy steels during cooling after austenitization. The model predicts the progress of ferrite, pearlite, and bainite transformations simultaneously during quenching and estimates the amount of martensite also by using Koistinen and Marburger equation. To verify the developed program, the calculated results are compared with experimental ones of casting product. Based on these results, newly designed heat treatment process is proposed and it was proved to be effective for industry.

• Structural Engineering and Mechanics
• /
• 제78권1호
• /
• pp.1-13
• /
• 2021

In this study, the porothermoelastic problem with the effect of the magnetic field and initial stress was investigated. We applied normal mode analysis to solve the resulting non-dimensional coupled equations. Numerical results for the displacements, temperature distribution, pore pressure, stresses, induced electric field and induced magnetic field distributions are presented graphically and discussed. The medium deformed because of thermal shock and due to the application of the magnetic field, there result an induced magnetic and an induced electric field in the medium. Numerical analyses are given graphically on the square (2D) and cubic (3D) domains to illustrate the effects of the porosity parameter, magnetic field and initial stress parameter on the physical variables.

• Structural Engineering and Mechanics
• /
• 제25권4호
• /
• pp.445-466
• /
• 2007

As the LNG (Liquefied Natural Gas) tank contains cryogenic liquid, realistic thermal analyses are of a primary importance for a successful design. The structural details of the LNG tank are so complicated that some strategies are necessary to reasonably predict its temperature distribution. The proposed heat transfer model can consider the beneficial effects of insulation layers and a suspended deck on temperature distribution of the outer concrete tank against cryogenic conditions simply by the boundary conditions of the outer tank model. To this aim, the equilibrium condition or heat balance in a steady state is utilized in a various way, and some aspects of heat transfer via conduction, convection and radiation are implemented as necessary. Overall thermal analysis procedures for the LNG tank are revisited to examine some unjustifiable assumptions of conventional analyses. Concrete and insulation properties under cryogenic condition and a reasonable conversion procedure of the temperature-induced nonlinear stress into the section forces are discussed. Numerical examples are presented to verify the proposed schemes in predicting the actual temperature and stress distributions of the tank as affected by the cryogenic LNG for the cases of normal operation and leakage from the inner steel tank. It is expected that the proposed schemes enable a designer to readily detect the effects of insulation layers and a suspended deck and, therefore, can be employed as a useful and consistent tool to evaluate the thermal effect in a design stage of an LNG tank as well as in a detailed analysis.