• Korean Journal of Computational Design and Engineering
• /
• v.18 no.6
• /
• pp.461-469
• /
• 2013

Computer simulation with FEM is very useful to analyze hypervelocity impact phenomena that are tremendously expensive or otherwise too impractical to analyze experimentally. Shock physics can be efficiently handled by mesh adaptation which allows finite element mesh to be locally optimized to resolve moving shock wave in explosion. In this paper, an adaptive meshing technique based upon quadtree data structure was applied to resolve ballistic impact phenomena. The technique can adaptively refine a mesh in the neighborhood of a shock and coarsen the mesh for the smooth flow behind the shock according to a criterion. The criterion for refinement and coarsening is based upon the standard deviation of the gradient of shock pressure on the associated field. Shock simulation starts with the rough mesh of the pressure field and mesh density is increased locally under the criterion at each time step. The results show that the mesh adaptation enables to minimize the global computation error of FEM and to increase storage and computational saving compared to the fixed resolution of the conventional static mesh approach.

• Korean Journal of Computational Design and Engineering
• /
• v.8 no.1
• /
• pp.10-18
• /
• 2003

FEM (Finite Element Method) is a fundamental numerical analysis technique in wide spread use in engineering application. As the solving time occupies small portion of entire FEM analysis time because of development of hardware, the relative lime to the whole analysis time to make mesh mod-els is growing. In particular, in the case of stiffeners such as features attached to plate in ship structure, the line constraints are imposed on mesh model together with other constraints such as holes. To auto-matically generate two dimensional quadrilateral mesh with the line constraints, an algorithm is pro-posed based on the constrained Delaunay triangulation and Q-Morph algorithm in which the line constraints are not considered. The performance of the proposed algorithm is evaluated. And some numerical results of our proposed algorithm ate presented.

• Proceedings of the KSME Conference
• /
• 2007.05a
• /
• pp.780-785
• /
• 2007

Accurate and fast haptic simulations of deformable objects are desired in many applications such as medical virtual reality. In haptic interactions with a coarse model, the number of nodes near the haptic interaction region is too few to generate detailed deformation. Thus, local refinement techniques need to be developed. Many approaches have employed purely geometric subdivision schemes, but they are not proper in describing the deformation behavior of deformable objects. This paper presents a continuum mechanics-based finite element adaptive method to perform haptic interaction with a deformable object. This method superimposes a local fine mesh upon a global coarse model, which consists of the entire deformable object. The local mesh and the global mesh are coupled by the s-version finite element method (s-FEM), which is generally used to enhance accurate solutions near the target points even more. The s-FEM can demonstrate a reliable deformation to users in real-time.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2006.04a
• /
• pp.373-380
• /
• 2006

Accurate and fast haptic simulations of deformable objects are desired in many applications such as medical virtual reality. In haptic interactions with a coarse model, the number of nodes near the haptic interaction region is too few to generate detailed deformation. Thus, local refinement techniques need to be developed. Many approaches have employed purely geometric subdivision schemes, but they are not proper in describing the deformation behavior of deformable objects. This paper presents a continuum mechanics-based finite element adaptive method to perform haptic interaction 'with a deformable object. This method superimposes a local fine mesh upon a global coarse model, which consists of the entire deformable object. The local mesh and the global mesh are coupled by the s-version finite element method (s-FEM), which is generally used to enhance accurate solutions near the target points even more. The s-FEM can demonstrate a reliable deformation to users in real-time.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.11
• /
• pp.1950-1956
• /
• 1997

An engineering goal in vibration and noise professionals is to develope quiet machines at the preliminary design stage, and various numerical techniques such as FEM, SEA or BEM are one of the schemes toward the goal. In this paper, the research has been focused on the sensitivity effect of mesh sizes for FEM application so that the optimum size of the mesh that leads to engineering solution within acceptable computing time could be generated. In order to evaluate the mesh size effect, three important parameters have been examined : natural frequencies, number of modes and driving point mobility. First, several lower modes including the fundamental frequency of a 2-D plate structure have been calculated as mesh size changes. Since theoretical values of natural frequencies for a simple structure are known, the deviation between the numerical and theoretical values is obtained as a function of mesh size. The result shows that the error is no longer decreased if the mesh size becomes a quarter wavelength or smaller than that. Second, the mesh size effect is also investigated for the number of modes. For the frequency band up to 1.4 kHz, the structure should have 38 modes in total. As the mesh size reaches to the quarter wavelength, the total count in modes approaches to the same values. Third, a mobility function at the driving point is compared between SEA and FEM result. In SEA application, the mobility function is determined by the modal density and the mass of the structure. It is independent of excitation frequencies. When the mobility function is calculated from a wavelength to one-tenth of it, the mobility becomes constant if the mesh becomes a quarter wavelength or smaller. We can conclude that dynamic parameters, such as eigenvalues, mode count, and mobility function, can be correctly estimated, while saving the computing burden, if a quarter wavelength (.lambda./4) mesh is used. Therefore, (.lambda./4) mesh is recommended in structural vibration analysis.

• Proceedings of the KIEE Conference
• /
• 1999.07e
• /
• pp.2044-2046
• /
• 1999

In this paper, the propagation of corona streamer was simulated using finite element method(FEM) combined with Flux-Corrected Transport(FCT) algorithm. To obtain more effective grid distribution, the adaptive mesh generation scheme was also employed. The employed adaptive scheme can refine the mesh where needed. In addition, it is possible to coarsen the unnecessarily dense mesh. Two models were analyzed with proposed method. And the results are shown.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.1870-1875
• /
• 2003

Optimization has been used in many engineering problems and must be one of the essential components during design process now. It is the process of maximizing the performance called an objective function of a system while satisfying some constraints, so finite element method is generally required in order to obtain these values during optimization. However, there are some difficulties to obtain them by means of FEM, where the changes of design variables cause the distortion and the regeneration of mesh that may result in inaccuracy and inefficiency. In order to overcome this problem, this paper proposed an alternative that the eXtended FEM introduced and developed by Ted Belytschko was applied to the optimization process because the key points of the X-FEM lie in that the discontinuity can be represented independently on the mesh by a function called in an enrichment function.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.21 no.2
• /
• pp.283-294
• /
• 2023

APro, developed in KAERI for the process-based total system performance assessment (TSPA) of deep geological disposal systems, performs finite element method (FEM)-based multiphysics analysis. In the FEM-based analysis, the mesh element quality influences the numerical solution accuracy, memory requirement, and computation time. Therefore, an appropriate mesh structure should be constructed before the mesh stability analysis to achieve an accurate and efficient process-based TSPA. A generic reference case of DECOVALEX-2023 Task F, which has been proposed for simulating stationary groundwater flow and time-dependent conservative transport of two tracers, was used in this study for mesh stability analysis. The relative differences in tracer concentration varying mesh structures were determined by comparing with the results for the finest mesh structure. For calculation efficiency, the memory requirements and computation time were compared. Based on the mesh stability analysis, an approach based on adaptive mesh refinement was developed to resolve the error in the early stage of the simulation time-period. It was observed that the relative difference in the tracer concentration significantly decreased with high calculation efficiency.

• The Mathematical Education
• /
• v.31 no.2
• /
• pp.121-132
• /
• 1992

The pointwise convergence define the relation-ship between the mesh-size and the tolerance. This will play an important role in improving quality of finite element approximate solution. In this paper, We evaluate the convergence on a certaon unknown point with a 1-irregular mesh refinement. This m that the degree of freedom is minimized within a tolerance.

• 제어로봇시스템학회:학술대회논문집
• /
• 2005.06a
• /
• pp.2410-2413
• /
• 2005

Finite element method (FEM) has been widely used as a useful numerical method that can analyze complex engineering problems in electro-magnetics, mechanics, and others. CEMTool, which is similar to MATLAB, is a command style design and analyzing package for scientific and technological algorithm and a matrix based computation language. In this paper, we present new 3D FEM package in CEMTool environment. In contrast to the existing CEMTool 2D FEM package and MATLAB PDE (Partial Differential Equation) Toolbox, our proposed 3D FEM package can deal with complex 3D models, not a cross-section of 3D models. In the pre-processor of 3D FEM package, a new 3D mesh generating algorithm can make information on 3D Delaunay tetrahedral mesh elements for analyses of 3D FEM problems. The solver of the 3D FEM package offers three methods for solving the linear algebraic matrix equation, i.e., Gauss-Jordan elimination solver, Band solver, and Skyline solver. The post-processor visualizes the results for 3D FEM problems such as the deformed position and the stress. Consequently, with our new 3D FEM toolbox, we can analyze more diverse engineering problems which the existing CEMTool 2D FEM package or MATLAB PDE Toolbox can not solve.