• Structural Engineering and Mechanics
• /
• v.76 no.4
• /
• pp.505-512
• /
• 2020

This paper presents an error estimation technique for 2-D crack analysis by an enriched natural element (more exactly, enriched Petrov-Galerkin NEM). A bare solution was approximated by PG-NEM using Laplace interpolation functions. Meanwhile, an accurate quasi-exact solution was obtained by a combined use of enriched PG-NEM and the global patch recovery. The Laplace interpolation functions are enriched with the near-tip singular fields, and the approximate solution obtained by enriched PG-NEM was enhanced by the global patch recovery. The quantitative error amount is measured in terms of the energy norm, and the accuracy (i.e., the effective index) of the proposed method was evaluated using the errors which obtained by FEM using a very fine mesh. The error distribution was investigated by calculating the local element-wise errors, from which it has been found that the relative high errors occurs in the vicinity of crack tip. The differences between the enriched and non-enriched PG-NEMs have been investigated from the effective index, the error distribution, and the convergence rate. From the comparison, it has been justified that the enriched PG-NEM provides much more accurate error information than the non-enriched PG-NEM.

• Nuclear Engineering and Technology
• /
• v.54 no.6
• /
• pp.2084-2093
• /
• 2022

Probabilistic safety assessment is widely used to quantify the risks of nuclear power plants and their uncertainties. When the lognormal distribution describes the uncertainties of basic events, the uncertainty of the top event in a fault tree is approximated with the sum of lognormal random variables after minimal cutsets are obtained, and rare-event approximation is applied. As handling complicated analytic expressions for the sum of lognormal random variables is challenging, several approximation methods, especially Monte Carlo simulation, are widely used in practice for uncertainty analysis. In this study, a theoretical approach for analyzing the sum of lognormal random variables using an efficient numerical integration method is proposed for uncertainty analysis in probability safety assessments. The change of variables from correlated random variables with a complicated region of integration to independent random variables with a unit hypercube region of integration is applied to obtain an efficient numerical integration. The theoretical advantages of the proposed method over other approximation methods are shown through a benchmark problem. The proposed method provides an accurate and efficient approach to calculate the uncertainty of the top event in probabilistic safety assessment when the uncertainties of basic events are described with lognormal random variables.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.28 no.9
• /
• pp.1314-1319
• /
• 2004

In this study, parameter optimization of micro-static mixer with a cantilever beam was accomplished for maximizing the mixing efficiency by using successive response surface approximations. Variables were chosen as the length of cantilever beam and the angle between horizontal and the cantilever beam. Sequential approximate optimization method was used to deal with both highly nonlinear and non-smooth characteristics of flow field in a micro-static mixer. Shape optimization problem of a micro-static mixer can be divided into a series of simple subproblems. Approximation to solve the subproblems was performed by response surface approximation, which does not require the sensitivity analysis. To verify the reliability of approximated objective function and the accuracy of it, ANOVA analysis and variables selection method were implemented, respectively. It was verified that successive response surface approximation worked very well and the mixing efficiency was improved very much comparing with the initial shape of a micro-static mixer.

• Proceedings of the KSR Conference
• /
• 2004.06a
• /
• pp.589-594
• /
• 2004

In design of suspension system for KHST, it was applied the approximated optimization method using meta-models which called Response Surface Model and Neural Network Model for 29 design variables and 46 performance index. These models was coded using correlation between design variables and performance indices that is made by the 66 times iterative execution through the design of experimental table consisted orthogonal array L32 and D-Optimal design table. The results show that the optimization process is very efficient and simply applicable for complex mechanical system such as railway vehicle system. Also it was compared with the sensitivity of some design variables in order to know the characteristics of two models. This paper describes the general method for dynamic analysis and design process of railway vehicle system applied to KHST development, and proposed the efficient methods for vibration mode analysis process dealing with test data and the function based approximation method using meta-model applicable for a complex mechanical system. This method will be able to apply to the other railway vehicle system in oder to systematize and generalize the design process of railway vehicle dynamic system.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.2
• /
• pp.123-131
• /
• 2005

According to ow previous study, we confirmed That the Petrov-Galerkin natural element method(PG-NEM) completely resolves the numerical integration inaccuracy in the conventional Bubnov-Galerkin natural element method(BG-NEM). This paper is an extension of PG-NEM to two-dimensional geometrically nonlinear problem. For the analysis, a linearized total Lagrangian formulation is approximated with the PS-NEM. At every load step, the grid points ate updated and the shape functions are reproduced from the relocated nodal distribution. This process enables the PG-NEM to provide more accurate and robust approximations. The representative numerical experiments performed by the test Fortran program, and the numerical results confirmed that the PG-NEM effectively and accurately approximates The large deformation problem.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.10 no.3
• /
• pp.440-448
• /
• 1999

This paper describes a full wave analysis of the scattering from electromagnetic absorbers which can be approximated as infinite periodic structure using hybrid finite element method. By introducing fictitious boundaries, equivalent finite region is defined and proper boundary conditions of each boundary are obtained by Floquet theorem. Since higher-order Floquet modes are employed, the method presented in this paper can be readily applied to the periodic structure haying a relatively long period. To reduce difficulty in evaluating the surface integral, the normal component to the surface were represented with the tangential component to the surface. Comparisons of calculated results with analytical or published ones show the validation of the method.

• Structural Engineering and Mechanics
• /
• v.17 no.1
• /
• pp.1-14
• /
• 2004

Numerical solution to linear bending analysis of circular plates is obtained by the method of harmonic differential quadrature (HDQ). In the method of differential quadrature (DQ), partial space derivatives of a function appearing in a differential equation are approximated by means of a polynomial expressed as the weighted linear sum of the function values at a preselected grid of discrete points. The method of HDQ that was used in the paper proposes a very simple algebraic formula to determine the weighting coefficients required by differential quadrature approximation without restricting the choice of mesh grids. Applying this concept to the governing differential equation of circular plate gives a set of linear simultaneous equations. Bending moments, stresses values in radial and tangential directions and vertical deflections are found for two different types of load. In the present study, the axisymmetric bending behavior is considered. Both the clamped and the simply supported edges are considered as boundary conditions. The obtained results are compared with existing solutions available from analytical and other numerical results such as finite elements and finite differences methods. A comparison between the HDQ results and the finite difference solutions for one example plate problem is also made. The method presented gives accurate results and is computationally efficient.

• Nuclear Engineering and Technology
• /
• v.52 no.2
• /
• pp.223-229
• /
• 2020

The present work proposes a solution to the static Boltzmann transport equation approximated by the simplified P₃ (SP₃) on angular, and the analytic coarse mesh finite difference (ACMFD) for spatial variables. Multi-group SP₃-ACMFD equations in 3D rectangular geometry are solved using the GMRES solution technique. As the core time dependent analysis necessitates the solution of an eigenvalue problem for an initial condition, this work is hence devoted to development and verification of the proposed static SP₃-ACMFD solver. A 3D multi-group static diffusion solver is also developed as a byproduct of this work to assess the improvement achieved using the SP₃ technique. Static results are then compared against transport benchmarks to assess the proximity of SP₃-ACMFD solutions to their full transport peers. Results prove that the approach can be considered as an acceptable interim approximation with outputs superior to the diffusion method, close to the transport results, and with the computational costs less than the full transport approach. The work would be further generalized to time dependent solutions in Part II.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.5
• /
• pp.976-983
• /
• 1988

An efficient unbalance response analysis method for rotor-bearing systems with spin speed dependent parameters is developed by utilizing a generalized modal analysis scheme. The spin speed dependent eigenvalue problem of the original system is transformed into the spin speed independent eigenvalue problem by introducing a lambda matrix, assuming the bearing dynamic coefficients are well approximated by polynomial functions of spin speed. This method features that it requires far less computational effort in unbalance response calculations and that the influence coefficients are readily available. In addition, the critical speeds and the corresponding logarithmic decrements can be readily identified from the resulting eigenvalues.

• Structural Engineering and Mechanics
• /
• v.64 no.4
• /
• pp.473-485
• /
• 2017

An isolated building, composed of superstructure and isolation system which have very different damping properties, is typically non-classical damping system. This results in inapplicability of traditional response spectrum method for isolated buildings. A multidimensional response spectrum method based on complex mode superposition is herein introduced, which properly takes into account the non-classical damping feature in the structure and a new method is developed to estimate velocity spectra from the commonly used displacement or pseudo-acceleration spectra based on random vibration theory. The error of forced decoupling method, an approximated approach, is discussed in the viewpoint of energy transfer. From the base-isolated benchmark model, as a numerical example, application of the procedure is illustrated companying with comparison study of time-history method, forced decoupling method and the proposed method. The results show that the proposed method is valid, while forced decoupling approach can't reflect the characteristics of isolated buildings and may lead to insecurity of structures.