• Proceedings of the Korean Society of Computational and Applied Mathematics Conference
• /
• 2003.09a
• /
• pp.1-1
• /
• 2003

In this talk, a reduced-order modeling methodology based on centroidal Voronoi tessellations (CVT's)is introduced. CVT's are special Voronoi tessellations for which the generators of the Voronoi diagram are also the centers of mass (means) of the corresponding Voronoi cells. The discrete data sets, CVT's are closely related to the h-means clustering techniques. Even with the use of good mesh generators, discretization schemes, and solution algorithms, the computational simulation of complex, turbulent, or chaotic systems still remains a formidable endeavor. For example, typical finite element codes may require many thousands of degrees of freedom for the accurate simulation of fluid flows. The situation is even worse for optimization problems for which multiple solutions of the complex state system are usually required or in feedback control problems for which real-time solutions of the complex state system are needed. There hava been many studies devoted to the development, testing, and use of reduced-order models for complex systems such as unsteady fluid flows. The types of reduced-ordered models that we study are those attempt to determine accurate approximate solutions of a complex system using very few degrees of freedom. To do so, such models have to use basis functions that are in some way intimately connected to the problem being approximated. Once a very low-dimensional reduced basis has been determined, one can employ it to solve the complex system by applying, e.g., a Galerkin method. In general, reduced bases are globally supported so that the discrete systems are dense; however, if the reduced basis is of very low dimension, one does not care about the lack of sparsity in the discrete system. A discussion of reduced-ordering modeling for complex systems such as fluid flows is given to provide a context for the application of reduced-order bases. Then, detailed descriptions of CVT-based reduced-order bases and how they can be constructed of complex systems are given. Subsequently, some concrete incompressible flow examples are used to illustrate the construction and use of CVT-based reduced-order bases. The CVT-based reduced-order modeling methodology is shown to be effective for these examples and is also shown to be inexpensive to apply compared to other reduced-order methods.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.24 no.5
• /
• pp.553-560
• /
• 2015

Fracture mechanics analysis for cracked pipes is essential for applying the leak-before-break (LBB) concept to nuclear piping design. For LBB assessment, crack instability and leak rate should be predicted accurately for through-wall cracked pipes. In a nuclear piping system, elbows are connected with straight pipes by circumferential welding; this weld region is often considered a critical location. Hence, accurate crack assessment is necessary for cracks in the interface between elbows and straight pipes. In this study, the stress intensity factor (SIF) and elastic crack opening displacement (COD) were estimated through detailed 3D elastic finite element (FE) analyses. Based on the results, closed-form solutions of shape factors for calculating the SIFs and elastic CODs were proposed for circumferential through-wall cracks in the abovementioned interfaces under internal pressure. In addition, the effect of the elbow on shape factors was investigated by comparing the results with the existing solutions for a straight pipe.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.29 no.10 s.241
• /
• pp.1361-1368
• /
• 2005

The structural and leakage integrity of steam generator tubes should be sustained all postulated loads with appropriate margin even if a crack is present. During the past three decades, for effective integrity evaluation, several limit load solutions have been used world-widely. However, to predict accurately load carrying capacities of specific components under different conditions, the solutions have to be modified by using lots of experimental data. The purpose of this paper is to propose a new burst pressure estimation scheme based on fracture mechanics analyses for steam generator tube with an axial or circumferential through-wall crack. A series of three dimensional elastic-plastic finite element analyses were carried out and, then, closed-form estimation equations with respect to both J-integral and crack opening displacement were derived through reference stress method. The developed engineering equations were utilized for structural integrity evaluation and the resulting data were compared to the corresponding ones fiom experiments as well as limit load solutions. Thereafter, since the effectiveness was proven by promising results, it is believed that the proposed estimation scheme can be used as an efficient tool for integrity evaluation of cracked steam generator tubes.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.37 no.1
• /
• pp.41-47
• /
• 2024

This paper utilizes computational asymptotic analysis to compute the boundary layer solution for composite beams and validates the findings through a comparison with ANSYS results. The boundary layer solution, presented as a sum of the interior solution and pure boundary layer effects, necessitates a mathematically rigorous formalization for both interior and boundary layer aspects. Computational asymptotic analysis emerges as a robust technique for addressing such problems. However, the challenge lies in connecting the boundary layer and interior solutions. In this study, we systematically separate the principles of virtual work and the principles of Saint-Venant to tackle internal and boundary layer issues. The boundary layer solution is articulated by calculating the Papkovich-Fadle eigenfunctions, representing them as linear combinations of real and imaginary vectors. To address warping functions in the interior solutions, we employed a least squares method. The computed solutions exhibit excellent agreement with 2D finite element analysis results, both quantitatively and qualitatively. This validates the effectiveness and accuracy of the proposed approach in capturing the behavior of composite beams.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.20 no.2
• /
• pp.138-149
• /
• 2000

Various modeling techniques for ultrasonic wave propagation and scattering problems in finite solid media are presented. Elastodynamic boundary value problems in inhomogeneous multi-layered plate-like structures are set up for modal analysis of guided wave propagation and numerically solved to obtain dispersion curves which show propagation characteristics of guided waves. As a powerful modeling tool to overcome such numerical difficulties in wave scattering problems as the geometrical complexity and mode conversion, the Boundary Element Method(BEM) is introduced and is combined with the normal mode expansion technique to develop the hybrid BEM, an efficient technique for modeling multi mode conversion of guided wave scattering problems. Time dependent wave forms are obtained through the inverse Fourier transformation of the numerical solutions in the frequency domain. 3D BEM program development is underway to model more practical ultrasonic wave signals. Some encouraging numerical results have recently been obtained in comparison with the analytical solutions for wave propagation in a bar subjected to time harmonic longitudinal excitation. It is expected that the presented modeling techniques for elastic wave propagation and scattering can be applied to establish quantitative nondestructive evaluation techniques in various ways.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.35 no.10
• /
• pp.1329-1335
• /
• 2011

On the basis of detailed 3D finite-element (FE) limit analyses, the plastic limit load solutions for pipes with slanted circumferential through-wall cracks (TWCs) subjected to axial tension, global bending, and internal pressure are reported. The FE model and analysis procedure employed in the present numerical study were validated by comparing the present FE results with existing solutions for plastic limit loads of pipes with idealized TWCs. For the quantification of the effect of slanted crack on plastic limit load, slant correction factors for calculating the plastic limit loads of pipes with slanted TWCs from pipes with idealized TWCs are newly proposed from extensive 3D FE calculations. These slant-correction factors are presented in tabulated form for practical ranges of geometry and for each set of loading conditions.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.2
• /
• pp.291-298
• /
• 1994

This paper presents a new stiffness matrix for the analysis of arbitrary thin-walled open beams in warp-restrained torsion. The element accounts for both flexural and warping torsional effects. To eliminate the ad hoc introduction of St. Venant stiffness in this $C^0$ element, the virtual work equation based on an orthogonal Cartesian coordinate system is used. The effectiveness of the derived block stiffness is addressed. The transformation matrix between two different reference systems is also shown. Numerical examples using the proposed matrix are compared with the classical solutions or other previous results in the literature.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.3 no.4
• /
• pp.71-82
• /
• 1999

An axisymmetric shell element which includes the effects of the meridional and circumferential cable prestresses is developed. The fluid-structure interaction is expressed as added mass effect which is in proportion to the acceleration of the structure in interface surface. The added mass is obtained by using finite element method under the assumption that the fluid is invicid, incompressible and irrotational. It is coded for personal computer by the maximum use of axisymmetic properties and the dynamic analysis are performed under seismic exitations. A ring element makes the characteristics of the axisymmetric shell to be fully utilized. The elgenvalue solutons under the initial prestresses and the internal fluid are well agreed with the exact solutions and references by using under 20 elements. The eigenvalues are decreased along the increasing the height of internal fluid and these effects are dominant under the lower wave numbers. The results of the seismic analysis show that the radial deflection under the meridional prestress is a little larger than that under the circumferential prestress.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.2
• /
• pp.87-95
• /
• 2018

This paper introduces a new explicit spectral element method for the simulation of SH-waves in semi-infinite domains. To simulate the wave motion in unbounded domains, it is necessary to reduce the infinite extent to a finite computational domain of interest. To prevent the wave reflection from the trunctated boundaries, perfectly matched layer(PML) wave-absorbing boundary is introduced. The forward problem for simulating SH-waves in PML-truncated domains can be formulated as second-order PDEs. The second-order semi-discrete form of the governing PDEs is constructed by using a mixed spectral elements with Legendre-gauss-Lobatto quadrature method, which results in a diagonalized mass matrix. Then the second-order semi-discrete form is transformed to a first-order, whose solutions are calculated by the fourth-order Runge-Kutta method. Numerical examples showed that solutions of SH-wave in the two-dimensional analysis domain resulted in stable and accurate, and reflections from truncated boundaries could be reduced by using PML boundaries. Elastic wave propagation analysis using explicit time integration method may be apt for solving larger domain problems such as three-dimensional elastic wave problem more efficiently.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.24 no.1
• /
• pp.55-60
• /
• 2011

In the present work, isogeometric analysis in linear elasticity problem is conducted using the basis functions from NURBS. The objectives of isogeometric analysis introduced is to integrate both geometric modeling(CAD) and computational analysis(CAE), and this can be accomplished from direct usage of geometric modeling by NURBS as the computational mesh. The merit of the isogeometry analysis is that NURBS surface are able to represent exact geometry from the control points and knot vectors, and also subsequent refinement is relatively simple relatively. In order to verify the computer codes developed in this study, it has been applied to two structural models of which geometry are simple ; 1) circular cylinder subjected to the constant internal pressure loading, 2) square plate with circular hole at center subjected to uniform tension. The exact solutions of these two models are available. Convergence of the approximate solutions by the present code for the isogeometry analysis are investigated by mesh refinement with inserting knots (h-refinement) and by mesh refinement with order elevation of the basis functions (p-refinement).