• 한국전산구조공학회논문집
• /
• 제36권5호
• /
• pp.295-305
• /
• 2023

본 연구에서는 평활화 유한요소법(Smoothed finite element method)을 도입한 위상분야법(Phase-field method)에 대해 소개하고자 한다. 위상분야법은 최근 균열 개시 및 전파 해석에 많이 사용되는 기법으로 균열 표면을 추적하기 위한 추가적인 처리기법이 필요하지 않는 특징이 있다. 위상분야법에서 복잡한 균열 전파를 포착하기 위해 높은 정확도의 변형률 에너지를 평활화 유한요소법을 도입하여 계산하였다. 평활화 유한요소법은 유한요소를 하위 셀로 나누고 각각의 하위 셀을 평활화 영역으로 재조립하여 변형률 에너지를 계산하게 된다. 또한 해석 시간 단축을 위하여 쿼드트리 요소망을 제안한 기법에 사용하였다. 수치 예제를 통하여 제안한 기법을 참조해 및 유한요소법과 비교하여 검증하였다.

• Structural Engineering and Mechanics
• /
• 제86권3호
• /
• pp.311-324
• /
• 2023

Herein, we present effective polygonal finite elements to which the strain-smoothed element (SSE) method is applied. Recently, the SSE method has been developed for conventional triangular and quadrilateral finite elements; furthermore, it has been shown to improve the performance of finite elements. Polygonal elements enable various applications through flexible mesh handling; however, further development is still required to use them more effectively in engineering practice. In this study, piecewise linear shape functions are adopted, the SSE method is applied through the triangulation of polygonal elements, and a smoothed strain field is constructed within the element. The strain-smoothed polygonal elements pass basic tests and show improved convergence behaviors in various numerical problems.

• 대한건축학회논문집:구조계
• /
• 제35권9호
• /
• pp.159-169
• /
• 2019

This work presents the three-dimensional extended strain smoothing approach in the framework of finite element method, so-called smoothed finite element method (S-FEM) for quasi-incompressible hyperelastic materials undergoing the large deformations. The proposed method is known that the incompressible limits, such as over-estimation of stiffness and distorted mesh sensitivity, can be overcome in two dimensions. Therefore, in this paper, the idea of Cell-based, Edge-based and Node-based strain smoothing approaches is extended to three-dimensions. The construction of subcells and smoothing domains for each methods are explained. The smoothed strain-displacement matrix and the stiffness matrix are obtained on each smoothing domain in the same manner with two-dimensional S-FEM. Various numerical tests are studied to demonstrate the validity and accuracy of 3D-S-FEM. The obtained results are compared with analytical solutions to express the efficacy of the methods.

• Structural Engineering and Mechanics
• /
• 제39권6호
• /
• pp.767-793
• /
• 2011

An efficient edge-based smoothed finite element method (ES-FEM) has been recently developed for solving solid mechanics problems. The ES-FEM uses triangular elements that can be generated easily for complicated domains. In this paper, the complexity study of the ES-FEM based on triangular elements is conducted in detail, which confirms the ES-FEM produces higher computational efficiency compared to the FEM. Therefore, the ES-FEM offers an excellent platform for adaptive analysis, and this paper presents an efficient adaptive procedure based on the ES-FEM. A smoothing domain based energy (SDE) error estimate is first devised making use of the features of the ES-FEM. The present error estimate differs from the conventional approaches and evaluates error based on smoothing domains used in the ES-FEM. A local refinement technique based on the Delaunay algorithm is then implemented to achieve high efficiency in the mesh refinement. In this refinement technique, each node is assigned a scaling factor to control the local nodal density, and refinement of the neighborhood of a node is accomplished simply by adjusting its scaling factor. Intensive numerical studies, including an actual engineering problem of an automobile part, show that the proposed adaptive procedure is effective and efficient in producing solutions of desired accuracy.

• Smart Structures and Systems
• /
• 제11권3호
• /
• pp.241-259
• /
• 2013

This paper presents the optimum load-speed diagram evaluation for a linear micromotor, including multitude cantilever piezoelectric bimorphs, briefly. Each microbeam in the mechanism can be actuated in both axial and flexural modes simultaneously. For this design, we consider quasi-static and linear conditions, and a relatively new numerical method called the smoothed finite element method (S-FEM) is introduced here. For this purpose, after finding an optimum volume fraction for piezoelectric layers through a standard numerical method such as quadratic finite element method, the relevant load-speed curves of the optimized micromotor are examined and compared by deterministic topology optimization (DTO) design. In this regard, to avoid the overly stiff behavior in FEM modeling, a numerical method known as the cell-based smoothed finite element method (CS-FEM, as a branch of S-FEM) is applied for our DTO problem. The topology optimization procedure to find the optimal design is implemented using a solid isotropic material with a penalization (SIMP) approximation and a method of moving asymptotes (MMA) optimizer. Because of the higher efficiency and accuracy of S-FEMs with respect to standard FEMs, the main micromotor characteristics of our final DTO design using a softer CS-FEM are substantially improved.

• 한국전산구조공학회논문집
• /
• 제37권3호
• /
• pp.187-195
• /
• 2024

본 논문에서는 평활화 유한요소법(Smoothed finite element method)을 도입한 응력 기반 구배 탄성론(Gradient elasticity)의 2차원 경계치 문제에 대한 연구를 수행하였다. 구배 탄성론은 기존 탄성론에서는 표현할 수 없는 미소규모의 크기 의존적인 기계적 거동을 설명하기 위해 제안되었다. 구배 탄성체론에서 고차 미분 방정식을 두 개의 2차 미분 방정식으로 분할하는 Ru-Aifantis 이론을 사용하기 때문에 평활화 유한요소법에 적용이 가능하게 된다. 본 연구에서 경계치 문제를 해결하기 위해 평활화 유한 요소 프레임워크에 스태거드 방식(Staggered scheme)을 사용하여 국부 변위장과 비국부 응력장을 평활화 영역 및 요소에서 각각 계산하였다. 구배 탄성에서 중요한 변수인 내부 길이 척도의 영향을 측정하기 위해 일련의 수치 예제를 수행하였다. 수치 해석 결과는 제안한 기법이 내부 길이 척도에 따라 균열 선단과 전위 선에 나타나는 응력 집중을 완화할 수 있음을 보여준다.

• Geomechanics and Engineering
• /
• 제11권1호
• /
• pp.1-39
• /
• 2016

The finite element method (FEM), discrete element method (DEM), and Discontinuous deformation analysis (DDA) are among the standard numerical techniques applied in computational geo-mechanics. However, in some cases there no possibility for modelling by traditional finite analytical techniques or other mesh-based techniques. The solution presented in the current study as a completely Lagrangian and mesh-free technique is smoothed particle hydrodynamics (SPH). This method was basically applied for simulation of fluid flow by dividing the fluid into several particles. However, several researchers attempted to simulate soil-water interaction, landslides, and failure of soil by SPH method. In fact, this method is able to deal with behavior and interaction of different states of materials (liquid and solid) and multiphase soil models and their large deformations. Soil indicates different behaviors when interacting with water, structure, instrumentations, or different layers. Thus, study into these interactions using the mesh based grids has been facilitated by mesh-less SPH technique in this work. It has been revealed that the fast development, computational sophistication, and emerge of mesh-less particle modeling techniques offer solutions for problems which are not modeled by the traditional mesh-based techniques. Also it has been found that the smoothed particle hydrodynamic provides advanced techniques for simulation of soil materials as compared to the current traditional numerical methods. Besides, findings indicate that the advantages of applying this method are its high power, simplicity of concept, relative simplicity in combination of modern physics, and particularly its potential in study of large deformations and failures.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제12권1호
• /
• pp.85-101
• /
• 2020

The optimal parameters for the fluid-structure interaction analysis using the Smoothed Particle Hydrodynamics (SPH) for fluids and finite elements for structures, respectively, are explored, and the effectiveness of the simulations with those parameters is validated by solving several open surface fluid problems. For the optimization of the Equation of State (EOS) and the simulation parameters such as the time step, initial particle spacing, and smoothing length factor, a dam-break problem and deflection of an elastic plate is selected, and the least squares analysis is performed on the simulation results. With the optimal values of the pivotal parameters, the accuracy of the simulation is validated by calculating the exerted force on a moving solid column in the open surface fluid. Overall, the SPH-FEM coupled simulation is very effective to calculate the fluid-structure interaction. However, the relevant parameters should be carefully selected to obtain accurate results.

• Interaction and multiscale mechanics
• /
• 제6권2호
• /
• pp.237-255
• /
• 2013

In this paper, a meshfree-enriched finite element method (ME-FEM) is introduced for the large deformation analysis of nonlinear path-dependent problems involving contact. In linear ME-FEM, the element formulation is established by introducing a meshfree convex approximation into the linear triangular element in 2D and linear tetrahedron element in 3D along with an enriched meshfree node. In nonlinear formulation, the area-weighted smoothing scheme for deformation gradient is then developed in conjunction with the meshfree-enriched element interpolation functions to yield a discrete divergence-free property at the integration points, which is essential to enhance the stress calculation in the stage of plastic deformation. A modified variational formulation using the smoothed deformation gradient is developed for path-dependent material analysis. In the industrial metal forming problems, the mortar contact algorithm is implemented in the explicit formulation. Since the meshfree-enriched element shape functions are constructed using the meshfree convex approximation, they pose the desired Kronecker-delta property at the element edge thus requires no special treatments in the enforcement of essential boundary condition as well as the contact conditions. As a result, this approach can be easily incorporated into a conventional displacement-based finite element code. Two elasto-plastic problems are studied and the numerical results indicated that ME-FEM is capable of delivering a volumetric locking-free and pressure oscillation-free solutions for the large deformation problems in metal forming analysis.

• Nuclear Engineering and Technology
• /
• 제48권1호
• /
• pp.246-258
• /
• 2016

The shield building of AP1000 was designed to protect the steel containment vessel of the nuclear reactor. Therefore, the safety and integrity must be ensured during the plant life in any conditions such as an earthquake. The aim of this paper is to study the effect of water in the water tank on the response of the AP1000 shield building when subjected to three-dimensional seismic ground acceleration. The smoothed particle hydrodynamics method (SPH) and finite element method (FEM) coupling method is used to numerically simulate the fluid and structure interaction (FSI) between water in the water tank and the AP1000 shield building. Then the grid convergence of FEM and SPH for the AP1000 shield building is analyzed. Next the modal analysis of the AP1000 shield building with various water levels (WLs) in the water tank is taken. Meanwhile, the pressure due to sloshing and oscillation of the water in the gravity drain water tank is studied. The influences of the height of water in the water tank on the time history of acceleration of the AP1000 shield building are discussed, as well as the distributions of amplification, acceleration, displacement, and stresses of the AP1000 shield building. Research on the relationship between the WLs in the water tank and the response spectrums of the structure are also taken. The results show that the high WL in the water tank can limit the vibration of the AP1000 shield building and can more efficiently dissipate the kinetic energy of the AP1000 shield building by fluid-structure interaction.