• Tunnel and Underground Space
• /
• v.24 no.2
• /
• pp.176-186
• /
• 2014

Roughness, aperture and filling material of rock joint are widely considered to affect the hydraulic characteristics of joint. Among these factors, in this study, the joint roughness was examined with artificial joint profiles generated by Monte Carlo simulating on the original profiles suggested by Barton and Choubey(1977). Original profiles and revised profiles were combined to establish flow channel models, in which the hydraulic characteristics were analyzed numerically on the basis of minimum aperture changes and flow channel shapes. Maximum flow rate was identified at the growing point of flow area after passing through minimum aperture generated by the two profiles, and it was resulted that maximum flow rate is inversely proportional to minimum aperture. Maximum flow rate per unit area showed different values because flow channel shapes and minimum aperture locations are different in each model. In flow channel, mechanical aperture showed approximately 1.07 ~ 3.00 times larger than hydraulic aperture. In this study, mechanical and hydraulic aperture were concluded to be closely related to $A_i$ value, and their relations can be denoted by $e_m=0.519A^{0.7169_i}$ and $e_h=0.6182A^{0.239}_i$, respectively.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.34 no.9
• /
• pp.1305-1310
• /
• 2010

In a practical design process, the method of extracting the design space information of the complex system for verifying, improving, and optimizing the design process by taking into account the design variables and their shape tolerance is very important. Finite element analysis has been successfully implemented and integrated with design of experiment such as D-Optimal array; thus, a response surface model and optimization tools have been obtained, and design variables can be optimized by using the model and these tools. Then, to guarantee the robustness of the design variables, a robust design should be additionally performed by taking into account the statistical variation of the shape tolerance of the optimized design variables. In this study, a new approach based on the use of the response surface model is proposed; in this approach, the standard normal distribution of the shape tolerance is considered. By adopting this approach, it is possible to simultaneously optimize variables and perform a robust design. This approach can serve as a means of efficiently modeling the trade-off among many conflicting goals in the applications of finite element analysis. A case study on the robust optimal design of disc brakes under thermal loadings was carried out to solve multiple objective functions and determine the constraints of the design variables, such as a thermal deformation and weight.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.3
• /
• pp.373-384
• /
• 2000

To assess the safety of nonlinear steel frame structures subjected to short duration dynamic loadings, especially seismic loading, a nonlinear time domain reliability analysis procedure is proposed in the context of the stochastic finite element concept. In the proposed algorithm, the finite element formulation is combined with concepts of the response surface method, the first order reliability method, and the iterative linear interpolation scheme. This leads to the stochastic finite element concept. Actual earthquake loading time-histories are used to excite structures, enabling a realistic representation of the loading conditions. The assumed stress-based finite element formulation is used to increase its efficiency. The algorithm also has the potential to evaluate the risk associated with any linear or nonlinear structure that can be represented by a finite element algorithm subjected to seismic loading or any short duration dynamic loading. The algorithm is explained with help of an example and verified using the Monte Carlo simulation technique.

• Journal of Advanced Navigation Technology
• /
• v.8 no.1
• /
• pp.26-37
• /
• 2004

Consistent and realistic error covariance information is important for position estimation, error analysis, fault detection, and integer ambiguity resolution for differential GNSS. In designing a position domain carrier-smoothed-code filter where incremental carrier phases are used for time-propagation, formulation of consistent error covariance information is not easy due to being bounded and temporal correlation of propagation noises. To provide consistent and correct error covariance information, this paper proposes two recursive filter algorithms based on carrier-smoothed-code techniques: (a) the stepwise optimal position projection filter and (b) the stepwise unbiased position projection filter. A Monte-Carlo simulation result shows that the proposed filter algorithms actually generate consistent error covariance information and the neglection of carrier phase noise induces optimistic error covariance information. It is also shown that the stepwise unbiased position projection filter is attractive since its performance is good and its computational burden is moderate.

• Journal of the Society of Naval Architects of Korea
• /
• v.57 no.2
• /
• pp.80-87
• /
• 2020

As the Arctic sea ice decreases due to various reasons such as global warming, the demand for ships and offshore structures operating in the Arctic region is steadily increasing. In the case of sea ice, the anisotropy is caused by the uncertainty inside the material. For most of the research, nevertheless, estimating the ice load has been treated deterministically. With regard to this, in this paper, a four-point bending strength analysis of an ice specimen was attempted using a stochastic finite element method. First, spatial distribution of the material properties used in the yield criterion was assumed to be a multivariate Gaussian random field. After that, a direct method, which is a sort of stochastic finite element method, and a sensitivity method using the sensitivity of response for random variables were proposed for calculating the probabilistic distribution of ice specimen strength. A parametric study was conducted with different mean vectors and correlation lengths for each material property used in the above procedure. The calculation time was about ten seconds for the direct method and about three minutes for the sensitivity methods. As the cohesion and correlation length increased, the mean value of the critical load and the standard deviation increased. On the contrary, they decreased as the friction angle increased. Also, in all cases, the direct and sensitivity methods yielded very similar results.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.37 no.6
• /
• pp.821-826
• /
• 2013

When the values of design variables change, mode switching can often occur. If the mode of interest is not tracked, the frequencies and modes for design optimization may be miscalculated owing to modes that differ from the intended ones. Thus, mode tracking must be employed to identify the frequencies and modes of interest whenever the values of design variables change during optimization. Furthermore, reliability-based design optimization (RBDO) must be performed for design problems with design variables containing uncertainty. In this research, we perform RBDO considering the mode tracking of a compressive coil spring, i.e., a component of the joint spring that supports a compressor, with design variables containing uncertainty by using only kriging metamodels based on multiple responses approach (MRA) without existing mode tracking methods. The reliability analyses for RBDO are employed using kriging metamodel-based Monte Carlo simulation.

• Journal of the Korean Electrochemical Society
• /
• v.5 no.2
• /
• pp.79-85
• /
• 2002

The present article is concerned with the application of the Monte Carlo simulation to electrochemistry of lithium intercalation from the thermodynamic view point. This article first introduced the fundamental concepts of the ensembles, and Ising and lattice gas models in statistical thermodynamics for the Monte Carlo simulation in brief. Finally the Monte Carlo method based upon the lattice gas model was employed to analyse thermodynamics of the lithium intercalation into the transition metal oxides. Especially we dealt with the thermodynamic properties as the electrode potential curve and the partial molar internal energy and entropy of lithium ion in the case of the $LiMn_2O_4$ electrode, and consequently confirmed the utility of the Monte Carlo method in the field of electrochemistry of the lithium intercalation.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.31 no.10
• /
• pp.10-20
• /
• 2003

In this study, the effective properties were numerically calculated for laminated plain weave textile composites with arbitrary s tacking orientation angles. A single-field macroelement with modified sub-domain integration was used in the analysis to reduce computer resource requirement while efficiently accounting for the internal microstructure. A sample calculation procedure based on the Monte Carlo method was employed to consider the random shift between the layers. Results showed that a significant deviation occurred when the orientation angles were near 0 deg for extensional modulus and Poisson's ratio and 45 deg for the shear modulus. It was also found that the average properties calculated by the 2-layer numerical specimen had large differences compared to the CLT results, which indicated that a caution must be needed when designig of thin plain weave composite structures.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.1
• /
• pp.147-158
• /
• 2000

For the large scale reinforced concrete cooling tower shells, the shape imperfection can be introduced due not only to mistakes in the process of construction but also to the long term behavior of concrete. The shape imperfection evokes the additional responses such as displacements and stresses in addition to the design values. In this study, the statistical behavior of the RC cooling tower shell due to the shape imperfection is investigated using the Monte Carlo simulation. The radius of cooling tower and the shell thickness are adopted as the parameters which cause the shape imperfection. The shape imperfection is modeled as a stochastic field rather than the local one of axisymmetric or bulge type of imperfection. The randomness in the radius is shown to be more affecting the structural responses than the randomness in the shell thickness. In addition to the geometrical randomness, the effect of randomness in the modulus of elasticity on the structural response is also investigated and compared with that of the geometrical ones.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.18 no.9
• /
• pp.1077-1082
• /
• 2007

This paper presents the analysis of the potential ultra-wideband(UWB) interference to WiMAX at 3.5 GHz bands and the mitigation method using transmit power control(TPC) of UWB system. UWB interference effect is evaluated with WiMax's outage probability over UWB density when multiple UWB systems and single WiMAX receiver distribute in unit area of 1$km^2$, When UWB distribution density is 20$devices/km^2$ and the dynamic range of TPC is 30 dB, UWB interference effect with TPC is decreased by 42 % rather than that without mitigation scheme. Finally, we describe that the proposed TPC is an effective method to mitigate UWB interference to WiMAX.