• The Transactions of the Korean Institute of Electrical Engineers B
• /
• v.52 no.6
• /
• pp.247-253
• /
• 2003

A novel finite element mesh regeneration method is presented for 3D shape optimization of electromagnetic devices. The method has its theoretical basis in the structural deformation of an elastic body. When the shape of the electromagnetic devices changes during the optimization process, a proper 3D finite element mesh can be easily obtained using the method from the initial mesh. For real engineering problems, the method guarantees a smooth shape with proper mesh quality, and maintains the same mesh topology as the initial mesh. Application of the optimum design of an electromagnetic shielding plate shows the effectiveness of the presented method.

• Transactions of Materials Processing
• /
• v.9 no.2
• /
• pp.140-151
• /
• 2000

The finite element mesh approach for tool surface description is applied effectively to analyze sheet metal stamping processes. To improve the mesh quality and the stability of the mesh generation process, a gybrid method based on the grid approach and the Delaunay triangulation is proposed in the present work. In the present study, a general method for the mathematical description of arbitrarily shaped tool surface is proposed by introducing the parametric surface approach. A polynomial function employed to describe the base parametric surface and the boundary curves are defined to describe arbitrary three-dimensional trimmed surfaces. To verify the validity of the proposed method, automatic mesh generation is carried out for some shosen complicated parts including actual automotive panel.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.12
• /
• pp.3159-3174
• /
• 1994

A simple octree encoding algorithm based on a tetrahedron root has been developed to be used for fully automatic generation of three dimensional finite element meshes. This algorithm starts octree decomposition from a tetrahedron root node instead of a hexahedron root node so that the terminal mode has the same topology as the final tetrahedral mesh. As a result, the terminal octant can be used as a tetrahedral finite element without transforming its topology. In this part(I) of the thesis, an efficient algorithm for the tetrahedron-based octree is proposed. For this development, the following problems have been solved, : (1) an efficient data structure for storing the octree and finite elements, (2) an encoding scheme of a tetrahedral octree, (3) a neighbor finding technique for the tetrahedron-based octree.

• Structural Engineering and Mechanics
• /
• v.76 no.3
• /
• pp.379-393
• /
• 2020

In this paper, we propose an automatic procedure to improve the accuracy of finite element solutions using enriched 2D solid finite elements (4-node quadrilateral and 3-node triangular elements). The enriched elements can improve solution accuracy without mesh refinement by adding cover functions to the displacement interpolation of the standard elements. The enrichment scheme is more effective when used adaptively for areas with insufficient accuracy rather than the entire model. For given meshes, an error for each node is estimated, and then proper degrees of cover functions are applied to the selected nodes. A new error estimation method and cover function selection scheme are devised for the proposed adaptive enrichment scheme. Herein, we demonstrate the proposed enrichment scheme through several 2D problems.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2003.04a
• /
• pp.288-295
• /
• 2003

An adaptive procedure in finite element analysis is presented by p-refinement of meshes in conjunction with a posteriori error estimator that is based on the recovery technique. In case of the recovery technique, the SPR(superconvergent patch recovery) approach has been modified for p-adaptive mesh refinement. The strategy of finding a nearly optimal distribution of polynomial degrees on a fixed finite element mesh is discussed such that a particular element has to be refined automatically to obtain an acceptable level of accuracy by increasing p-levels non-uniformly. To verify the proposed algorithm, the limit value approach is proposed which utilizes the exact strain energy computed from the extrapolation equation. A new pre-processor is developed for the p-version finite element program in which the vector graphic editor is used for the automatic generation of node connection and coordinate by halfedge solid data structure according to uniform or nonuniform p-distribution. The general 2-D algorithm is also developed to generate face modes and internal modes in accordance with different mesh types. The quality of the error estimator is investigated with the help of two mumerical examples. The results show that the sequences of p-distributions obtained by the proposed error indicator closely follow the optimal trajectory.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.202-207
• /
• 2004

An elastic-plastic finite element analysis of fatigue crack closure is performed for plane strain conditions. The stabilization behavior of crack opening level and the effect of mesh size on the crack opening stress are investigated. In order to obtain a stabilized crack opening level for plane strain conditions, the crack must be advanced through approximately four times the initial monotonic plastic zone. The crack opening load tends to increase with the decrease of mesh size. The mesh size nearly equal to the theoretical plane strain cyclic plastic zone size may provide reasonable numerical results comparable with experimental crack opening data. The crack opening behavior is influenced by the crack growth increment and discontinuous opening behavior is observed. A procedure to predict the most appropriate mesh size for different stress ratio is suggested. Crack opening loads predicted by the FE analysis based on the procedure suggested resulted in good agreement with experimental ones within the error of 5 %. Effect of the distance behind the crack tip on the crack opening load determined by the ASTM compliance offset method based on the load-displacement relation and by the rotational offset method based on the load-differential displacement relation is investigated. Optimal gage location and method to determine the crack opening load is suggested.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.28 no.9
• /
• pp.1359-1367
• /
• 2004

The purpose of this study is to predict the dynamic transmission error of the helical gear system. To do so, the equations of motion in the helical gear system which consists of motor, coupling, gear, torque sensor, and brake are derived. As the input parameters, the mass moment of inertia by a 3D CAD software and the equivalent stiffness of the bearings and shaft are calculated and the coupling stiffness is measured. The static transmission error as an excitation is calculated by in-house program. Dynamic transmission error is predicted by solving the equations of motion. Mode shape, the dynamic mesh force and the bearing force are also calculated. In this analysis, the relationship between the dynamic mesh force and the bearing force and mode shape behavior in gear mesh are checked. As a result, the magnitude of mesh force is highly related with the gear mesh behavior in mode shape. The finite element analysis is conducted to find out the natural frequency of gear system. The natural frequencies by finite element analysis have a good agreement with the results by equation of motion. Finally, dynamic transmission error is measured by the specially designed experiment and the results by equation of motion are validated.

• Computational Structural Engineering
• /
• v.9 no.3
• /
• pp.145-155
• /
• 1996

This paper suggests a new algorithm of automatic mesh generation over planar domains with arbitrarily shaped boundaries and control curves. The algorithm is based on the method of recursively subdividing the domain by the path connecting, with minimum penalty value, two points on the super-loop, which consists of the boundaries and the control curves, The algorithm is not subject to any limitation on the shape of the domain, and its process can be fully automated. Therefore, this algorithm can be implemented into computer programs which require minimal user intervention while generating finite element meshes over complicated domains. This algorithm can also be easily extended for application to the generation of meshes over curved surfaces, or to the adaptive mesh generation.

• Journal of Advanced Marine Engineering and Technology
• /
• v.22 no.4
• /
• pp.490-497
• /
• 1998

An elast-biscoplastic finite element analysis is performed to investigate detailed growth behavior of creep cracks and the numerical results are compared with experimental results. The results of mesh translation method are compared with those of node release method. Load line displancement curve obtained from the crack growth analysis by mesh translation shows the improved results than that obtained from the crack growth by node release method when the secondary creep rate is only used as creep material property. The results of accounting for primary creep rate and instantaneous plasticity shows a good agreement with the experimental result.

• Wind and Structures
• /
• v.3 no.4
• /
• pp.291-302
• /
• 2000

A new finite element technique to solve the problem of wind and structure interactions is presented. Conventionally, wind analysis is performed on the Eulerian description in which the finite element mesh would not move in accordance with the wind flow. However, it is not the case in wind-structure interaction problems because nodes attached to the surface of structure should move with the displacement of structure. The arbitrary Lagrangian-Eulerian (ALE) method treats the mesh and flow independently, and allow the mesh to move. In this study, the analysis domain is divided into regions of the structure, air around the structure and the interface of two regions. To satisfy the compatibility and equilibrium conditions between separated regions and to carry out the efficient analysis, the rigid link is used. Also the equation of wind and that of structure are arranged in a single matrix equation.