• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.1
• /
• pp.85-94
• /
• 2002

In this study an adaptive node generation procedure in the Element-free Galerkin (EFG) method using bubble-meshing technique is Proposed. Since we construct the initial configuration of nodes by subdivision of background cell, abrupt changes of inter-nodal distance between higher and lower error regions are unavoidable. This unpreferable nodal spacing induces additional errors. To obtain the smooth nodal configuration the nodal configurations are regenerated by bubble-meshing technique. This bubble meshing technique was originally developed to generate a set of well-shaped triangles and tetrahedra. In odder to evaluate the effect of abrupt changes of nodal spacing, one-dimensional problems with various nodal configurations mere investigated. To demonstrate the performance of proposed scheme, the sequences of making optimal nodal configuration with bubble meshing technique are investigated for several problems.

• Journal of Ocean Engineering and Technology
• /
• v.13 no.1 s.31
• /
• pp.39-50
• /
• 1999

특이 가중함수로 표현된 shepard interpolant와 일관조건을 사용하여 무요소법 형성함수를 도출하였다. 따라서 통상의 EFGM(Element Free Galerkin Method)과는 달리 변위로 주어지는 경계조건을 자연스럽게 부과할 수 있다. 수치계산 예로서 외팔보 문제를 다루었는데 보이론과 비교하여 매우 잘 맞는 결과를 보여주고, 유한요소법과의 결합도 자연스럽게 이루어짐을 보인다. 또 penny-shaped 균열을 다루는데, 응력확대계수는 균열 표면의 변위로부처 직접 계산하여 해석해와 비교한다.

• Steel and Composite Structures
• /
• v.22 no.2
• /
• pp.277-299
• /
• 2016

In this paper, free vibration, forced vibration, resonance and stress wave propagation behavior in nanocomposite plates reinforced by wavy carbon nanotube (CNT) are studied by a mesh-free method based on first order shear deformation theory (FSDT). The plates are resting on Winkler-Pasternak elastic foundation and subjected to periodic or impact loading. The distributions of CNTs are considered functionally graded (FG) or uniform along the thickness and their mechanical properties are estimated by an extended rule of mixture. In the mesh-free analysis, moving least squares (MLS) shape functions are used for approximation of displacement field in the weak form of motion equation and the transformation method is used for imposition of essential boundary conditions. Effects of CNT distribution, volume fraction, aspect ratio and waviness, and also effects of elastic foundation coefficients, plate thickness and time depended loading are examined on the vibrational and stresses wave propagation responses of the nanocomposite plates reinforced by wavy CNT.

• Transactions of Materials Processing
• /
• v.13 no.7
• /
• pp.614-622
• /
• 2004

A new approach for the automatic elimination of free edges in the finite element model for the analysis of sheet metal forming processes is presented. In general, the raw finite element model constructed from an automatic mesh generator is not well suited for the direct use in the downstream forming analysis due to the many free edges which requires tedious time consuming interactive graphic operations of the users. In the present study, a general method for the automatic elimination of free edges is proposed by introducing a CAD/CAE hybrid method. In the method a trimmed parametric surface is generated to fill the holes which are orginated from the free edges by using the one step elastic finite element analysis. In addition, mesh generation algorithm is suggested which can be used in the general trimmed surface. In order to verify the validity of the proposed method, various examples including actual automobile sheet metal parts are given and discussed.

• Transactions of Materials Processing
• /
• v.12 no.4
• /
• pp.401-407
• /
• 2003

This paper newly proposes a mesh regularization method for the enhancement of the efficiency in sheet metal forming analysis. The regularization method searches for distorted elements with appropriate searching criteria and constructs patches including the elements to be modified. Each patch is then extended to a three-dimensional surface in order to obtain the information of the continuous coordinates. In constructing the surface enclosing each patch, NURBS(Non-Uniform Rational B-Spline) surface is employed to describe a three-dimensional free surface. On the basis of the constructed surface, each node is properly arranged to form unit elements as close as to a square. The state variables calculated from its original mesh geometry are mapped into the new mesh geometry for the next stage or incremental step of a forming analysis. The analysis results with the proposed method are compared to the results from the direct forming analysis without mesh regularization in order to confirm the validity of the method.

• Computers and Concrete
• /
• v.16 no.6
• /
• pp.933-961
• /
• 2015

Damage by high-speed impact fracture is a dominant mode of failure in several applications of concrete structures. Numerical modelling can play a crucial role in understanding and predicting complex fracture processes. The commonly used mesh-based Finite Element Method has difficulties in accurately modelling the high deformation and disintegration associated with fracture, as this often distorts the mesh. Even with careful re-meshing FEM often fails to handle extreme deformations and results in poor accuracy. Moreover, simulating the mechanism of fragmentation requires detachment of elements along their boundaries, and this needs a fine mesh to allow the natural propagation of damage/cracks. Smoothed Particle Hydrodynamics (SPH) is an alternative particle based (mesh-less) Lagrangian method that is particularly suitable for analysing fracture because of its capability to model large deformation and to track free surfaces generated due to fracturing. Here we demonstrate the capabilities of SPH for predicting brittle fracture by studying a slender concrete structure (column) under the impact of a high-speed projectile. To explore the effect of the projectile material behaviour on the fracture process, the projectile is assumed to be either perfectly-elastic or elastoplastic in two separate cases. The transient stress field and the resulting evolution of damage under impact are investigated. The nature of the collision and the constitutive behaviour are found to considerably affect the fracture process for the structure including the crack propagation rates, and the size and motion of the fragments. The progress of fracture is tracked by measuring the average damage level of the structure and the extent of energy dissipation, which depend strongly on the type of collision. The effect of fracture property (failure strain) of the concrete due to its various compositions is found to have a profound effect on the damage and fragmentation pattern of the structure.

• Korean Journal of Computational Design and Engineering
• /
• v.12 no.3
• /
• pp.220-232
• /
• 2007

Current geometric modeling and analysis are commonly based on B-Rep modeling and a finite elements method respectively. Furthermore, it is difficult to represent an object whose material property is heterogeneous using the B-Rep method because the B-Rep is basically used for homogeneous models. In addition, meshes are required to analyze a property of a model when the finite elements method is applied. However, the process of generating meshes from B-Rep is cumbersome and sometimes difficult especially when the model is deformed as time goes by because the topology of deforming meshes are changed. To overcome those problems in modeling and analysis including homogeneous and heterogeneous materials, we suggest a unified modeling and analysis method based on implicit representation of the model using R-function which is suggested by Rvachev. For implicit modeling of an object a distance field is approximated and blended for a complex object. Using the implicit function mesh-free analysis is possible where meshes are not necessary. Generally mesh-free analysis requires heavy computational cost compared to a finite elements method. To improve the computing time of function evaluation, we utilize GPU programming. Finally, we give an example of a simple pipe design problem and show modeling and analysis process using our unified modeling and analysis method.

• Ocean Systems Engineering
• /
• v.1 no.4
• /
• pp.317-336
• /
• 2011

A time domain finite element based method is employed to analyze wave radiation by multiple cylinders. The nonlinear free surface and body surface boundary conditions are satisfied based on the perturbation method up to the second order. The first- and second-order velocity potential problems at each time step are solved through a finite element method (FEM). The matrix equation of the FEM is solved through an iteration and the initial solution is obtained from the result at the previous time step. The three-dimensional (3D) mesh required is generated based on a two-dimensional (2D) hybrid mesh on a horizontal plane and its extension in the vertical direction. The hybrid mesh is generated by combining an unstructured grid away from cylinders and two structured grids near the cylinder and the artificial boundary, respectively. The fluid velocity on the free surface and the cylinder surface are calculated by using a differential method. Results for various configurations including two-cylinder and four-cylinder cases are provided to show the mutual influence due to cylinders on the first and second waves and forces.

• Structural Engineering and Mechanics
• /
• v.50 no.2
• /
• pp.181-199
• /
• 2014

Locking phenomenon is a mesh problem and can be staved off with mesh refinement. If the studier is not preferred going to the solution with increasing mesh size or the computer memory can stack over flow than using higher order plate finite element or using integration techniques is a solution for this problem. The purpose of this paper is to show the shear locking phenomenon can be avoided by increase low order finite element mesh size of the plates and to study shear locking-free analysis of thick plates using Mindlin's theory by using higher order displacement shape function and to determine the effects of various parameters such as the thickness/span ratio, mesh size on the linear responses of thick plates subjected to uniformly distributed loads. A computer program using finite element method is coded in C++ to analyze the plates clamped or simply supported along all four edges. In the analysis, 4-, 8- and 17-noded quadrilateral finite elements are used. It is concluded that 17-noded finite element converges to exact results much faster than 8-noded finite element, and that it is better to use 17-noded finite element for shear-locking free analysis of plates.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2001.04a
• /
• pp.335-342
• /
• 2001

In this paper, an improved Element-Free Galerkin (EFG) method is proposed by adding enrichment function to the standard EFG approximation and a discontinuity function is implemented in constructing the shape function across the crack surface. In this method, the singularity and the discontinuity of the crack are efficiently modeled by using initial node distribution to evaluate reliable stress intensity factor, though the standard EFG method requires placing additional nodes near the crack tip. The proposed method enables the initial node distribution to be kept without any additional nodal d.o.f. and expresses the asymptotic stress field near the crack tip successfully. Numerical example verifies the improvement and the effectiveness of the method.