• 한국전산구조공학회논문집
• /
• 제20권3호
• /
• pp.321-327
• /
• 2007

본 연구는 탄성균열문제를 신속하고 정확하게 해석할 수 있는 새로운 개념의 그리드(grid) 없는 유한차분법을 제시한다. 이동최소제곱법을 이용한 Taylor 전개식 구성을 통해 직접적인 미분계산 없이 근사함수와 그 미분을 손쉽게 계산한다. 그리드로 인한 절점 간의 종속성이 없어 해석영역 내의 불연속면 모델링이 용이하여 차분식 구성시 균열로 인한 불연속 효과를 고려하는 과정도 자연스럽다. 유한차분법에 근간을 두고 있어 지배 미분방정식을 직접 이산화하기 때문에 수치적분이 필요한 수치기법에 비해 계산속도도 빠르다. 모드 I과 모드 II 균열문제 해석을 통해 본 해석기법이 정확하고 효율적으로 응력확대계수를 계산할 수 있음을 보였다.

• Wind and Structures
• /
• 제26권1호
• /
• pp.35-43
• /
• 2018

This paper describes the effect of different testing parameters (configuration of infrastructure and truck position on road) on truck aerodynamic coefficients under cross wind conditions, by means of a numerical approach known as Large Eddy Simulation (LES). In order to estimate the air flow behaviour around both the infrastructure and the truck, the filtered continuity and momentum equations along with the Smagorinsky-Lilly model were solved. A solution for these non-linear equations was approached through the finite volume method (FVM) and using temporal and spatial discretization schemes. As for the results, the aerodynamic coefficients acting on the truck model exhibited nearly constant values regardless of the Reynolds number. The flat ground is the infrastructure where the rollover coefficient acting on the truck model showed lowest values under cross wind conditions (yaw angle of $90^{\circ}$), while the worst infrastructure studied for vehicle stability was an embankment with downward-slope on the leeward side. The position of the truck on the road and the value of embankment slope angle that minimizes the rollover coefficient were determined by successfully applying the Response Surface Methodology.

• Journal of Mechanical Science and Technology
• /
• 제17권12호
• /
• pp.1844-1854
• /
• 2003

This paper discusses details of modeling and robust control of an AMB (active magnetic bearing) spindle, and part I presents a modeling and validation process of the AMB spindle. There are many components in AMB spindle : electromagnetic actuator, sensor, rotor, power amplifier and digital controller. If each component is carefully modeled and evaluated, the components have tight structured uncertainty bounds and achievable performance of the system increases. However, since some unknown dynamics may exist and the augmented plant could show some discrepancy with the real plant, the validation of the augmented plant is needed through measuring overall frequency responses of the actual plant. In addition, it is necessary to combine several components and identify them with a reduced order model. First, all components of the AMB spindle are carefully modeled and identified based on experimental data, which also render valuable information in quantifying structured uncertainties. Since sensors, power amplifiers and discretization dynamics can be considered as time delay components, such dynamics are combined and identified with a reduced order. Then, frequency responses of the open-loop plant are measured through closed-loop experiments to validate the augmented plant. The whole modeling process gives an accurate nominal model of a low order for the robust control design.

• 한국기계가공학회지
• /
• 제21권1호
• /
• pp.102-110
• /
• 2022

This paper proposes a solution to the problems of manufacturing cost and processability by applying discrete material and thickness optimization (DMTO) and minimizing the use of high-strength, lightweight materials in the optimization process. A simple infant pop-up seat model was selected as the application target, and the weight reduction effect and variation in strength according to the optimization results were observed. In this study, a simplified finite element model of an infant pop-up seat frame was first constructed. The model was used to perform a static structural analysis to verify the weight and strength of each part. The D-optimal design of the experimental method was then used to observe the influence of each part on the weight and strength. This process was applied using discrete thickness optimization (DTO) (which applies high-strength, lightweight materials and optimizes only the thickness) and DMTO (which considers both the material and thickness). The DTO and DMTO results were compared to verify the design method that determines the major parts and simultaneously considers the material and thickness. Accordingly, in this study, an optimal lightweight design that satisfied the strength standards of the seat frame was derived. Furthermore, discretization parameters were used to minimize the application of high-strength, lightweight materials.

• Nuclear Engineering and Technology
• /
• 제52권8호
• /
• pp.1626-1637
• /
• 2020

This work covers an important point of the benchmark released by the expert group on Uncertainty Analysis in Modeling of Light Water Reactors. This ambitious benchmark aims to determine the uncertainty in light water reactors systems and processes in all stages of calculation, with emphasis on multi-physics (coupled) and multi-scale simulations. The Gesellschaft für Anlagen und Reaktorsicherheit methodology is used to propagate the thermal-hydraulic uncertainty of macroscopic parameters through TRACE5.0p3/PARCSv3.0 coupled code. The main innovative points achieved in this work are i) a new thermal-hydraulic model is developed with a highly-accurate 3D core discretization plus an iterative process is presented to adjust the 3D bypass flow, ii) a control rod insertion occurrence -which data is obtained from a real PWR test- is used as a transient simulation, iii) two approaches are used for the propagation process: maximum response where the uncertainty and sensitivity analysis is performed for the maximum absolute response and index dependent where the uncertainty and sensitivity analysis is performed at each time step, and iv) RESTING MATLAB code is developed to automate the model generation process and, then, propagate the thermal-hydraulic uncertainty. The input uncertainty information is found in related literature or, if not found, defined based on expert judgment. This paper, first, presents the Gesellschaft für Anlagen und Reaktorsicherheit methodology to propagate the uncertainty in thermal-hydraulic macroscopic parameters and, then, shows the results when the methodology is applied to a PWR reactor.

• Interaction and multiscale mechanics
• /
• 제2권1호
• /
• pp.45-68
• /
• 2009

An aim of the study is to develop an efficient numerical simulation technique that can handle the two-scale analysis of fluid permeation filters fabricated by the partial sintering technique of small spherical ceramics. A solid-fluid mixture homogenization method is introduced to predict the mechanical characters such as rigidity and permeability of the porous ceramic filters from the micro-scale geometry and configuration of partially-sintered particles. An extended finite element (X-FE) discretization technique based on the enriched interpolations of respective characteristic functions at fluid-solid interfaces is proposed for the non-interface-fitted mesh solution of the micro-scale analysis that needs non-slip condition at the interface between solid and fluid phases of the unit cell. The homogenization and localization performances of the proposed method are shown in a typical two-dimensional benchmark problem whose model has a hole in center. Three-dimensional applications to the body-centered cubic (BCC) and face-centered cubic (FCC) unit cell models are also shown in the paper. The 3D application is prepared toward the computer-aided optimal design of ceramic filters. The accuracy and stability of the X-FEM based method are comparable to those of the standard interface-fitted FEM, and are superior to those of the voxel type FEM that is often used in such complex micro geometry cases.

• 수산해양기술연구
• /
• 제35권1호
• /
• pp.54-59
• /
• 1999

High-speed electronic digital computers have enabled engineers to employ various numerical discretization techniques for solutions of complex problems. The Finite Element Method is one of the such technique. The Finite Element Method is one of the numerical analysis based on the concepts of fundamental mathematical approximation. Three dimensional plate structures used often in partition of ship, box girder and frame are analyzed by Finite Element Method. In design of structures, the static deflections, stress concentrations and dynamic deflections must be considered. However, these problem belong to geometrically nonlinear mechanical structure analysis. The analysis of each element is independent, but coupling occurs in assembly process of elements. So, to overcome such a difficulty the shell theory which includes transformation matrix and a fictitious rotational stiffness is taken into account. Also, the Mindlin's theory which is considered the effect of shear deformation is used. The Mindlin's theory is based on assumption that the normal to the midsurface before deformation is "not necessarily normal to the midsurface after deformation", and is more powerful than Kirchoff's theory in thick plate analysis. To ensure that a small number of element can represent a relatively complex form of the type which is liable to occur in real, rather than in academic problem, eight-node quadratic isoparametric elements are used. are used.

• Computers and Concrete
• /
• 제5권4호
• /
• pp.343-358
• /
• 2008

This article focuses on concrete structures submitted to impact loading and is aimed at predicting local damage in the vicinity of an impact zone as well as the global response of the structure. The Discrete Element Method (DEM) seems particularly well suited in this context for modeling fractures. An identification process of DEM material parameters from macroscopic data (Young's modulus, compressive and tensile strength, fracture energy, etc.) will first be presented for the purpose of enhancing reproducibility and reliability of the simulation results with DE samples of various sizes. The modeling of a large structure by means of DEM may lead to prohibitive computation times. A refined discretization becomes required in the vicinity of the impact, while the structure may be modeled using a coarse FE mesh further from the impact area, where the material behaves elastically. A coupled discrete-finite element approach is thus proposed: the impact zone is modeled by means of DE and elastic FE are used on the rest of the structure. The proposed approach is then applied to a rock impact on a concrete slab in order to validate the coupled method and compare computation times.

• Structural Engineering and Mechanics
• /
• 제61권2호
• /
• pp.193-207
• /
• 2017

In this study, the vibration of an electrostatically actuated micro cantilever beam is analyzed in which a viscoelastic layer covers a portion of the micro beam length. This proposed model is considered as the main element of mass and pollutant micro sensors. The nonlinear motion equation is extracted by means of Hamilton principle, considering nonlinear shortening effect for Euler-Bernoulli beam. The non-linear effects of electrostatic excitation, geometry and inertia have been taken into account. The viscoelastic model is assumed as Kelvin-Voigt model. The motion equation is discretized by Galerkin approach. The linear free vibration mode shapes of non-uniform micro beam i.e. the linear mode shape of the system by considering the geometric and inertia effects of viscoelastic layer, have been employed as comparison function in the process of the motion equation discretization. The discretized equation of motion is solved by the use of multiple scale method of perturbation theory and the results are compared with the results of numerical Runge-Kutta approach. The frequency response variations for different lengths and thicknesses of the viscoelastic layer have been founded. The results indicate that if a constant volume of viscoelastic layer is to be deposited on the micro beam for mass or gas sensor applications, then a modified configuration may be found by using the analysis of this paper.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 2004년도 추계학술대회논문집
• /
• pp.75-78
• /
• 2004

A mesh generation algorithm adapted to the mesh density map using the Delaunay mesh generation technique is developed. In the finite element analyses of the forging processes, the numerical error increases as the process goes on because of discrete property of the finite elements or severe distortion of elements. Especially, in the region where stresses and strains are concentrated, the numerical discretization error will be highly increased. However, it is too time consuming to use a uniformly fine mesh in the whole domain to reduce the expected numerical error. Therefore, it is necessary to construct locally refined mesh at the region where the error is concentrated such as at the die corner. In this study, the point insertion algorithm is used and the mesh size is controlled by moving nodes to optimized positions according to a mesh density map constructed with a posteriori error estimation. An optimization technique is adopted to obtain a good position of nodes. And optimized smoothing techniques are also adopted to have smooth distribution of the mesh and improve the mesh element quality.