• 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.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제8권2호
• /
• pp.153-168
• /
• 2016

Sloshing phenomenon is a complicated free surface flow problem that increases the dynamic pressure on the sidewalls and the bottom of the storage tanks. When the storage tanks are partially filled, it is essential to be able to evaluate the fluid dynamic loads on the tank's perimeter. In this paper, a numerical code was developed to determine the pressure distribution on the rectangular and trapezoidal storage tanks' perimeters due to liquid sloshing phenomenon. Assuming the fluid to be inviscid, the Laplace equation and the nonlinear free surface boundary conditions were solved using coupled boundary element - finite element method. The code performance for sloshing modeling was validated using Nakayama and Washizu's results. Finally, this code was used for partially filled rectangular and trapezoidal storage tanks and free surface displacement, pressure distribution and horizontal and vertical forces exerted on the tanks' perimeters due to liquid sloshing phenomenon were estimated and discussed.

• 대한기계학회논문집B
• /
• 제21권9호
• /
• pp.1230-1243
• /
• 1997

A numerical technique for simulating incompressible viscous flow with free surface is presented. The flow field is obtained by penalty finite element formulation. In this work, a modified volume of fluid (VOF) method which is compatible with 4-node element is proposed to track the moving free surface. This scheme can be applied to irregular mesh system, and can be easily extended to three dimensional geometries. Numerical analyses were done for two benchmark examples, namely the broken dam problem and the solitary wave propagation problem. The numerical results were in close agreement with the existing data. Illustrative examples were studied to show the effectiveness of the proposed numerical scheme.

• 한국소음진동공학회논문집
• /
• 제13권9호
• /
• pp.663-670
• /
• 2003

As the micro-technologies in the high precision manufacturing processes are developed, the demand for micro actuating device is increasing. But, it is difficult to achieve high resolution and wide operating range simultaneously with the conventional actuating systems which are contacting and type of dual servo system. So, the contact-free surface actuators whose movers are suspended or levitated were proposed. These systems can be applied to high precision stages and alignment apparatuses. The suspended mover can be assumed to be rigid body, but the mover is a structure in this study, therefore the vibration caused during the operating process has a serious adverse effect on the performance and it is very important to identify the vibrational characteristics. In this paper, a contact-free surface actuator is modeled in finite element method and updated by using the experimental modal data. Finally, the static and dynamic characteristics of the finite element model are predicted and then discussed.

• Structural Engineering and Mechanics
• /
• 제86권6호
• /
• pp.765-778
• /
• 2023

This study investigates the free vibration and buckling analyses of isotropic curved plate structures fixed at all ends. The Kirchhoff-Love Plate Theory (KLPT) and Finite Element Method (FEM) are employed to model the curved structure. In order to perform the finite element analysis, a four-node quadrilateral element with 5 degrees of freedom (DOF) at each node is utilized. Additionally, the drilling effect (θ_z) is considered as minimal to satisfy the DOF of the structure. Lagrange's equation of motion is used in order to obtain the first ten natural frequencies and the critical buckling values of the structure. The effects of various radii of curvatures and aspect ratio on the natural frequency and critical buckling load values for the single-bay and two-bay curved frames are investigated within this scope. A computer code based on finite element analysis is developed to perform free vibration and buckling analysis of curved plate frames. The natural frequency and critical buckling load values of the present study are compared with ANSYS R18.2 results. It has been concluded that the results of the present study are in good agreement with ANSYS results for different radii of curvatures and aspect ratio values of both single-bay and two-bay structures.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제12권4호
• /
• pp.271-287
• /
• 2008

This paper deals with American call and put options. Determining the fair price and the free boundary of an American option is a very difficult problem since they depends on each other. This paper presents numerical algorithms of finite element method based on the three-level scheme to compute both the price and the free boundary. One algorithm is designed for American call options and the other one for American put options. These algorithms are formulated on the system of the Jamshidian equation for the option price and the free boundary. Here, the Jamshidian equation is of a kind of the nonhomogeneous Black-Scholes equations. We prove the existence and uniqueness of the numerical solution by the Lax-Milgram lemma and carried out extensive numerical experiments to compare with various methods.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2003년도 추계학술대회논문집
• /
• pp.977-983
• /
• 2003

Free vibration analysis of helical springs was performed by the use of the commercial finite element analysis program, ANSYS. The investigation of national frequency was focused on the effect of various parameters such as boundary conditions, spring indices, number of coil turns and helix angles which are considered to affect the free vibration of a spring. The finite element method was validated by comparison with the result of a previouosly published literature. The similarity of frequency trend was shown among three boundary conditions: clamped-clamped, free-free, simpliy supported-simply supported but there was no similarity in light of mode shapes among them. Several modes showed similar frequencies on and near the frequencies identified by the natural frequency formula of Wahl. Natural frequencies increased with spring indices and number of turns decreasing and with helix angles increasing. The results investigated by finiete element method were compared with the experemental result and theoretical result and showed a good agreement among them.

• Coupled systems mechanics
• /
• 제5권1호
• /
• pp.59-86
• /
• 2016

This paper deals with the free vibration analysis of a dynamical coupled system: flexible gravity dam- compressible rectangular reservoir. The finite element method is used to compute the natural frequencies and modal shapes of the system. Firstly, the reservoir and subsequently the dam is modeled by classical 8-node elements and the natural frequencies plus modal shapes are calculated. Afterwards, a new 21-node element is introduced and the same procedure is conducted in which an efficient method is employed to carry out the integration operations. Finally, the coupled dam-reservoir system is modeled by solely one 21-node element and the free vibration of dam-reservoir interaction system is investigated. As an important result, it is clearly concluded that the one high-order element treats more precisely than the eight-node elements, since the first one utilizes fifth-degree polynomials to construct the shape functions and the second implements polynomials of degree two.

• Structural Engineering and Mechanics
• /
• 제25권4호
• /
• pp.467-480
• /
• 2007

Two methods were developed for analyzing problems with two adjacent sub-domains modeled by different kinds of elements in finite element method. Each sub-domain can be defined independently without the consideration of equivalent division with common nodes used for the interface. These two methods employ an individual interface to accomplish the compatibility. The MLSA method uses the moving least square approximation which is the basic formulation for Element Free Galerkin Method to formulate the interface. The displacement field assumed by this method does not pass through nodes on the common boundary. Therefore, nodes can be chosen freely for this method. The results show that the MLSA method has better approximation than traditional methods.

• Structural Engineering and Mechanics
• /
• 제27권3호
• /
• pp.311-331
• /
• 2007

The purpose of this paper is to study shear locking-free analysis of thick plates using Mindlin's theory and to determine the effects of the thickness/span ratio, the aspect ratio and the boundary conditions on the linear responses of thick plates subjected to uniformly distributed loads. Finite element formulation of the equations of the thick plate theory is derived by using higher order displacement shape functions. 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, 8- and 17-noded quadrilateral finite elements are used. Graphs and tables are presented that should help engineers in the design of thick plates. 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. It is also concluded, in general, that the maximum displacement and bending moment increase with increasing aspect ratio, and that the results obtained in this study are better than the results given in the literature.