• Journal of Theoretical and Applied Mechanics
• /
• 제1권1호
• /
• pp.1-27
• /
• 1995

In this paper a finite element method for free-surface problems is described. the method is based on two different forms of Hamilton's principle. To test the present computational method two specific wave problems are investigated; the dispersion relations and the nonlinear effect for the well-known solitary waves are treated. The convergence test shows that the present scheme is more efficient than other existing methods, e.g. perturbation scheme.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제20권3호
• /
• pp.243-259
• /
• 2016

The Richards equation for water movement in unsaturated soil is highly nonlinear partial differential equations which are not solvable analytically unless unrealistic and oversimplifying assumptions are made regarding the attributes, dynamics, and properties of the physical systems. Therefore, conventionally, numerical solutions are the only feasible procedures to model flow in partially saturated porous media. The standard Finite element numerical technique is usually coupled with an Euler time discretizations scheme. Except for the fully explicit forward method, any other Euler time-marching algorithm generates nonlinear algebraic equations which should be solved using iterative procedures such as Newton and Picard iterations. In this study, lumped mass and distributed mass in the frame of Picard and Newton iterative techniques were evaluated to determine the most efficient method to solve the Richards equation with finite element model. The accuracy and computational efficiency of the scheme and of the Picard and Newton models are assessed for three test problems simulating one-dimensional flow processes in unsaturated porous media. Results demonstrated that, the conventional mass distributed finite element method suffers from numerical oscillations at the wetting front, especially for very dry initial conditions. Even though small mesh sizes are applied for all the test problems, it is shown that the traditional mass-distributed scheme can still generate an incorrect response due to the highly nonlinear properties of water flow in unsaturated soil and cause numerical oscillation. On the other hand, non oscillatory solutions are obtained and non-physics solutions for these problems are evaded by using the mass-lumped finite element method.

• Journal of applied mathematics & informatics
• /
• 제19권1_2호
• /
• pp.385-395
• /
• 2005

A finite element scheme is considered for the viscous Cahn-Hilliard equation with the nonconstant gradient energy coefficient. The scheme inherits energy decay property and mass conservation as for the classical solution. We obtain the corresponding error estimate using the extended Lax-Richtmyer equivalence theorem.

• Structural Engineering and Mechanics
• /
• 제20권4호
• /
• pp.405-420
• /
• 2005

In this paper, a hybrid/mixed nonlinear shell element is developed in polar coordinate system based on Hellinger/Reissner variational principle and the large-deflection theory of plate. A numerical solution scheme is formulated using the hybrid/mixed finite element method (HMFEM), in which the nodal values of bending moments and the deflection are the unknown discrete parameters. Stability of the present element is studied. The large-deflection analyses are performed for simple supported and clamped circular plates under uniformly distributed and concentrated loads using HMFEM and the traditional displacement finite element method. A parametric study is also conducted in the research. The accuracy of the shell element is investigated using numerical computations. Comparisons of numerical solutions are made with theoretical results, finite element analysis and the available numerical results. Excellent agreements are shown.

• 한국전산구조공학회논문집
• /
• 제23권6호
• /
• pp.659-665
• /
• 2010

The finite element method(FEM) is proven to be an effective approximate method of structural analysis if proper element types and meshes are chosen, and recently, the method is often applied to solve complex dynamic and nonlinear problems. A properly chosen element type and mesh yields reliable results for dynamic finite element structural analysis. However, dynamic behavior of a structure may include unpredictably large strains in some parts of the structure, and using the initial mesh throughout the duration of a dynamic analysis may include some elements to go through strains beyond the elements' reliable limits. Thus, the finite element mesh for a dynamic analysis must be dynamically adaptive, and considering the rapid process of analysis in real time, the dynamically adaptive finite element mesh generating schemes must be computationally efficient. In this paper, a computationally efficient dynamically adaptive finite element mesh generation scheme for dynamic analyses of structures is described. The concept of representative strain value is used for error estimates and the refinements of meshes use combinations of the h-method(node movement) and the r-method(element division). The shape coefficient for element mesh is used to correct overly distorted elements. The validity of the scheme is shown through a cantilever beam example under a concentrated load with varying values. The example shows reasonable accuracy and efficient computing time. Furthermore, the study shows the potential for the scheme's effective use in complex structural dynamic problems such as those under seismic or erratic wind loads.

• 한국지진공학회논문집
• /
• 제23권1호
• /
• pp.83-88
• /
• 2019

Structural analysis remains as an essential part of any integrated civil engineering system in today's rapidly changing computing environment. Even with enormous advancements in capabilities of computers and mobile tools, enhancing computational efficiency of algorithms is necessary to meet the changing demands for quick real time response systems. The finite element method is still the most widely used method of computational structural analysis; a robust, reliable and automated finite element structural analysis module is essential in a modern integrated structural engineering system. To be a part of an automated finite element structural analysis, an efficient adaptive mesh generation scheme based on R-H refinement for the mesh and error estimates from representative strain values at Gauss points is described. A coefficient that depends on the shape of element is used to correct overly distorted elements. Two simple case studies show the validity and computational efficiency. The scheme is appropriate for nonlinear and dynamic problems in earthquake engineering which generally require a huge number of iterative computations.

• 대한조선학회논문집
• /
• 제53권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.

• Journal of Mechanical Science and Technology
• /
• 제14권6호
• /
• pp.666-674
• /
• 2000

Large deflection dynamic responses of laminated composite cylindrical shells under impact are analyzed by the geometrically nonlinear finite element method based on a generalized Sander's shell theory with the first order transverse shear deformation and the von-Karman large deflection assumption. A modified indentation law with inelastic indentation is employed for the contact force. The nonlinear finite element equations of motion of shell and an impactor along with the contact laws are solved numerically using Newmark's time marching integration scheme in conjunction with Akay type successive iteration in each step. The ply failure region of the laminated shell is estimated using the Tsai- Wu quadratic interaction criteria. Numerical results, including the contact force histories, deflections and strains are presented and compared with the ones by linear analysis. The effect of the radius of curvature on the composite shell behaviors is investigated and discussed.

• Structural Engineering and Mechanics
• /
• 제47권1호
• /
• pp.27-44
• /
• 2013

A finite element computational procedure for the accurate analysis of quasistatic thermorheological complex structures response is developed. The geometrical nonlinearity, arising from large displacements and rotations (but small strains), is accounted for by the total Lagrangian description of motion. The Schapery's nonlinear single-integral viscoelastic constitutive model is modified for a time-stress-temperature-dependent behavior. The nonlinear thermo-viscoelastic constitutive equations are incrementalized leading to a recursive relationship and thereby the resulting finite element equations necessitate data storage from the previous time step only, and not the entire deformation history. The Newton-Raphson iterative scheme is employed to obtain a converged solution for the non-linear finite element equations. The developed numerical model is verified with the previously published works and a good agreement with them is found. The applicability of the developed model is demonstrated by analyzing two examples with different thermal/mechanical loading histories.

• 한국전산구조공학회논문집
• /
• 제26권5호
• /
• pp.385-392
• /
• 2013

유한요소법은 구조해석법으로 가장 많이 사용되는 방법으로 자리잡고 있으며, 근래에는 다소 복잡한 동적 및 비선형 문제에도 사용이 일반화되고 있다. 이러한 거동 예측이 어려운 구조해석에도 구조물을 적절한 유한요소와 요소망으로 표현하면 신뢰있는 해석 결과를 얻을 수 있다. 구조물의 동적 또는 비선형 거동에는 예상하지 않은 부분에서 큰 변형이 일어날 수 있으며, 유한요소해석 과정에서 같은 요소망을 계속 사용하면 요소의 모양이 신뢰 범위 밖으로 변형될 수 있으므로 요소망 역시 동적으로 적응할 필요가 있다. 또한, 유한요소 프로그램의 사용자 요구 사항 중 하나가 실시간으로 빠르게 진행되는 것이므로 연산면에서 효율적이어야 한다. 본 연구는 시간영역 동적해석에서 전 단계 해석 결과를 사용하여 계산된 대표 변형률값을 오차 평가에 사용하여 절점 이동인 r-법과 요소 분할인 h-법의 조합으로 요소 세분화를 진행하여 동적으로 적응하는 요소망 형성 과정을 기술한다. 해석 중 과대하게 변형되는 요소는 모양계수 개념으로 방지한다. 간단한 프레임의 동적 유한요소해석을 예제로 정확성과 연산 효율성을 보여준다. 본 연구에서 제시하는 적응적 유한요소망 형성 전략은 복잡한 동적 및 비선형 해석에 일반적으로 적용될 수 있다.