• Journal of Theoretical and Applied Mechanics
• /
• v.1 no.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
• /
• v.20 no.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
• /
• v.19 no.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
• /
• v.20 no.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.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.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.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.23 no.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.

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

• Journal of Mechanical Science and Technology
• /
• v.14 no.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
• /
• v.47 no.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.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.5
• /
• pp.385-392
• /
• 2013

The finite element method has become the most widely used method of structural analysis and recently, the method has often been applied to complex dynamic and nonlinear structural analyses problems. Even for these complex problems, where the responses are hard to predict, finite element analyses yield reliable results if appropriate element types and meshes are used. However, the dynamic and nonlinear behaviors of a structure often include large deformations in various portions of the structure and if the same mesh is used throughout the analysis, some elements may deform to shapes beyond the reliable limits; thus dynamically adapting finite element meshes are needed in order for the finite element analyses to be accurate. In addition, to satisfy the users requirement of quick real run time of finite element programs, the algorithms must be computationally efficient. This paper presents an adaptive finite element mesh generation scheme for dynamic analyses of structures that may adapt at each time step. Representative strain values are used for error estimates and combinations of the h-method(node movement) and the r-method(element division) are used for mesh refinements. A coefficient that depends on the shape of an element is used to limit overly distorted elements. A simple frame example shows the accuracy and computational efficiency of the scheme. The aim of the study is to outline the adaptive scheme and to demonstrate the potential use in general finite element analyses of dynamic and nonlinear structural problems commonly encountered.