• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.36 no.3
• /
• pp.203-211
• /
• 2023

A Markov chain Monte Carlo (MCMC) simulation is proposed for probabilistic full waveform inversion (FWI) in a layered half-space. Dynamic responses on the half-space surface are estimated using the thin-layer method when a harmonic vertical force is applied. Subsequently, a posterior probability distribution function and the corresponding objective function are formulated to minimize the difference between estimations and observed data as well as that of model parameters from prior information. Based on the gradient of the objective function, a proposal distribution and an acceptance probability for MCMC samples are proposed. The proposed MCMC simulation is applied to several layered half-space examples. It is demonstrated that the proposed MCMC simulation for probabilistic FWI can estimate probabilistic material properties such as the shear-wave velocities of a layered half-space.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.40 no.2
• /
• pp.229-236
• /
• 2016

For this research, the trajectory of a projectile was simulated via the multibody dynamics analysis of a self-propelled mortar. The dynamic model was composed of a mortar model and a vehicle model, and was simulated using the RecurDyn program. Interior ballistic was applied to the mortar model, and exterior ballistic was conducted by Matlab using the simulation results of the interior trajectory. Through repetitive Monte-Carlo simulations, the accuracy of the mortar was analyzed by considering variations in the aiming angle and vehicle dynamic response.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.38 no.7
• /
• pp.647-656
• /
• 2014

To supplement the ISO-GUM method for the evaluation of measurement uncertainty, a simulation program using the Monte Carlo method (MCM) was developed, and the MCM and GUM methods were compared. The results are as follows: (1) Even under a non-normal probability distribution of the measurand, MCM provides an accurate coverage interval; (2) Even if a probability distribution that emerged from combining a few non-normal distributions looks as normal, there are cases in which the actual distribution is not normal and the non-normality can be determined by the probability distribution of the combined variance; and (3) If type-A standard uncertainties are involved in the evaluation of measurement uncertainty, GUM generally offers an under-valued coverage interval. However, this problem can be solved by the Bayesian evaluation of type-A standard uncertainty. In this case, the effective degree of freedom for the combined variance is not required in the evaluation of expanded uncertainty, and the appropriate coverage factor for 95% level of confidence was determined to be 1.96.

• Journal of the Korean Society of Radiology
• /
• v.14 no.5
• /
• pp.503-510
• /
• 2020

Concrete is one of the most widely used materials as the shielding structures of a nuclear facilities. It is also the most generated radioactive waste in quantity while dismantling facilities. Since the concrete captures neutrons and generates various radionuclides, radiation measurement and analysis of the sample was fulfilled prior to dismantle facilities. An HPGe detector is used in general for the radiation measurement, and effective correction factors such as geometrical correction factor, self-absorption correction, and absolute detector efficiency have to be applied to the measured data to decide exact radioactivity of the sample. Correction factors are obtained by measuring data using a standard source with the same geometry and chemical states as the sample under the same measurement conditions. However, it is very difficult to prepare standard concrete sources because concrete is limited in pretreatment due to various constituent materials and high density. In addition, the concrete sample obtained by core drill is a volumetric source, which requires geometric correction for sample diameter and self absorption correction for sample density. Therefore in recent years, many researchers are working on the calculation of effective correction factors using Monte carlo simulation instead of measuring them using a standard source. In this study we calculated, using Geant4, one of the Monte carlo codes, the correction factors for the various diameter and density of the concrete core sample at the gamma ray energy emitted from the nuclides ¹⁵²Eu and ⁶⁰Co, which are the most generated in radioactive concrete.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.40 no.4
• /
• pp.336-345
• /
• 2012

In crack growth life, uncertainties are caused by variance of geometry, applied loads and material properties. Therefore, the reliability estimation for these uncertainties is required to keep the robustness of calculated life. The stress intensity factors are the most important variable in crack growth life calculation, but its equation is hard to know for complex geometry, therefore they are processed by the finite element analysis which takes long time. In this paper, the response surface is considered to increase efficiency of the reliability analysis for crack growth life of a turbine wheel. The approximation model of the stress intensity factors is obtained by the regression analysis for FEA data and the response surface of crack growth life is generated for selected factors. The reliability analysis is operated by the Monte Carlo Simulation for the response surface. The results indicate that the response surface could reduce computations that need for reliability analysis for the turbine wheel, which is hard to derive stress intensity factor equation, successfully.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.29 no.2C
• /
• pp.247-254
• /
• 2004

We introduce a new space-time signal processing scheme that uses both transmitter diversity technique and transmitter beamforming technique for code-division multiple access (CDMA) systems. Over complex Gaussian Rayleigh channel, the introduced scheme achieves the diveristy gain through the transmitter diversity technique. and the SNR gain by th transmitter beamforming technique. Bit error rate (BER) analyses are given to each of the three cases in which the transmitter diversity scheme, the transmitter beamforming scheme and the introduced scheme are used, in the slowly varying Rayleigh frequency nonselective fading channel. The Monte-Carlo simulation results are shown to match to the analytic results. When the channels between distant antennas are independent, analytic results show that the introduced scheme achieves the lowest $E_{b/}$ $N_{0}$ at target BER 10$^{-6}$ . When the channels between distant antennas are correlated, analytic and simulation results show that the introduced scheme is more robust to the change of channel correlation.n.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.2
• /
• pp.153-158
• /
• 2017

In this paper, we take into account topology optimization problems considering spatial randomness in the material property of elastic modulus. Based on 88 lines MATLAB Code, Monte Carlo analysis has been performed for MBB(messerschmidt-$b{\ddot{o}}lkow$-blohm) model using 5,000 random sample fields which are generated by using the spectral representation scheme. The random elastic modulus is assumed to be Gaussian in the spatial domain of the structure. The variability of the volume fraction of the material, which affects the optimum topology of the given problem, is given in terms of correlation distance of the random material. When the correlation distance is small, the randomness in the topology is high and vice versa. As the correlation distance increases, the variability of the volume fraction of the material decreases, which comply with the feature of the linear static analysis. As a consequence, it is suggested that the randomness in the material property is need to be considered in the topology optimization.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.24 no.4
• /
• pp.383-389
• /
• 2011

In this paper, a statistical calibration method is proposed in order to identify the variability of complex modulus for a rubber material due to operational temperature and experimental/model errors. To describe temperature- and frequency-dependent material properties, a fractional derivative model and a shift factor relationship are used. A likelihood function is defined as a product of the probability density functions where experimental values lie on the model. The variation of the fractional derivative model parameters is obtained by maximizing the likelihood function. Using the proposed method, the variability of a synthetic rubber material is estimated and applied to a rubber mount problem. The dynamic characteristics of the rubber mount are calculated using a finite element model of which material properties are sampled from Monte Carlo simulation. The calculated dynamic stiffnesses show very large variation.

• Journal of the Korean Geotechnical Society
• /
• v.34 no.5
• /
• pp.37-51
• /
• 2018

The slope failure due to the rainfall infiltration occurs frequently in Korea, since the depth of the weathered residual soil layer is shallow in mountainous region. Depth of the failure surface is shallow and tends to pass near the interface between impermeable bedrock and soil layer. Soil parameters that have a significant impact on the instability of unsaturated slopes due to rainfall infiltration inevitably include large uncertainties. Therefore, this study proposes a probabilistic analysis procedure by Monte Carlo Simulation which considers the hydraulic characteristics and strength characteristics of soil as random variables in order to predict slope failure due to rainfall infiltration. The Green-Ampt infiltration model was modified to reflect the boundary conditions on the slope surface according to the rainfall intensity and the boundary condition of the shallow impermeable bedrock was introduced to predict the stability of unsaturated soil slope with shallow bedrock under constant rainfall intensity. The results of infiltration analysis were used as inputs of infinite slope analysis to calculate the safety factor. The proposed analysis method can be used to calculate the time-dependent failure probability of soil slope due to rainfall infiltration.

• Proceedings of the KSME Conference
• /
• 2007.05a
• /
• pp.921-926
• /
• 2007

A dynamic system under random disturbance is considered in the study. In order to control the system efficiently, proper reduction of system dimension is indispensible in design stage. The reduction method using component cost analysis in conjunction with stochastic analysis is proposed for the control of a system. System response is obtained in terms of dynamic moment equation via Fokker-Plank-Kolmogorov(F-P-K) equation. The dynamic moment response of the system under random disturbance are reduced by using of deterministic version of component cost analysis. The reduced system via proposed "stochastic component cost analysis" is successfully implemented for dynamic response and shows remarkable control performance effectively utilizing "stochastic controller" in physical time domain.