• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.4
• /
• pp.1051-1061
• /
• 1995

In order perform leak-before-break design of nuclear piping systems and integrity evaluation of reactor vessels, full stress-strain (.sigma. - .epsilon.) curves and fracture resistance (J-R) curves are required. However it is time-consuming and expensive to obtain J-R curves experimentally. The objective of this paper is to develop two methods for J-R curve prediction. In the first method, elastic-plastic finite element analyses for a series of crack length / specimen width ratio were performed. Accordingly the load versus load line displacement (P .delta.) curve corresponding to the fracture strain is obtained and the J-R curve based on the generalized locus method is obtained. In the second method, the correlation between .sigma.-.epsilon. curves and J-R curves was statistically analyzed and an empirical equation to predict the J-R curve from the .sigma.-.epsilon. test result is proposed. A good correlation between the predicted results based on the proposed methods and the experimental ones is obtained.

• Journal of Advanced Research in Ocean Engineering
• /
• v.3 no.1
• /
• pp.15-24
• /
• 2017

This paper explores a series of numerical simulations of dynamic responses of multi-piles (dolphin) type substructures for 2.5MW class offshore wind turbine. Firstly computational fluid dynamics (CFD) simulation was performed to evaluate wave loads on the dolphin type substructures with the design wave condition for the west-south region of Korea. Numerical wave tank (NWT) based on CFD was adopted to generate numerically a progressive regular wave using a virtual piston type wave maker. It was found that the water-piercing area of piles of the substructure is a key parameter determining the wave load exerted in horizontal direction. In the next the dynamic structural responses of substructure members under the wave load were calculated using finite element analysis (FEA). In the FEA approach, the dynamic structural responses were able to be calculated including a deformable body effect of substructure members when wave load on each member was determined by Morison's formula. The paper numerically identifies dynamic response characteristics of dolphin type substructures for 2.5MW class offshore wind turbine.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.23 no.6 s.165
• /
• pp.988-1000
• /
• 1999

An experimental study is presented for the effect of number of terms of a pewee series type stress function on stress analysis around a hole in tensile loaded plate. The hybrid method coupling photoelastsic data inputs and complex variable formulations involving conformal mappings and analytical continuity is used to calculate tangential stress on the boundary of the hole in uniaxially loaded, finite width tensile plate. In order to measure isochromatic data accurately, actual photoelastic fringe patterns are two times multiplied and sharpened by digital image processing. For qualitative comparison, actual fringes are compared with calculated ones. For quantitative comparison, percentage errors and standard deviations with respect to percentage errors are caculated for all measured points by changing the number of terms of stress function. The experimental results indicate that stress concentration factors analyzed by the hybrid method are accurate within three percent compared with ones obtained by theoretical and finite element analysis.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.1
• /
• pp.132-144
• /
• 1987

This paper presents the theroetical analysis of the crack tip stress intensity factor for a center crack in a finite width plate with variable thickness. The analyses were based on Laurent's expansions of complex stress potentials where the expansion coefficients are determined from the boundary conditions. The perturbation method was employed in numerical calculations. The correction factor F(.lambda.)is given in the form of power series of .lambda. [a numerical formula] where .lambda.=a/w$^{1}$; Dimensionless crack length, .betha.=t$_{2}$/t; Thickness ratio .omega.=w$_{2}$/w$_{1}$; width ratio The correction factor values vary with the width ratio .omega. and the maximum variation occurs around .betha.=1. For the case of .betha.=1 or .betha.=0 (uniform thickness plate0, the correction factor values agree well with Feddersen's formula. In all cases, as .lambda. approaches to 1 (thickness interface), the correction factor values are decreased rapidly for .betha.>1, and increased rapidly for .betha.<1.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2004.10a
• /
• pp.855-859
• /
• 2004

We develop a new four-node flat shell element which is accurate, efficient, and suitable to be used on general purpose. The new element has a hybrid Trefftz element with drilling degrees of freedom as a membrane part. We define the two independent displacement field: the internal displacement field that satisfies governing equations in the domain a priori and the boundary displacement field that is usually used as a conventional finite element method. The hybrid Trefftz variational formulation connects these two displacement fields on the boundary of the domain. To add drilling degrees of freedom, we introduce the Allman's quadratic displacement field to the boundary displacement field. As a result, our flat shell element has 6 degrees of freedom per a node. We also use the well-known DKMQ plate bending element for the plate part of the proposed element. The DKMQ element satisfies Mindlin-Reissner‘s plate theory along the edge of the element and gives proper behavior regardless of the thickness. A series of numerical experiments shows that the performance of the new element such as accuracy, rate of convergence, robustness to mesh quality, and so on.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.1
• /
• pp.59-65
• /
• 2004

The stress intensity factors have been widely used in numerical studies of crack growth direction. However in many cases, omissive terms of the series expansion are quantitatively significant, so we consider the computation of such terms. For this purpose, we used the finite element method with isoparametric quadratic quarter-point elements. For examples, infinite square plate with a slant crack subjected to a uniaxial load is analyzed. The numerical analysis were performed for the wide range of crack tip element lengths and inclined angles. The numerical results obtained are compared with the theoretical solutions. Also they were accurate and efficient.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.1
• /
• pp.109-118
• /
• 2003

The recent trend towards miniaturization causes an increased demand for parts with very small dimensions. Milli-structure components are classified as a component group whose size is between macro- and micro-scale. The manufacturing process of these components of thin sheet metal forming has a microscopic properties in addition to a typical phenomenon of bulk deformation because of the forming size. Also, the material properties and the deformation behavior change with miniaturization, which means that, a coarse grained materials show a higher resistance against deformation, when the grain size is in the range of the sheet thickness. In this study, a new numerical approach is proposed to simulate intergranular milli-structure in forming by the finite element method. The grain element and grain boundary element are introduced to simulate the milli-structure in the bending. The grain element is used to analyze the deformation of individual grain while the grain boundary element is for the investigation on the movement of the grain boundary. Also, the result of the finite element analysis is confirmed by a series of milli-sized forming experiments.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.27 no.4
• /
• pp.197-201
• /
• 2015

Experimental and numerical analysis has been carried out in this paper to understand correlation between volume and yield of cellular bulkhead. It was firstly confirmed from these results that the conditions and parameters considered in the finite element analysis were reasonable and realistic due to the fact that the yield loads determined by the two different methods were equivalent in actual. Based on this results, a series of intensive numerical analysis has been further performed and revealed that the yield load varied in direct proportion to the change in the volume of the cellular bulkhead.

• Journal of the Korean Institute of Gas
• /
• v.12 no.4
• /
• pp.58-62
• /
• 2008

The purpose of this paper is to establish a standard finite element analysis model of a piston by carrying out three dimensional modeling of a series six-cylindered CNG engine's piston to forecast temperature distribution at stationary state and the following thermal stress and variation, and cross checking it with existing analysis. Also, in order to evaluate the affects of the cooling system to the piston's heat load, the paper analyzed piston's temperature and thermal stress distribution according to the cooling water temperature changes and the following variations. As a result, the maximum temperature was found at the center of the crown in the piston and the maximum thermal stress occurred from the lower part of the piston.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.1
• /
• pp.69-81
• /
• 2015

The present paper deals with the determination of temperature, displacement and thermal stresses in a semi-infinite hollow circular disk due to internal heat generation within it. Initially the disk is kept at arbitrary temperature F(r, z). For times t > 0 heat is generated within the circular disk at a rate of g(r, z, t) $Btu/hr.ft^3$. The heat flux is applied on the inner circular boundary (r = a) and the outer circular boundary (r = b). Also, the lower surface (z = 0) is kept at temperature $Q_3(r,t)$ and the upper surface ($Z={\infty}$) is kept at zero temperature. Hollow circular disk extends in the z-direction from z = 0 to infinity. The governing heat conduction equation has been solved by using finite Hankel transform and the generalized finite Fourier transform. As a special case mathematical model is constructed for different metallic disk have been considered. The results are obtained in series form in terms of Bessel's functions. These have been computed numerically and illustrated graphically.