• 제어로봇시스템학회:학술대회논문집
• /
• 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.

• International Journal of Fuzzy Logic and Intelligent Systems
• /
• v.5 no.2
• /
• pp.119-123
• /
• 2005

This paper describes a fuzzy-based system for analyzing the stress intensity factors (SIFs) of three-dimensional (3D) cracks. 3D finite element method(FEM) was used to obtain the SIF for subsurface cracks and surface cracks existing in inhomogeneous materials. A geometry model, i.e. a solid containing one or several 3D cracks is defined. Several distributions of local node density are chosen, and then automatically superposed on one another over the geometry model by using the fuzzy theory. Nodes are generated by the bucketing method, and ten-noded quadratic tetrahedral solid elements are generated by the Delaunay triangulation techniques. The singular elements such that the mid-point nodes near crack front are shifted at the quarter-points, and these are automatically placed along the 3D crack front. The complete FE model is generated, and a stress analysis is performed. The SIFs are calculated using the displacement extrapolation method. The results were compared with those surface cracks in homogeneous materials. Also, this system is applied to analyze cladding effect of surface cracks in inhomogeneous materials.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.57 no.5
• /
• pp.43-49
• /
• 2015

Most of agricultural structures located in seashore could not avoid rapid deterioration of concrete because chloride-ion and $CO_2$ gradually penetrate into concrete. However, since most of models can be able to describe the phenomenon of penetration by using one or two dimensional models based on finite difference method (FDM), those modes can not simulate the real geometry and it takes a lot of computational time to complete even the calculation. To overcome those weaknesses, three dimensional numerical model considering time dependent variables such as surface concentration of chloride and diffusion coefficient of domain based on finite element method (FEM) was suggested. This model also included the neutralization occurred by the penetration of $CO_2$. Because the model used various sizes of tetrahedral mesh instead of equivalent rectangular mesh, it reduced the computational time to compare with FDM. As this model is based on FEM, it will be easily extended to execute multi-physics simulation including water evaporation and temperature change of concrete.

• Tribology and Lubricants
• /
• v.30 no.4
• /
• pp.205-211
• /
• 2014

Recent studies emphasize the importance of pivot stiffness in the analysis of tilting pad bearings (TPBs). The present paper develops a finite element model of the pad pivot and compares the predicted pivot stiffness to the results of Hertzian contact model calculations. Specifically, a finite element analysis generates tetrahedral mesh models with ~40,000 nodes for a ball-socket pivot and ~50,000 nodes for a rocker-back pivot. These models assume a frictionless boundary condition in the contact area. Increasing the applied loads on the pad in conjunction with increasing time steps ensures rapid convergence during the nonlinear numerical analysis. Predictions are performed using the developed finite element model for increasing the differential diameters between the pad pivot (or ball) and the bearing housing (or socket). The predictions show that the pivot contact area increases with decreasing differential diameters and increasing applied loads. Further, the maximum deformation occurring at the pivot center increases with increasing differential diameters and increasing applied loads. The pivot stiffness increases nonlinearly with decreasing differential diameters and increasing applied loads. Comparisons of results of the developed finite element model to those of Hertzian contact model calculations assuming a small contact area show that the latter model underestimates the pivot stiffnesses predicted by the finite element models of the ball-socket and rocker-back pivots, particularly for small differential diameters. This result implies the need for cautionduring the design of pivot stiffness by the Hertzian contact model.

• Journal of the Korean Society for Precision Engineering
• /
• v.25 no.2
• /
• pp.148-155
• /
• 2008

The aim of the study is to interpret the distribution of occlusal force by 3-dimensional finite element analysis of ISP(Implant Supported Prosthesis) supported by minimum number of implant to restore the edentulous patients. For this study, the Astra Tech implant system is used. Geometric modeling for 6 and 4 fixture ISP group is performed with respect to the bone, implant and one piece superstructure, respectively. Implants are arbitrarily placed according to the anatomical limit of lower jaw and for the favorable distribution of occlusal force, which is applied at the end of cantilever extension of ISP with 30mm. Element type is tetrahedral for finite element model and the typical mechanical properties, Young's modulus and Poisson's ratio of each material, cortical, cancellous bone and implant material are utilized for the finite element analysis. From this study, we can see the distribution of equivalent stress equal to real situation and speculate the difference in the stress distribution in the whole model and at each implant fixture, From the analysis, the area of maximum stress is distributed on distal contact area between bone and fixture in the crestal bone. The maximum stress is 53MPa at the 0.2mm area from the bone-implant interface in the maximum side for 300N load condition for 4 fixture case, which is slightly less than the stress calculated from allowable strain. This stress has not been deduced to directly cause the loss of crestal bone around implant fixture, but the stress can be much reduced as the old peoples may have lower chewing force. Thus, clinical trial may be performed with this treatment protocol to use 4 fixtured ISP for old patients.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2006.05a
• /
• pp.660-665
• /
• 2006

This paper presents a finite element method to analyze the free and forced vibration of a hard disk drive (HDD) considering the flexibility of a spinning disk-spindle with fluid dynamic bearings (FDBs), an actuator with pivot bearings, an air bearing between head-disk interface and the base with complicated geometry. Finite element equation of each component is consistently derived with the satisfaction of the geometric compatibility of the internal boundary between each component. The spinning disk, hub and FDBs are modeled by annular sector elements, beam elements and stiffness and damping elements, respectively. The actuator am, E-block, suspension and base plate are modeled by tetrahedral elements. The pivot bearing in the actuator and the air bearing between head-disk interfaces are modeled by the stiffness element with five degrees of freedom and the axial stiffness, respectively. A global matrix equation obtained by assembling the finite element equations of each substructure is transformed to a state-space matrix-vector equation, and both damped natural frequencies and modal damping ratios are calculated by solving the associated eigenvalue problem with the restarted Arnoldi iteration method. Modal and shock testing are performed to show that the proposed method well predicts the vibration characteristics of a HDD.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1997.04a
• /
• pp.496-500
• /
• 1997

This paper describes an automated system for analyzing the stress intensity factors(SIFs) of three-dimensional (3D) cracks. A geometry model, i.e.a solid containing one or several 3D cracks is defined. Several distributions of local node density are chosen, and then automatically superposed on one another over the geometry model by using the fuzzy knowledge processing. Nodes are generated by the bucketing method, and ten-noded quadratic tetrahedral solid elements are generated by the Delauuay triangulation techniques. The singular elements such that the mid-point nodes near crack fornt are shifted at the quarter-points are automatically placed along the 3D crack front. THe complete finite element (FE) model generated, i.e the mesh with material properties and boundary conditions is given to one of the commercial FE codes, and a stress analysis is performed. The SIFs are calculated using the displacement extrapolation method. To demonstrate practical performance of the present system, a semi- elliptical surface crack in a plate subjected to tension is solved.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.3
• /
• pp.139-170
• /
• 2013

We extend a p-hierarchical decomposition of the second degree finite element space of N$\acute{e}$d$\acute{e}$lec for tetrahedral meshes in three dimensions given in [1] to meshes with hexahedral elements, and derive p-hierarchical decompositions of the second degree finite element space of Raviart-Thomas in two dimensions for triangular and quadrilateral meshes. After having proved stability of these subspace decompositions and requiring certain saturation assumptions to hold, we construct a local a posteriori error estimator for fem and bem coupling of a time-harmonic electromagnetic eddy current problem in $\mathbb{R}^3$. We perform some numerical tests to underline reliability and efficiency of the estimator and test its usefulness in an adaptive refinement scheme.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.13 no.8
• /
• pp.3295-3300
• /
• 2012

Three Dimensional finite element method (FEM) was used to obtain the stress intensity factor for subsurface cracks and surface cracks existing in inhomogeneous materials. A geometry model, i.e. a solid containing one or several 3D cracks is defined. Several distributions of local node density are chosen, and then automatically superposed on one another over the geometry model. Nodes are generated by the bubble packing, and ten-noded quadratic tetrahedral solid elements are generated by the Delaunay triangulation techniques. To examine accuracy and efficiency of the present system, the stress intensity factor for a semi-elliptical surface crack in a plate subjected to uniform tension is calculated, and compared with Raju-Newman's solutions. Then the system is applied to analyze interaction effects of two dissimilar semi-elliptical cracks in a plate subjected to uniform tension.

• Journal of Computational Design and Engineering
• /
• v.1 no.2
• /
• pp.88-95
• /
• 2014

Because the cost of performance testing using actual products is expensive, manufacturers use lower-cost computer-aided design simulations for this function. In this paper, we propose using hexahedral meshes, which are more accurate than tetrahedral meshes, for finite element analysis. We propose automatic hexahedral mesh generation with sharp features to precisely represent the corresponding features of a target shape. Our hexahedral mesh is generated using a voxel-based algorithm. In our previous works, we fit the surface of the voxels to the target surface using Laplacian energy minimization. We used normal vectors in the fitting to preserve sharp features. However, this method could not represent concave sharp features precisely. In this proposal, we improve our previous Laplacian energy minimization by adding a term that depends on multi-normal vectors instead of using normal vectors. Furthermore, we accentuate a convex/concave surface subset to represent concave sharp features.