• Proceedings of the Korea Concrete Institute Conference
• /
• 1996.04a
• /
• pp.247-252
• /
• 1996

In the concrete structures, cracks occur in various causes and the cracks seriously affect the functions of structures. The analysis techniques of progressive crack in the concrete have been improved with the advance of numerical techniques. The discrete crack model used in finite element program for the analysis of progressive failure is very effective, but it can not be easily implemented into numerical procedures because of difficult handing of nodal points in finite element meshes for crack growth. This paper introduces one of the techniques which skips the difficulty. In this paper, the modeling of progressive failure using finite element formulation is explained for the analysis of concrete fracture. The discontinuous element using the discontinuous shape function and the dual mapping technique in the numerical integration are implemented into finite element code for this purpose. It is shown that developed finite element program can predict the quasi-brittle behavior of concrete including ultimate load. The comparisons of the analysis results with other data are also shown.

• Computational Structural Engineering
• /
• v.5 no.2
• /
• pp.105-118
• /
• 1992

Static performance was compared for the triangular plate elements through some numerical experiments. Four Kirchhoff elements and six Mindlin elements were selected for the comparison. Numerical tests were executed for the problems of rectangular plates with regular and distorted meshes, rhombic plates, circular plates and cantilever plates. Among the Kirchhoff 9 DOF elements, the discrete Kirchhoff theory element was the best. Element distortion and the aspect ratio were shown to have negligible effects on the displacement behaviour. The Specht's element resulted in better results than the Bergan's but it was sensitive to the aspect ratio. The element based on the hybrid stress method also resulted in good results but it assumed to be less reliable. Among the linear Mindlin elements, the discrete shear triangle was the best in view of reliability, accuracy and convergence. Since the thin plate behaviour of it was as good as the DKT element, it can be used effectively in the finite element code regardless of the thickness. As a quadratic Mindlin element, the MITC7 element resulted in best results in almost all cases considered. The results were at least as good as those of doubly refined meshes of linear elements.

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

We present applications of MLS-based finite elements, which enable us to easily treat highly complex nonmatching finite element meshes and discontinuities. The shape functions of MLS-based finite element can be easily generated with the aid of Moving Least Square approximation on the parental domain. The major advantage includes that the position of element nodes as well as the number of the element nodes can be conveniently adjusted according to the nature of the problems under consideration, so that finite-element mesh is straightforwardly adapted to evolving discontinuities and. interfaces. Furthermore, we show that the present MLS-based finite elements are efficiently applied for elastic-plastic deformations, wherein the implicit constraint of incompressibility should be properly handled.

• Biomaterials and Biomechanics in Bioengineering
• /
• v.3 no.2
• /
• pp.71-84
• /
• 2016

The finite element method is wide used in simulation in the biomechanical structures, but a lack of studies concerning finite element mesh quality in biomechanics is a reality. The present study intends to analyze the importance of the mesh quality in the finite element model results from humeral structure. A sensitivity analysis of finite element models (FEM) is presented for the humeral bone and cartilage structures. The geometry of bone and cartilage was acquired from CT scan and geometry reconstructed. The study includes 54 models from same bone geometry, with different mesh densities, constructed with tetrahedral linear elements. A finite element simulation representing the glenohumeral-joint reaction force applied on the humerus during $90^{\circ}$ abduction, with external load as the critical condition. Results from the finite element models suggest a mesh with 1.5 mm, 0.8 mm and 0.6 mm as suitable mesh sizes for cortical bone, trabecular bone and humeral cartilage, respectively. Relatively to the higher minimum principal strains are located at the proximal humerus diaphysis, and its highest value is found at the trabecular bone neck. The present study indicates the minimum mesh size in the finite element analyses in humeral structure. The cortical and trabecular bone, as well as cartilage, may not be correctly represented by meshes of the same size. The strain results presented the critical regions during the $90^{\circ}$ abduction.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.5
• /
• pp.385-392
• /
• 2013

The finite element method has become the most widely used method of structural analysis and recently, the method has often been applied to complex dynamic and nonlinear structural analyses problems. Even for these complex problems, where the responses are hard to predict, finite element analyses yield reliable results if appropriate element types and meshes are used. However, the dynamic and nonlinear behaviors of a structure often include large deformations in various portions of the structure and if the same mesh is used throughout the analysis, some elements may deform to shapes beyond the reliable limits; thus dynamically adapting finite element meshes are needed in order for the finite element analyses to be accurate. In addition, to satisfy the users requirement of quick real run time of finite element programs, the algorithms must be computationally efficient. This paper presents an adaptive finite element mesh generation scheme for dynamic analyses of structures that may adapt at each time step. Representative strain values are used for error estimates and combinations of the h-method(node movement) and the r-method(element division) are used for mesh refinements. A coefficient that depends on the shape of an element is used to limit overly distorted elements. A simple frame example shows the accuracy and computational efficiency of the scheme. The aim of the study is to outline the adaptive scheme and to demonstrate the potential use in general finite element analyses of dynamic and nonlinear structural problems commonly encountered.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.700-705
• /
• 2004

Shell finite elements are widely used for the analysis of thin section objects such as sheet metal parts, automobile bodies and et al. due to their computational efficiency. Since many of input data for finite element analysis are given as solid models or triangulated surface models, one should extract midsurface information from these input data initially and then construct shell meshes on the extracted midsurfaces. In this paper, a method of generating shell elements on midsurfaces directly from input models have been proposed. In order to construct shell meshes, the input models should be triangulated on surfaces first, and then tetrahedral elements are generated by using an advancing front method, and finally mid shell surfaces are obtained from tetrahedral meshes. Some examples are given to demonstrate the efficiency of the proposed method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.31 no.7 s.262
• /
• pp.739-746
• /
• 2007

Recently, models in STL files are widely used in reverse engineering processes, CAD systems and analysis systems. However the models have poor geometric quality and include only triangles, so the models are not suitable for the finite element analysis. This paper presents a general method that generates finite element mesh from STL file models. Given triangular meshes, the method estimates triangles and makes clusters which consist of triangles. The clusters are merged by some geometric indices. After merging clusters, the method applies plane meshing algorithm, based on domain decomposition method, to each cluster and then the result plane mesh is projected into the original triangular set. Because the algorithm uses general methods to generate plane mesh, we can obtain both tri and quad meshes unlike previous researches. Some mechanical part models are used to show the validity of the proposed method.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.12
• /
• pp.2995-3006
• /
• 1993

An automatic mesh generation scheme has been developed for finite element analysis with two-dimensional, quadrilateral elements. The basic strategies of the method are to transform the analysis domain into loops with key nodes and the loops are recursively subdivided into subloops with the use of best split lines. Finally by using the basic loop operators, the meshes are completed. In this algorithm an eight-node loop operator is proposed, which is useful in the area where the change of element size is large and the splitting criteria for subdividing the loops have also been modified to the existing algorithms. Lines, arcs, and cubic spline curves are used to define the boundaries of analysis domain. Sample meshes for several geometries are presented to demonstrate the robustness of the algorithm.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1994.10a
• /
• pp.64-73
• /
• 1994

An updated Lagrangian finite element analysis with automatic remeshing scheme is applied to the three-dimensional hot extrusion through landless square dies. In the remeshing procedure, it is very difficult that the meshes are generated automatically with consideration of physical characteristics. In the presented study, the mesh generation is accomplished by modular concept. The generated meshes by modular concept have advantages, especially for three-dimensional problems, such as economized computational time and consideration of physical characteristic. In the problem, orifice shapes of square die are divided into two for the extrusion of solid sections. The orifice adaptive modules are developed for each type and the numerical examples are carried out for each type.

• Proceedings of the KIEE Conference
• /
• 1998.07a
• /
• pp.241-243
• /
• 1998

This paper describes a new method for the faster shape optimization of the electromagnetic devices. In a conventional iterative method of shape design optimization using design sensitivity based on a finite element method, meshes for a new shape of the model are generated and a discretized system equation is solved using the meshes in each iteration. They cause much design time. To save this time, a polynomial approximation of the finite element solution with respect to the geometric design parameters using Taylor expansion is constructed. This approximate state variable expressed explicitly in terms of design parameters is employed in a gradient-based optimization method. The proposed method is applied to the shape design of quadrupole magnet.