• Nuclear Engineering and Technology
• /
• 제25권1호
• /
• pp.51-62
• /
• 1993

입자수송방정식의 각분할해석법에서 입자흐름방향을 묘사하는 각방향에 대한 각분할집합이 고정된 값이 아니고 각입자속으로 조종되는 각요소법에 근거하여 방사선차폐해석 목적의 전산프로그램(TDET)를 2차원직각좌표계에 대해서 개발하였다. 산란선원항의 각의존성은 Spherical Harmonics Series Expansion으로 해석하였고 입자의 에너지 의존성은 다중에너지군으로 처리하였다. 3 종류의 Benchmark 시험을 통해서 TDET 프로그램을 검증하였다. 평판형 등방적 선원을 가진 사방흡수체에 대한 해석에서 TDET 해석결과가 MORSE-CG 해석결과에 잘 일치하고. 각분할법으로 해석할때 나타나는 Ray effect를 DOT 4.3 보다 잘 치 유하고 있다. 차폐체 내에 좁은 Vacuum duct가 있는 문제의 해석에서 TDET는 MORSE-CG와 마찬가지로 duct를 통한 streaming leak-age에 대해서 분명하고 현저한, 그리고 DOT 4.3보다 매우 훌륭한 해석결과를 보여주고 있다. 원자로차폐구조물 규모에 대해, 2 에너지군, 등방적 산란 및 선형 비등방 산란의 경우에 대해서 해석한 결과 각분할법으로 계산하여 제시하고 있는 기준치에서는 여러 mesh들에 걸쳐 넓은 규모의 Ray effect를 나타내고 있는데, TDET에 의한 해석결과에서는 이웃 mesh들 간에 미세한 Ray effect를 나타내고 있다.

• 한국도로학회논문집
• /
• 제17권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
• /
• 제35권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.

• 대한기계학회논문집A
• /
• 제29권9호
• /
• 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
• /
• 제7권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.

• 한국산학기술학회논문지
• /
• 제21권6호
• /
• pp.35-40
• /
• 2020

도상 자갈층은 열차하중을 분산하고 지지하는 중요한 역할을 수행하나, 큰 공극으로 인하여 열차 동하중 작용시 침하가 발생하여 지속적인 유지보수가 필요하다. 본 논문에서는 특정한 형상의 인공자갈 3종류 (정사각 형상, 테트라 포드 형상, 육각 형상) 를 혼합한 도상 자갈층을 2 차원 개별 요소 해석법을 이용하여 모델링하고, 열차 동하중에 의한 도상 층의 침하 특성을 평가하기 위하여 침목에 의한 반복 재하 시험을 실시하였다. 동 해석 결과 침목 하부에서 작은 힘의 전달경로가 형성된 테트라 포드형의 인공 자갈이 포함된 도상 층에서 침하가 가장 작게 발생하였다. 또한, 도상층 내 접촉력 분포로부터 오목 형상을 갖는 테트라 포드 형상을 포함한 인공 자갈 혼합 도상에서 힘의 다중 전달 경로를 형성하기 쉬우며, 이를 통하여 도상 층으로 전달된 열차 하중이 보다 다양한 경로로 균등하게 분산될 수 있는 것을 확인하였다. 이와 같은 결과로부터 반복 열차 하중 재하 시 오목 형상의 인공 자갈이 혼합된 도상 층에서 훨씬 효과적으로 침하 발생을 저감시킬 수 있다는 것을 알 수 있었다. 도상 층에서의 이와 같은 침하 저감 효과는 유지보수를 절감할 수 있는 자갈궤도 구조에 대한 가능성을 보여준다.

• 한국지반공학회논문집
• /
• 제30권9호
• /
• pp.41-55
• /
• 2014

일반적으로 암반 터널 지진해석 시 절리의 영향을 무시하며 암반을 연속체라 가정한다. 하지만 이와 같은 가정의 적절성에 관해서는 연구된 바 없다. 본 연구에서는 응답변위법을 적용하여 지진동에 의하여 유발되는 암반 절리의 변위가 터널의 응답에 미치는 영향을 일련의 개별요소법을 이용하여 수치적으로 규명하였다. 해석에는 절리와 터널의 교차위치, 절리의 강성, 경사, 간격, 그리고 암반과 터널 사이의 인터페이스 강성이 고려되었다. 해석 결과, 절리의 강성이 터널에 가장 지배적인 영향을 미치는 것으로 나타났다. 절리의 간격은 클수록 변위가 국부적으로 집중되어 응답이 증가하는 것으로 나타났다. 절리의 경사는 정적인 경우와는 상반되게 $45^{\circ}$에서 응답이 가장 작게 계산되었다. 이는 전단변형으로 인하여 수직면과 수평면에 큰 전단응력이 작용하여 주응력면이 $45^{\circ}$로 회전하기 때문인 것으로 판단된다. 결과적으로 절리의 경사, 간격, 강성 모두 터널의 응답에 큰 영향을 미치며 특히 간격이 큰 미끈한 수직 또는 수평방향 절리가 터널을 통과하는 경우가 가장 위험한 것으로 나타났다. 또한 터널 라이닝에는 휨응력보다는 전단응력이 큰 폭으로 증가하므로 내진설계 시 이를 반드시 검토해야 할 것으로 판단된다.

• 대한기계학회논문집
• /
• 제17권1호
• /
• pp.90-100
• /
• 1993

본 연구에서는 압전 감지기/작동기를 이용한 복합적층판의 진동제어해석을 위 하여 판요소를 사용한 능률적인 유한요소코드 개발에 있다. 운동방정식은 고전 적층 판이론과 Hamilton의 법칙을 이용하여 유도하며 압전방정식으로부터 전기적-기계적 연 계를 고려한 감지식과 작동식을 구한다.각식들은 유한요소 보간함수에 의하여 절점 변위에 대한 행렬방정식으로 변환된다. 요소마다 하나의 전기적 자유도를 가진 4-절 점 12-자유도 판요소를 사용하여 효율적인 계산을 가능하게 하였다. 압전 감지기/작 동기를 도입함에 있어 하나의 전극에 대해 압전 감지기/작동기는 하나의 감지/작동전 압을 갖는다. 각 요소에 전극번호를 부가함으로써 다양한 형상의 전극을 쉽게 모델 링하였으며 전극의 특성도 충분히 고려하였다. 전기적 하중에 의한 압전보의 변형과 변형에 대한 감지전압에 대한 계산을 수행하여 기존의 연구와 비교함으로써 본 프로그 램의 타당성을 확인하였다. 나아가 여러가지 전극형상에 대한 복합재료 평판의 시간 영역과 주파수영역에서 응답을 계산하였다.

• 한국도로학회논문집
• /
• 제16권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
• /
• 제12권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.