• Proceedings of the Korean Geotechical Society Conference
• /
• 2002.03a
• /
• pp.651-656
• /
• 2002

In general, the most important problem in slope stability analysis is that there is no definite way to describe the natural three-dimensional Joint network. Therefore, the many approaches were tried to anlayze the slope stability. Numerical modeling approach is one of the branch to resolve the complexity of natural system. UDEC, FLAC, and SWEDGE are widely used commercial code for the purpose on stability analysis. For the purpose on the more appropriate application of these kind of code, however, three-dimensional distribution of joint network must be identified in more explicit way. Remaining problem is to definitely describe the three dimensional network of joint and bedding, but it is almost impossible in practical sense. Three dimensional joint generation method with random number generation and the results of generation to UDEC have been applied to settle the refered problems in field site. However, this approach also has a important problem, and it is that joint network is generated only once. This problem lead to the limitation on the application to field case, in practical sense. To get rid of this limitation, Monte Carlo Simulation is proposed in this study 1) statistical analysis of input values and definition of the applied system with statistical parameter, 2) instead of the consideration of generated network as a real system, generated system is just taken as one reliable system, 3) present the design parameters, through the statistical analysis of ouput values Results of this study are not only the probability of failure, but also area of failure block, shear strength, normal strength and failure pattern, and all of these results are described in statistical parameters. The results of this study, shear strength, failure area, pattern etc, can provide the direct basement on the design, cutoff angle, support pattern, support strength and etc.

• Computational Structural Engineering
• /
• v.4 no.2
• /
• pp.111-117
• /
• 1991

The analysis method calculating the mean and standard deviation for the eigenvalue of complicated structures in which the limit state equation is implicitly expressed is formulated and applied to the buckling analysis by combining probabilistic finite element method with direct differential method which is a kind of sensitivity analysis technique. Also, the probability of buckling failure is calculated by combining classical reliability techniques such a MVFOSM and AFOSM. As random variables external load, elastic modulus, sectional moment of inertia and member length are chosen and Parkinson's iteration algorithm in AFOSM is used. The accuracy of the results by this study is verified by comparing the results with the crude Monte Carlo simulation and Importance Sampling Method. Through the case study of some structures the important aspects of buckling reliability analysis are discussed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.10
• /
• pp.1406-1418
• /
• 2002

The direct simulation Monte Carlo(DSMC) method is employed to calculate the etch rate on Al wafer. The etchant is assumed to be Cl$_2$. The etching process of an Al wafer in a helicon plasma etcher is examined by simulating molecular collisions of reactant and product. The flow field inside a plasma etch reactor is also simulated by the DSMC method fur a chlorine feed gas flow. The surface reaction on the Al wafer is simply modelled by one-step reaction: 3C1$_2$+2Allongrightarrow1 2AIC1$_3$. The gas flow inside the reactor is compared for six different nozzle locations. It is found that the flow field inside the reactor is affected by the nozzle locations. The Cl$_2$ number density on the wafer decreases as the nozzle location moves toward the side of the reactor. Also, the present numerical results show that the nozzle location 1, which is at the top of the reactor chamber, produces a higher etch rate.

• Journal of Energy Engineering
• /
• v.14 no.1
• /
• pp.46-53
• /
• 2005

This paper presents a new approach for determining an accurate incentive levels of Direct Load Control (DLC) program using probabilistic techniques. The economic analysis of DLC resources needs to identify the hourly-by-hourly expected energy-not-served resulting from the random outage characteristics of generators as well as to reflect the availability and duration of DLC resources, which results the computational explosion. Therefore, the conventional methods are based on the scenario approaches to reduce the computation time as well as to avoid the complexity of economic studies. In this paper, we have developed a new technique based on the sequential Monte Carlo simulation to evaluate the required expected load control amount in each hour and to decide the incentive level satisfying the economic constraints. In addition, we have applied the variance reduction technique to enhance the efficiency of the simulation. To show the efficiency and effectiveness of the suggested method, the numerical studies have been performed for the modified IEEE 24-bus reliability test system.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.2
• /
• pp.1-6
• /
• 2004

• Nuclear Engineering and Technology
• /
• v.53 no.12
• /
• pp.4142-4149
• /
• 2021

In the radiation protection application, the metal-polymer composites have been developed for their radiation shielding properties. In this research, the elastomer composites doped by 10 ㎛ and 100nm size of lead, bismuth and tungsten particles as filler with 30 and 60 wt percentages were prepared. To survey the shielding properties of the polymer composites using gamma-ray emitted from 152Eu and 137Cs sources, the gamma flux was measured by using NaI(Tl) detector, then the linear attenuation coefficient was calculated. Also, the Monte Carlo simulation (MCs) method was used. The results showed a direct relationship between the linear attenuation coefficients of the absorbent and filler ratio. Also, the decrease in the particle size of the shielding material in each weight percentage improved the radiation shielding features. When the dimension of the particles was in the order of nano-size, more attenuation was achieved. At low energies used for medical diagnostic X-ray applications due to the predominance of the photoelectric effect, bismuth and lead were suitable selection as filler.

• Journal of Electrical Engineering and Technology
• /
• v.11 no.3
• /
• pp.602-608
• /
• 2016

Because of the random characteristics of lightning, the Monte Carlo method is applied to estimate the flashover rate due to lightning, however, the simulations using previous methods are difficult to both beginner and expert in power corporations. Therefore, this paper proposes the new and easy method to simulate the flashover rate by connection of electromagnetic transients program (EMTP) and MATLAB. The magnitude of a lightning strike is based on a curve measured in the field, while the classification of direct and indirect lightning depends on the striking distance. In a Korean distribution system, the flashover rate induced by lightning is simulated using proposed method. Simulations of the footing resistance according to the existence of an overhead ground wire (OHGW) are performed and the simulation results are discussed. The simulation results are compared with findings obtained with the IEEE Flash 2.0 program.

• Advances in aircraft and spacecraft science
• /
• v.4 no.5
• /
• pp.585-611
• /
• 2017

This article covenants with the post buckling witticism of carbon nanotube reinforced composite (CNTRC) beam supported with an elastic foundation in thermal atmospheres with arbitrary assumed random system properties. The arbitrary assumed random system properties are be modeled as uncorrelated Gaussian random input variables. Unvaryingly distributed (UD) and functionally graded (FG) distributions of the carbon nanotube are deliberated. The material belongings of CNTRC beam are presumed to be graded in the beam depth way and appraised through a micromechanical exemplary. The basic equations of a CNTRC beam are imitative constructed on a higher order shear deformation beam (HSDT) theory with von-Karman type nonlinearity. The beam is supported by two parameters Pasternak elastic foundation with Winkler cubic nonlinearity. The thermal dominance is involved in the material properties of CNTRC beam is foreseen to be temperature dependent (TD). The first and second order perturbation method (SOPT) and Monte Carlo sampling (MCS) by way of CO nonlinear finite element method (FEM) through direct iterative way are offered to observe the mean, coefficient of variation (COV) and probability distribution function (PDF) of critical post buckling load. Archetypal outcomes are presented for the volume fraction of CNTRC, slenderness ratios, boundary conditions, underpinning parameters, amplitude ratios, temperature reliant and sovereign random material properties with arbitrary system properties. The present defined tactic is corroborated with the results available in the literature and by employing MCS.

• Structural Engineering and Mechanics
• /
• v.64 no.4
• /
• pp.413-426
• /
• 2017

Probabilistic methods are used in engineering where a computational model contains random variables. The proposed method under development: Direct Optimized Probabilistic Calculation (DOProC) is highly efficient in terms of computation time and solution accuracy and is mostly faster than in case of other standard probabilistic methods. The novelty of the DOProC lies in an optimized numerical integration that easily handles both correlated and statistically independent random variables and does not require any simulation or approximation technique. DOProC is demonstrated by a collection of deliberately selected simple examples (i) to illustrate the efficiency of individual optimization levels and (ii) to verify it against other highly regarded probabilistic methods (e.g., Monte Carlo). Efficiency and other benefits of the proposed method are grounded on a comparative case study carried out using both the DOProC and MC techniques. The algorithm has been implemented in mentioned software applications, and has been used effectively several times in solving probabilistic tasks and in probabilistic reliability assessment of structures. The article summarizes the principles of this method and demonstrates its basic possibilities on simple examples. The paper presents unpublished details of probabilistic computations based on this method, including a reliability assessment, which provides the user with the probability of failure affected by statistically dependent input random variables. The study also mentions the potential of the optimization procedures under development, including an analysis of their effectiveness on the example of the reliability assessment of a slender column.

• Advances in aircraft and spacecraft science
• /
• v.4 no.5
• /
• pp.503-514
• /
• 2017

The attitude control of an aircraft is usually fulfilled by means of thrusters at high altitudes. Therefore, the possibility of using also aerodynamic surfaces would produce the advantage of reducing the amount of fuel for the thrusters to be loaded on board. For this purpose, Zuppardi already considered some aerodynamic problems linked to the use of a wing flap in a previous paper. A NACA 0010 airfoil with a trailing edge flap of 35% of the chord, in the range of angle of attack 0-40 deg and flap deflections up to 30 deg was investigated. Computer tests were carried out in hypersonic, rarefied flow by a direct simulation Monte Carlo code at the altitudes of 65 and 85 km of Earth Atmosphere. The present work continues this subject, considering the same airfoil and free stream conditions but two flap extensions of 45% and 25% of the chord and two flap deflections of 15 and 30 deg. The main purpose is to compare the influence of the flap dimension with that of the flap deflection. The present analysis is carried out in terms of: 1) percentage variation of the global aerodynamic coefficients with respect to the no-flap configuration, 2) increment of pressure and heat flux on the airfoil lower surface due to the Shock Wave-Shock Wave Interaction (SWSWI) with respect to the same quantities with no SWSWI or in no-flap configuration, 3) flap hinge moment. Issues 2) and 3) are important for the design of the mechanical and thermal protection system and of the flap actuator, respectively. Under the above mentioned test and geometrical conditions, the flap deflection is aerodynamically more effective than the flap extension, because it involves higher variation of the aerodynamic coefficients. However, tests verify that a smaller deflection angle involves the advantage of a smaller increment of pressure and heat flux on the airfoil lower surface, due to SWSWI, as well as a smaller hinge moment.