• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.597-600
• /
• 2008

Preconditioners for a two-dimensional combined finite element formulation have been devised and tested for fluid-structure interaction (FSI) problems. The FSI code simulating the interaction of a elastic body with an unsteady flow is based on P2P1 finite element combined formulation. It has been shown that two preconditioners among them perform well with respect to computational memory and convergence for a bench-mark problem. Based on the verification of the preconditioners for the two-dimensional combined formulation, four preconditioners are proposed for the problem of an elastic body interacting with a flow.

• Interaction and multiscale mechanics
• /
• v.6 no.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.

• Advances in aircraft and spacecraft science
• /
• v.1 no.2
• /
• pp.145-175
• /
• 2014

The results of a series of numerical experiments are presented to verify some of the important developments made in the first part of this paper. Firstly, the static solution of an algebraic system obtained through Strong Formulation Finite Element Method (SFEM) is presented. Secondly, the stress and strain recovery procedure is descripted for the present technique. It will be clear that the present approach is suitable for any strong formulation finite element methodology, due to the presented general approach based on the unknown displacements and on the elasticity equations. Thirdly, the numerical solutions for some classical and other numerical results found in literature are exposed. Finally, an arbitrarily shaped composite plate is solved and good agreement is observed for all the presented cases.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.5 no.6
• /
• pp.141-147
• /
• 1997

For the effective finite element analysis of the structures including material interfaces or contact surfaces, interface elements are proposed. Most of early works in this problem require not only iterative computation but also complex formulation because of the kinematic nonlinearities caused from the discontinuous behavior and the stress concentration phenomena. The proposed elements, however, are consistently formulated using relative displacements and tractions between top and bottom regular finite elements. The effectiveness of these elements are shown by solving various numerical sample problems including an leaf spring and comparing with results of general finite element analysis. As a result, more stable solutions are conveniently obtaines using interface elements than regular finite elements.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.83-105
• /
• 2013

In this paper, a multi-scale meshfree-enriched finite element formulation is presented for the analysis of acoustic wave propagation problem. The scale splitting in this formulation is based on the Variational Multi-scale (VMS) method. While the standard finite element polynomials are used to represent the coarse scales, the approximation of fine-scale solution is defined globally using the meshfree enrichments generated from the Generalized Meshfree (GMF) approximation. The resultant fine-scale approximations satisfy the homogenous Dirichlet boundary conditions and behave as the "global residual-free" bubbles for the enrichments in the oscillatory type of Helmholtz solutions. Numerical examples in one dimension and two dimensional cases are analyzed to demonstrate the accuracy of the present formulation and comparison is made to the analytical and two finite element solutions.

• Journal of the Korean Geotechnical Society
• /
• v.16 no.4
• /
• pp.191-199
• /
• 2000

In this research, the finite element analysis of piezocone penetration and dissipation tests has been conducted using the anisotropic elastoplastic-viscoplastic bounding surface model, virtual work equation, and theory of mixtures formulated in the Up[dated Lagrangian reference frame for the large deformation and finite strain nature of piezocone penetration. The formulated equations have been implemented into a finite element program. The cone resistance, excess pore water pressure, and dissipation of excess pore water pressure from the finite element analysis have been compared and investigated. An effective simulation could be performed with the use of the anisotropic and viscous soil model. The finite element formulations and the results are described in part 'I' and part 'II' respectively.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.5
• /
• pp.820-830
• /
• 1989

The tire inflation and contact problem has been solved by a finite element method. The finite element formulation is derived from the equilibrium equations by the principle of virtual work in the form of an updated Lagrangian formulation for incremental analysis. Then, a contact formulation is added to the finite element formulation to calculate stress state of tire in contact with flat rigid road under the load due to the self-weight of a vehicle. In the finite element analysis, equations of effective material properties are introduced to analyze a plane strain model of the shell-like tire by considering the bending effect of reinforced steel cords. The proposed equations of effective material properties produced stress concentration around the edge of belt layers, which does not appear when other well-known equations of material properties are adopted. The result from the above algorithm demonstrates the validity of the formulation and the proposed equations for the effective elastic constants. The result fully interprets the cause of separation between belt layers by showing the stress concentration.

• Korean Journal of Computational Design and Engineering
• /
• v.15 no.6
• /
• pp.411-417
• /
• 2010

In this paper, the boom of a floating crane is modeled as a 3-dimensional elastic beam in order to analyze the dynamic response of the crane and its cargo. The boom is divided into more than two elements based on finite element formulation, and deformation of each element is expressed in terms of shape matrix and nodal coordinates. The equations of motion for the elastic boom consist of a mass matrix, a stiffness matrix, and a quadratic velocity vector that contains the gyroscopic and Coriolis forces. The size and complicity of the matrices increase in proportion with the number of elements. Therefore, it is not possible to derive the equations of motion explicitly for different number of elements. To overcome this difficulty, matrices for one 3-dimensional element are expressed with elementary sub-matrices. In particular, the quadratic velocity vector is derived as a product of a shape matrix and a 3-dimensional rotation matrix. By using the derived matrices, the equations of motion for the multi-element boom are automatically constructed. To verify the implementation of the elastic boom based on finite element formulation, we simulated a simple vibration of the elastic boom and compared the average deformation with the analytic solution. Finally, heave motion of the floating crane and surge motion of the cargo are presented as application examples of the elastic boom.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1992.10a
• /
• pp.117-123
• /
• 1992

In recent years, the finite element method has become one of the most popular numerical technique for obtaining solutions of engineering science problems. However, there exist various uncertainties in modeling the problems, such as the dimensions(geometry shape), the material properties, boundary conditions, etc. The consideration for the uncertainties inherent in the problems can be made by understanding the influences of uncertain parameters[1]. Determining the influences of uncertainties as statistical quantities using the standard finite element method requires enormous computing time, while the probabilistic finite element method is realized as an efficient scheme[2,3] yielding statistical solution with just a few direct computations. In this paper, a formulation of the probabilistic fluid-structure interaction problem accounting for the first order perturbation of geometric shape is derived, and especially probabilistical acoustic pressure scattering from the structure with surrounding fluid is focused on. In Section 2, governing equations for the fluid-structure problems are given. In Section 3, a finite element formulation, based on the functional, is presented. First order perturbation of geometric shape with randomness is incorporated into the finite element formulation in conjunction with discretization of the random fields in Section 4 and 5. Finally, the proposed formulation is applied to a acoustic pressure scattering problem from an infinitely long cylindrical shell structure with randomness of radial perturbation.

• Structural Engineering and Mechanics
• /
• v.69 no.6
• /
• pp.615-626
• /
• 2019

The 3-noded metric Timoshenko beam element with an offset of the internal node from the element centre is used here to demonstrate the best-fit paradigm using function space formulation under locking and mesh distortion. The best-fit paradigm follows from the projection theorem describing finite element analysis which shows that the stresses computed by the displacement finite element procedure are the best approximation of the true stresses at an element level as well as global level. In this paper, closed form best-fit solutions are arrived for the 3-noded Timoshenko beam element through function space formulation by combining field consistency requirements and distortion effects for the element modelled in metric Cartesian coordinates. It is demonstrated through projection theorems how lock-free best-fit solutions are arrived even under mesh distortion by using a consistent definition for the shear strain field. It is shown how the field consistency enforced finite element solution differ from the best-fit solution by an extraneous response resulting from an additional spurious force vector. However, it can be observed that when the extraneous forces vanish fortuitously, the field consistent solution coincides with the best-fit strain solution.