• Structural Engineering and Mechanics
• /
• v.49 no.1
• /
• pp.111-128
• /
• 2014

The stochastic response surface method (SRSM) and the response surface method (RSM) are often used for structural reliability analysis, especially for reliability problems with implicit performance functions. This paper aims to compare these two methods in terms of fitting the performance function, accuracy and efficiency in estimating probability of failure as well as statistical moments of system output response. The computational procedures of two response surface methods are briefly introduced first. Then their capabilities are demonstrated and compared in detail through two examples. The results indicate that the probability of failure mainly reflects the accuracy of the response surface function (RSF) fitting the performance function in the vicinity of the design point, while the statistical moments of system output response reflect the accuracy of the RSF fitting the performance function in the entire space. In addition, the performance function can be well fitted by the SRSM with an optimal order polynomial chaos expansion both in the entire physical and in the independent standard normal spaces. However, it can be only well fitted by the RSM in the vicinity of the design point. For reliability problems involving random variables with approximate normal distributions, such as normal, lognormal, and Gumbel Max distributions, both the probability of failure and statistical moments of system output response can be accurately estimated by the SRSM, whereas the RSM can only produce the probability of failure with a reasonable accuracy.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.15 no.3
• /
• pp.44-51
• /
• 2016

This paper presents the characteristics of the energy absorption in an expansion tube type impact absorber that is applied to a high weight drop tester and the use of a response surface methodology to predict the impact energy absorption. In order to identify the characteristics of the energy absorption, a set of finite element analysis was conducted with Abaqus Explicit. Moreover, the ISCD-II sampling method and a first order polynomial were used to build a response surface. As a result, we demonstrated that the impact energy could be controlled by four main design variables, namely an expansion pipe's thickness, inner radius, pressing die's expansion angle and expansion ratio. Additionally, we observed the relationship between the four main design variables and the impact energy absorbing time, displacement, and maximum impact force.

• Earthquakes and Structures
• /
• v.8 no.4
• /
• pp.915-934
• /
• 2015

Within the context of Structural Health Monitoring (SHM), it is often the case that structural systems are described by uncertainty, both with respect to their parameters and the characteristics of the input loads. For the purposes of system identification, efficient modeling procedures are of the essence for a fast and reliable computation of structural response while taking these uncertainties into account. In this work, a reduced order metamodeling framework is introduced for the challenging case of nonlinear structural systems subjected to earthquake excitation. The introduced metamodeling method is based on Nonlinear AutoRegressive models with eXogenous input (NARX), able to describe nonlinear dynamics, which are moreover characterized by random parameters utilized for the description of the uncertainty propagation. These random parameters, which include characteristics of the input excitation, are expanded onto a suitably defined finite-dimensional Polynomial Chaos (PC) basis and thus the resulting representation is fully described through a small number of deterministic coefficients of projection. The effectiveness of the proposed PC-NARX method is illustrated through its implementation on the metamodeling of a five-storey shear frame model paradigm for response in the region of plasticity, i.e., outside the commonly addressed linear elastic region. The added contribution of the introduced scheme is the ability of the proposed methodology to incorporate uncertainty into the simulation. The results demonstrate the efficiency of the proposed methodology for accurate prediction and simulation of the numerical model dynamics with a vast reduction of the required computational toll.

• Nuclear Engineering and Technology
• /
• v.54 no.9
• /
• pp.3478-3487
• /
• 2022

In this study, we incorporate an anisotropic scattering scheme involving spherical harmonics into the method of characteristics (MOC). The neutron transport solution in a light water reactor can be significantly improved because of the impact of an anisotropic scattering source with the MOC flat source approximation. Several problems are selected to verify the proposed scheme and investigate its effects and accuracy. The MOC anisotropic scattering source is based on the expansion of spherical harmonics with Legendre polynomial functions. The angular flux, scattering source, and cross section are expanded in terms of the surface spherical harmonics. Later, the polynomial is expanded to achieve the odd and even parity of the source components. Ultimately, the MOC angular and scalar fluxes are calculated from a combination of two sources. This paper presents various numerical examples that represent the hot and cold conditions of a reactor core with boron concentration, burnable absorbers, and control rod materials, with and without a reflector or baffle. Moreover, a small critical core problem is considered which involves significant neutron leakage at room temperature. We demonstrate that an anisotropic scattering source significantly improves solution accuracy for the small core high-leakage problem, as well as for practical large core analyses.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.291-291
• /
• 2021

모형이 갖는 불확실성의 정량화나 매개변수의 최적화는 계산시간의 기하급수적인 증가를 가져온다. 계산시간의 효율성을 극대화할 수 있는 기법으로 최근 대체모형이 개발되었으며, 다양한 분야에서 적용되고 있다. 그러나 대체모형은 훈련된 데이터 공간에서 크게 벗어난 극한 사상를 정확하게 모의하기는 어려운 단점이 있다. 본 연구는 이와 같은 대체모형의 단점을 개선할 수 있는 새로운 PCK(polynomial chaos-krigging) 기법을 제시한다. PCK는 PCE(polynomial chaos expansion) 기법과 OK(ordinary krigging) 기법을 결합한 것이며, PCK의 효과는 기존의 PCE 및 OK 모형의 결과와 비교하여 입증하였다. 본 연구의 분석 결과는 다음과 같다. (1) PCK는 더 적은 수의 훈련 샘플만으로도 원래 모형을 더 정확하게 대체할 수 있다. (2) 원래 훈련 샘플보다 약 3배 더 큰 극한사상을 모의했을 때, PCE와 OK는 예측이 실패하였지만, PCK의 예측은 정확하였다. (3) 민감도 분석 결과 PCK의 매개변수 특성과 거동이 PCE 및 OK보다 원래 모형의 특성과 거동에 더 일치한다. 본 연구에서는 3개의 대체모형의 결과를 원래모형의 결과와 비교하였으며 그 적용성을 극한강우에 대해 검토하였다. 일반적으로 훈련 샘플의 범위와 비슷한 강우사상에 대해서는 모든 대체모형의 결과가 우수하였으나, 훈련 샘플의 범위에서 벗어난 극한 사상의 모의는 PCK만 적용이 가능하였다. 제안된 대체모형은 극한사상의 예측에 있어 기존 대체모형보다 매우 향상된 정확도를 제공함을 확인할 수 있었다.

• Structural Engineering and Mechanics
• /
• v.28 no.2
• /
• pp.129-152
• /
• 2008

This paper considers the solution of the stochastic differential equations (SDEs) with random operator and/or random excitation using the spectral SFEM. The random system parameters (involved in the operator) and the random excitations are modeled as second order stochastic processes defined only by their means and covariance functions. All random fields dealt with in this paper are continuous and do not have known explicit forms dependent on the spatial dimension. This fact makes the usage of the finite element (FE) analysis be difficult. Relying on the spectral properties of the covariance function, the Karhunen-Loeve expansion is used to represent these processes to overcome this difficulty. Then, a spectral approximation for the stochastic response (solution) of the SDE is obtained based on the implementation of the concept of generalized inverse defined by the Neumann expansion. This leads to an explicit expression for the solution process as a multivariate polynomial functional of a set of uncorrelated random variables that enables us to compute the statistical moments of the solution vector. To check the validity of this method, two applications are introduced which are, randomly loaded simply supported reinforced concrete beam and reinforced concrete cantilever beam with random bending rigidity. Finally, a more general application, randomly loaded simply supported reinforced concrete beam with random bending rigidity, is presented to illustrate the method.

• The Transactions of the Korean Institute of Electrical Engineers C
• /
• v.53 no.12
• /
• pp.626-633
• /
• 2004

In this paper, we present a time domain combined field integral equation (TD-CFIE) formulation to analyze the transient electromagnetic response from three-dimensional dielectric objects. The solution method in this paper is based on the method of moments (MoM) that involves separate spatial and temporal testing procedures. A set of the RWG (Rao, Wilton, Glisson) functions Is used for spatial expansion of the equivalent electric and magnetic current densities and a combination of RWG and its orthogonal component is used as spatial testing. We also investigate spatial testing procedures for the TD-CFIE to select the proper testing functions that are derived from the Laguerre polynomials. These basis functions are also used for temporal testing. Use of this temporal expansion function characterizing the time variable enables one to handle the time derivative terms in the integral equation and decouples the space-time continuum in an analytic fashion. Numerical results computed by the proposed formulation are presented and compared with the solutions of the frequency domain combined field integral equation (FD-CFIE).

• Journal of the Institute of Electronics Engineers of Korea SD
• /
• v.39 no.2
• /
• pp.27-38
• /
• 2002

A new self-consistent mathematical model for semiconductor quantum well device was developed. The model was based on the direct solution of the Boltzmann transport equation, coupled to the Schrodinger and Poisson equations. The solution yielded the distribution function for a two-dimensional electron gas(2DEG) in quantum well devices. To solve the Boltzmann equation, it was transformed into a tractable form using a Legendre polynomial expansion. The Legendre expansion facilitated analytical evaluation of the collision integral, and allowed for a reduction of the dimensionality of the problem. The transformed Boltzmann equation was then discretized and solved using sparce matrix algebra. The overall system was solved by iteration between Poisson, Schrodinger and Boltzmann equations until convergence was attained.

• ETRI Journal
• /
• v.29 no.1
• /
• pp.8-17
• /
• 2007

In this paper, we present a time domain combined field integral equation formulation (TD-CFIE) to analyze the transient electromagnetic response from dielectric objects. The solution method is based on the method of moments which involves separate spatial and temporal testing procedures. A set of the RWG functions is used for spatial expansion of the equivalent electric and magnetic current densities, and a combination of RWG and its orthogonal component is used for spatial testing. The time domain unknowns are approximated by a set of orthonormal basis functions derived from the Laguerre polynomials. These basis functions are also used for temporal testing. Use of this temporal expansion function characterizing the time variable makes it possible to handle the time derivative terms in the integral equation and decouples the space-time continuum in an analytic fashion. Numerical results computed by the proposed formulation are compared with the solutions of the frequency domain combined field integral equation.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.16 no.2
• /
• pp.211-221
• /
• 2018

Nodal transport methods are proposed for solving the simplified even-parity neutron transport (SEP) equation. These new methods are attributed to the success of existing nodal diffusion methods such as the Polynomial Expansion Nodal and the Analytic Function Expansion Nodal Methods, which are known to be very effective for solving the neutron diffusion equation. Numerical results show that the simplified even-parity transport equation is a valid approximation to the transport equation and that the two nodal methods developed in this study also work for the SEP transport equation, without conflict. Since accuracy of methods is easily increased by adding node unknowns, the proposed methods will be effective for coarse mesh calculation and this will also lead to computation efficiency.