• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.46 no.3
• /
• pp.210-218
• /
• 2018

In order to generate meshes automatically for finite element analysis of complex structures such as aircraft, a large number of triangular elements are typically created. However, triangular elements are less accurate than rectangular elements, so it is difficult to obtain a reliable solution. This problem can be improved through the meshfree method using the back cell integration. However, this method also causes some problems such as over-use of the integration points and inefficiency of the integral domain. In order to improve these problems, a method of performing integration by setting the integral area based on a node basis has been proposed, but in the case of incompressible material problems, the numerical accuracy deteriorates due to the vibration phenomenon of the solution. Therefore, in this paper, the modified meshfree method is proposed which sets the integral domain as an element domain instead of the nodal domain, and the proposed method improves the numerical instability caused by the conventional meshfree method without decreasing the accuracy regardles of the shape of integral domain. The effectiveness of the modified meshfree method is verified by using 2-D examples.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.37 no.4
• /
• pp.555-567
• /
• 2013

A finite element formulation for a Rayleigh-damped Bernoulli-Euler beam subjected to support motions, which accompanies quasi-static rigid-body motion, is presented by using the quasi-static decomposition method. Moving support beam elements, one of whose nodes is coincident with the moving support, are developed to represent the effect of a moving support. Statically determinate beams subjected to support motions can be modeled successfully by using moving support elements. Examples of cantilever and simply-supported beams subjected to support motions are illustrated, and the numerical results are compared with the analytical solutions. The comparison shows good agreement.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.29 no.4 s.235
• /
• pp.534-539
• /
• 2005

Nodal displacements are referred to the initial configuration in the total Lagrangian formulation and to the last converged configuration in the updated Lagrangian furmulation. This research proposes a relative nodal displacement method to represent the position and orientation for a node in truss structures. Since the proposed method measures the relative nodal displacements relative to its adjacent nodal reference frame, they are still small for a truss structure undergoing large deformations for the small size elements. As a consequence, element formulations developed under the small deformation assumption are still valid for structures undergoing large deformations, which significantly simplifies the equations of equilibrium. A structural system is represented by a graph to systematically develop the governing equations of equilibrium for general systems. A node and an element are represented by a node and an edge in graph representation, respectively. Closed loops are opened to form a spanning tree by cutting edges. Two computational sequences are defined in the graph representation. One is the forward path sequence that is used to recover the Cartesian nodal displacements from relative nodal displacement sand traverses a graph from the base node towards the terminal nodes. The other is the backward path sequence that is used to recover the nodal forces in the relative coordinate system from the known nodal forces in the absolute coordinate system and traverses from the terminal nodes towards the base node. One open loop and one closed loop structure undergoing large deformations are analyzed to demonstrate the efficiency and validity of the proposed method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.2
• /
• pp.151-158
• /
• 2005

This research analyzes eigenvalue analysis of stiffened plates on the Pasternak foundations using the finite clement method. For analyzing the stiffened plates, both the Mindlin plate theory and Timoshenko beam-column theory were applied. In application of the finite element method, 8-nodes serendipity clement system and 3-nodes finite element system were used for plate and beam elements, respectively. Elastic foundations were modeled as the Pasternak foundations in which the continuity effect of foundations is considered. In order to verify the theory of this study, solutions obtained by this analysis were compared with the classical solutions in reference, experimental solutions and solutions obtained by SAP 2000. The natural frequency of stiffened plates on Pasternak foundations were determined according to changes or foundation parameters and dimensions of stiffener.

• Journal of the Society of Naval Architects of Korea
• /
• v.53 no.5
• /
• pp.388-399
• /
• 2016

A motion analysis of an underwater towed cable is a complex task due to its nonlinear nature of the problem. The major source of the nonlinearity of the underwater cable analysis is that the motion of the cable involves large rigid-body motion. This large rigid-body motion makes difficult to use standard displacement-based finite element method. In this paper, the authors apply recently developed nodal position-based finite element method which can deal with the geometric nonlinearity due to the large rigid-body motion. In order to enhance the stability of the large-scale nonlinear cable motion simulation, an efficient time-integration scheme is proposed, namely predictor/multi-corrector Newmark scheme. Three different predictors are introduced, and the best predictor in terms of stability and robustness for impulsive cable motion analysis is proposed. As a result, the nonlinear motion of underwater cable is predicted in a very efficient manner compared to the classical finite element of finite difference methods. The efficacy of the method is demonstrated with several test cases, involving static and dynamic motion of a single cable element, and also under water towed cable composed of multiple cable elements.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.1
• /
• pp.24-34
• /
• 1990

In this paper, the formulation of a new simplified finite element is made to analyze the ultimate strength of damaged tubular members subjected to combined axial force and end moment. A damaged tubular member that has the bending deformation and the local dent is modeled by beam elements. Tangent elastic stiffness matrix of a beam element which contains the effect of the geometric nonlinearity is derived by using the updated Lagrangian approach. Here the contribution of the stiffness in the dented area is neglected since its resistance against the external loads is considered to be small. A fully plastic interaction curve of the element under combined loads taking account of the local dent effect is selected as a yielding criterion at each nodal point. Also tangent elasto-plastic stiffness matrix of the element is formulated by plastic node method. Comparison with the present solution and the existing experimental results is made showing that the present method gives quite an accurate solution.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2005.11a
• /
• pp.788-791
• /
• 2005

Nodal displacements are referred to the Initial configuration in the total Lagrangian formulation and to the last converged configuration in the updated Lagrangian formulation. This research proposes a relative nodal displacement method to represent the position and orientation for a node in truss structures. Since the proposed method measures the relative nodal displacements relative to its adjacent nodal reference frame, they are still small for a truss structure undergoing large deformations for the small size elements. As a consequence, element formulations developed under the small deformation assumption are still valid fer structures undergoing large deformations, which significantly simplifies the equations of equilibrium. A structural system is represented by a graph to systematically develop the governing equations of equilibrium for general systems. A node and an element are represented by a node and an edge in graph representation, respectively. Closed loops are opened to form a spanning tree by cutting edges. Two computational sequences are defined in the graph representation. One is the forward path sequence that is used to recover the Cartesian nodal displacements from relative nodal displacements and traverses a graph from the base node towards the terminal nodes. The other is the backward path sequence that is used to recover the nodal forces in the relative coordinate system from the known nodal forces in the absolute coordinate system and traverses from the terminal nodes towards the base node. One closed loop structure undergoing large deformations is analyzed to demonstrate the efficiency and validity of the proposed method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.3
• /
• pp.321-327
• /
• 2007

This study presents a new gridless finite difference method for solving elastic crack problems. The method constructs the Taylor expansion based on the MLS(Moving Least Squares) method and effectively calculates the approximation and its derivatives without differentiation process. Since no connectivity between nodes is required, the modeling of discontinuity embedded in the domain is very convenient and discontinuity effect due to crack is naturally implemented in the construction of difference equations. Direct discretization of the governing partial differential equations makes solution process faster than other numerical schemes using numerical integration. Numerical results for mode I and II crack problems demonstrates that the proposed method accurately and efficiently evaluates the stress intensity factors.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2009.05a
• /
• pp.815-821
• /
• 2009

하천의 2차원 흐름 및 하상변동, 오염확산 해석을 위한 유체의 수치해석법에는 유한요소법, 유한차분법, 유한차분법의 변형인 유한체적법, 경계적분법 등이 있으며, 국내의 경우 비구조적 요소망(unstructured mesh)을 이용하여 복잡한 형상을 표현하기가 상대적으로 용이한 유한요소법이 널리 사용되고 있다. 하천을 유한 요소화 하는 전처리 과정은 전체 해석 과정을 자동화 하는데 있어 필수적인 요소이며, 주로 삼각 요소망 또는 사각 요소망을 이용하여 해석을 수행하게 된다. 삼각 요소망의 경우 상대적으로 자동화하기 쉬운 반면 사각 요소망의 생성은 절점 생성 자체가 삼각 요소망 보다 더 많은 기하학적 제한 요소를 가지고 있기 때문에 상대적으로 완성도 높은 알고리즘을 구현하기가 어렵다 할 수 있다. 이에 따라 본 연구에서는 2차원 상에서 사각 요소망(quadrilateral elements)을 생성할 수 있는 Paving method를 중심으로 한 요소망 생성 알고리즘에 대해 고찰하고, 국내 최초의 범용 수치해석 모형인 RAMS(River Analysis and Modeling System)에 적용하였다. Paving method는 1990년에 Blacker and Stephenson에 의해 제안되었으며, Sandia National Laboratories에 의해 완성되었다. Paving Method는 advancing front style의 요소망을 생성하게 되고, 바깥쪽에서 안쪽으로 element layer를 생성하면서 채워나간다. 본 연구에서는 기존의 요소망 생성 프로세스에서 element 삽입 전의 검증 기능을 강화한 새로운 버전의 paving method를 적용하엿다.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.33 no.4
• /
• pp.310-320
• /
• 2009

In this paper, we estimate area-averaged solution of peening residual stress using a 3-D multi-impact symmetry-cell FE model. The symmetry-cell model includes factors reflecting peening phenomena and plastic shot. Area-averaged solution is much closer to XRD experimental solution than 4-node-averaged solution in plastic shot FE model. We then obtain FE Almen saturation curve corresponding to experimental Almen curve based on area-averaged solution. Using the curve, we obtain FE area-averaged solution in major peening materials, and compare the FE solution with experimental solution. In peening materials, surface, maximum compressive residual stress and deformation depth reach experimental solutions. Thus, FE Almen curve is useful for estimation of residual stress solution and could improve the efficiency of peening process. Consequently, it is confirmed that concept of area-averaged solution is the realistic analytical method for evaluation of peening residual stress.