• Nuclear Engineering and Technology
• /
• v.25 no.1
• /
• pp.51-62
• /
• 1993

A computational program ［TDET］ of the particle transport equation is developed on radiation shielding problem in two-dimensional cartesian geometry based on the discrete element method. Not like the ordinary discrete ordinates method, the quadrature set of angles is not fixed but steered by the spatially dependent angular fluxes. The angular dependence of the scattering source term in the particle transport equation is described by series expansion in spherical harmonics, and the energy dependence of the particles is considered as well. Three different benchmark tests are made for verification of TDET : For the ray effect analysis on a square absorber with a flat isotropic source, the results of TDET calculation are quite well conformed to those of MORSE-CG calculation while TDET ameliorates the ray effect more effectively than S$_{N}$ calculation. In the analysis of the streaming leakage through a narrow vacuum duct in a shield, TDET shows conspicuous and remarkable results of streaming leakage through the duct as well as MORSE-CG does, and quite better than S$_{N}$ calculation. In a realistic reactor shielding situation which treats in two cases of the isotropic scattering and of linearly anisotropic scattering with two groups of energy, TDET calculations show local ray effect between neighboring meshes compared with S$_{N}$ calculations in which the ray effect extends broadly over several meshes.eshes.

• International Journal of Highway Engineering
• /
• v.17 no.4
• /
• pp.25-31
• /
• 2015

PURPOSES: Implementation and verification of the simple linear cohesive viscoelastic contact model that can be used to simulate dynamic behavior of sticky aggregates. METHODS: The differential equations were derived and the initial conditions were determined to simulate a free falling ball with a sticky surface from a ground. To describe this behavior, a combination of linear contact model and a cohesive contact model was used. The general solution for the differential equation was used to verify the implemented linear cohesive viscoelastic API model in the DEM. Sensitivity analysis was also performed using the derived analytical solutions for several combinations of damping coefficients and cohesive coefficients. RESULTS : The numerical solution obtained using the DEM showed good agreement with the analytical solution for two extreme conditions. It was observed that the linear cohesive model can be successfully implemented with a linear spring in the DEM API for dynamic analysis of the aggregates. CONCLUSIONS: It can be concluded that the derived closed form solutions are applicable for the analysis of the rebounding behavior of sticky particles, and for verification of the implemented API model in the DEM. The assumption of underdamped condition for the viscous behavior of the particles seems to be reasonable. Several factors have to be additionally identified in order to develop an enhanced contact model for an asphalt mixture.

• Structural Engineering and Mechanics
• /
• v.35 no.1
• /
• pp.19-35
• /
• 2010

This paper proposes a hybrid heuristic and criteria-based method of optimum design which combines the advantages of both the iterated simulated annealing (SA) algorithm and the rigorously derived optimality criteria (OC) for structural optimum design of reinforced concrete (RC) buildings under multi-load cases based on the current Chinese design codes. The entire optimum design procedure is divided into two parts: strength optimum design and stiffness optimum design. A modified SA with the strategy of adaptive feasible region is proposed to perform the discrete optimization of RC frame structures under the strength constraints. The optimum stiffness design is conducted using OC method with the optimum results of strength optimum design as the lower bounds of member size. The proposed method is integrated into the commercial software packages for building structural design, SATWE, and for finite element analysis, ANSYS, for practical applications. Finally, two practical frame-shear-wall structures (15-story and 30-story) are optimized to illustrate the effectiveness and practicality of the proposed optimum design method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.29 no.9 s.240
• /
• pp.1262-1269
• /
• 2005

Thrust of linear motor is one of the important factor to specify motor performance. Maximum thrust can be obtained by increasing the current in conductor and is relative to the sizes of conductor and magnet. But, the current and the size of conductor have an effect on temperature of linear motor. Therefore, it is practically important to find design results that can effectively maximize the thrust of linear motor within limited range of temperature. Finite element analysis was applied to calculate thrust and the temperature of the conductor was calculated by the thermal resistance. The diameter of copper wire among design variables has discrete value and number of turns must be integer. Considering these facts, special techinque for optimum design is presented. To reduce excessive computation time of thrust in optimization, the design variables was redefined by analysis of variance and second order regression model for thrust was determined by response surface metheodology. As a result, it is shown that the proposed method has an advantage in optimum design of linear motor.

• Current Optics and Photonics
• /
• v.7 no.5
• /
• pp.529-536
• /
• 2023

The Morlet wavelet transform method is proposed to analyze a single interferogram with spatial carrier frequency that is captured by an optical interferometer. The method can retain low frequency components that contain the phase information of a measured optical surface, and remove high frequency disturbances by wavelet decomposition and reconstruction. The key to retrieving the phases from the low-frequency wavelet components is to extract wavelet ridges by calculating the maximum value of the wavelet transform amplitude. Afterwards, the wrapped phases can be accurately solved by multiple iterative calculations on wavelet ridges. Finally, we can reconstruct the wave-front of the measured optical element by applying two-dimensional discrete cosine transform to those wrapped phases. Morlet wavelet transform does not need to remove the spatial carrier frequency components manually in the processing of interferogram analysis, but the step is necessary in the Fourier transform algorithm. So, the Morlet wavelet simplifies the process of the analysis of interference fringe patterns compared to Fourier transform. Consequently, wavelet transform is more suitable for automated programming analysis of interference fringes and avoiding the introduction of additional errors compared with Fourier transform.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.21 no.6
• /
• pp.35-40
• /
• 2020

The ballast layers play a key role in distributing and supporting a trainload. On the other hand, it settles down by dynamic train loading due to large void ratios. Consequently, it requires continuous maintenance. In this paper, ballast layers mixed with three types of specially shaped artificial ballast (AB) (Rectangular, Tetrapod, Hexagonal) were modeled by using a two dimensional DEM (Discrete Element Method). Repeated loading tests were performed to evaluate the settlement behavior of the ballast layers. The smallest settlement was observed in the case of the ballast layer mixed with Tetrapod AB than in other cases, according to an analysis of the force transfer routes. In addition, contact force analysis showed that the Tetrapod AB, which has a concave shape, could easily make small and multi-channel force-transfer routes. This means that the stress in the ballast layer by the train loading transferred through the sleeper uniformly was distributed well by the AB. Therefore, the settlement of the ballast layer mixed with the concave-shaped Tetrapod AB could be reduced effectively under a repeated train loading. The effects of a decrease in settlement of the ballast layer highlight the possibility of a maintenance-free ballasted track.

• Journal of the Korean Geotechnical Society
• /
• v.30 no.9
• /
• pp.41-55
• /
• 2014

In performing seismic analysis of tunnels, it is a common practice to ignore the rock joints and to assume that the rock mass surrounding the tunnel is continuous. The applicability of this assumption has not yet been validated in detail. This study performs a series of pseudo-static discrete element analyses to evaluate the effect of rock joint on the seismic response of tunnels. The parameters considered are joint intersection location, joint spacing, joint stiffness, joint dip, and interface stiffness. The results show that the joint stiffness has the most critical influence on the tunnel response. The tunnel response increases with the spacing, resulting in localized concentration of moment and shear stress. The response of the tunnel is the lowest for joints dipping at $45^{\circ}$. This is because large shear stresses result in rotation of the principal planes by $45^{\circ}$. In summary, the weathered and smooth, vertical or horizontal, and widely spaced joint set will significantly increase the tunnel response under seismic loading. The tunnel linings are shown to be most susceptible to damage due to induced shear stress, and therefore should be checked in the seismic design.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.1
• /
• pp.90-100
• /
• 1993

A finite element formulation of vibration control of a laminated plate with piezoelectric sensor/ actuators is presented. Classical lamination theory with the induced strain actuation and Hamilton's principle are used to formulate the equations of motion of the system. The total charge developed on the sensor layer is calculated from the direct piezoelectric equation. The equations of motion and the total charge are discretized with 4 node, 12 degrees of freedom quadrilateral plate bending elements with one electrical degree of freedom. The mass and stiffness of the piezoelectric layer are introduced by treating them as another layer in laminated plate. Piezoelectric sensor/actuators are distributed, but discrete due to the geometry of electrodes. By defining an i.d. number of electrode for each element, modelling of electrodes with variable geometry can be achieved. The static response of a piezoelectric bimorph beam to electrical loading and sensor voltage to given displacement are calculated. For a laminated plate under the negative velocity feedback control, the direct time response by the Newmark-.betha. method and damped frequencies and modal damping ratios by modal state space analysis are derived.

• International Journal of Highway Engineering
• /
• v.16 no.2
• /
• pp.43-52
• /
• 2014

PURPOSES : To solve problems in current compaction control DCPT(Dynamic Cone Penetrometer Test), highly correlated with various testing methods, simple, and economic is being applied. However, it、s hard to utilize DCPT results due to the few numerical analyses for DCPT have been performed and the lack of data accumulation. Therefore, this study tried to verify the validation of numerical modeling for DCPT by comparing and analyzing the results of numerical analyses with field tests. METHODS: The ground elastic modulus and PR(Penetration Rate) value were estimated by using PFC(Particle Flow Code) 3D program based on the discrete element method. Those values were compared and analyzed with the result of field tests. Also, back analysis was conducted to describe ground elastic modulus of field tests. RESULTS : Relative errors of PR value between the numerical analyses and field tests were calculated to be comparatively low. Also, the relationship between elastic modulus and PR value turned out to be similar. CONCLUSIONS : Numerical modeling of DCPT is considered to be suitable for describing field tests by carrying out numerical analysis using PFC 3D program.

• Geomechanics and Engineering
• /
• v.12 no.4
• /
• pp.639-655
• /
• 2017

The binary mixture consists of two types of granular media with different physical attributes and sizes, which can be characterized by the percentage of large granules by weight (P) and the particle size ratio (${\alpha}$). Researchers determine that two thresholds ($P_S$ and $P_L$) exist for the peak shear strength of binary mixtures, i.e., at $P{\leq}P_S$, the peak shear strength is controlled by the small granules; at $P{\leq}P_L$, the peak shear strength is controlled by the large granules; at $P_S{\leq}P{\leq}P_L$, the peak shear strength is governed by both the large and small granules. However, the thresholds of binary mixtures with different ${\alpha}$ values, and the explanation related to the inner details of binary mixtures to account for why these thresholds exist, require further confirmation. This paper considers the mechanical behavior of binary mixtures with DEM analysis. The thresholds of binary mixtures are found to be strongly related to their coordination numbers $Z_L$ for all values of ${\alpha}$, where $Z_L$ denotes the partial coordination number only between the large particles. The arrangement structure of the large particles is examined when P approaches the thresholds, and a similar arrangement structure of large particles is formed in both 2D and 3D particle systems.