• Korean Journal of Hydrosciences
• /
• v.3
• /
• pp.45-59
• /
• 1992

The flow characteristics varying with rate of the radius of curvature to width (Rc/B) in open channel bends are investigated with a simplified numerical model. Secondary flow velocity and transverise bed slope are formulated from the equations of momentum and force balance analysis, respectively. The conservation equations of mass and streamwise momentum are simplified by depth integration and its solution could be obtained from the explicit finite difference method. Three sets of computer simulation are executed. The rates of Rc/B adopted in simulations are 2.7, 5.4 and 8.1. The terms analyzed in this paper secondary flow velocity, streamwise velocity, the path of maximum steamwise velocity, deviation angle, and mass-shift velocity.

• Structural Engineering and Mechanics
• /
• v.13 no.4
• /
• pp.387-410
• /
• 2002

Formulation of an 8 nodes assumed strain shell element is presented for the analysis of shells. The stiffness matrix based on the Mindlin-Reissner theory is analytically integrated through the thickness. The element is free of membrane and shear locking behavior by using the assumed strain method such that the element performs very well in modeling of thin shell structures. The material is assumed to be isotropic and laminated composite. The element has six degrees of freedom per node and can model the stiffened plates and shells. A great number of numerical testing carried out for the validation of present 8 node shell element are in good agreement with references.

• Structural Engineering and Mechanics
• /
• v.87 no.4
• /
• pp.363-373
• /
• 2023

Since the first family of structure-dependent methods can simultaneously integrate unconditional stability and explicit formulation in addition to second order accuracy, it is very computationally efficient for solving inertial problems except for adopting auto time-stepping techniques due to no nonlinear iterations. However, an unusual stability property is first found herein since its unconditional stability interval is drastically different for zero and nonzero damping. In fact, instability might occur for solving a damped stiffness hardening system while an accurate result can be obtained for the corresponding undamped stiffness hardening system. A technique of using a stability factor is applied to overcome this difficulty. It can be applied to magnify an unconditional stability interval. After introducing this stability factor, the formulation of this family of structure-dependent methods is changed accordingly and thus its numerical properties must be re-evaluated. In summary, a large stability factor can result in a large unconditional stability interval but also lead to a large relative period error. As a consequence, a stability factor must be appropriately chosen to have a desired unconditional stability interval in addition to an acceptable period distortion.

• Earthquakes and Structures
• /
• v.9 no.1
• /
• pp.173-194
• /
• 2015

Site response analysis is an important topic in earthquake engineering. A time-domain numerical method called as one-dimensional (1D) finite element artificial boundary method is proposed to simulate the homogeneous plane elastic wave propagation in a layered half space subjected to the obliquely incident plane body wave. In this method, an exact artificial boundary condition combining the absorbing boundary condition with the inputting boundary condition is developed to model the wave absorption and input effects of the truncated half space under layer system. The spatially two-dimensional (2D) problem consisting of the layer system with the artificial boundary condition is transformed equivalently into a 1D one along the vertical direction according to Snell's law. The resulting 1D problem is solved by the finite element method with a new explicit time integration algorithm. The 1D finite element artificial boundary method is verified by analyzing two engineering sites in time domain and by comparing with the frequency-domain transfer matrix method with fast Fourier transform.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.2
• /
• pp.87-95
• /
• 2018

This paper introduces a new explicit spectral element method for the simulation of SH-waves in semi-infinite domains. To simulate the wave motion in unbounded domains, it is necessary to reduce the infinite extent to a finite computational domain of interest. To prevent the wave reflection from the trunctated boundaries, perfectly matched layer(PML) wave-absorbing boundary is introduced. The forward problem for simulating SH-waves in PML-truncated domains can be formulated as second-order PDEs. The second-order semi-discrete form of the governing PDEs is constructed by using a mixed spectral elements with Legendre-gauss-Lobatto quadrature method, which results in a diagonalized mass matrix. Then the second-order semi-discrete form is transformed to a first-order, whose solutions are calculated by the fourth-order Runge-Kutta method. Numerical examples showed that solutions of SH-wave in the two-dimensional analysis domain resulted in stable and accurate, and reflections from truncated boundaries could be reduced by using PML boundaries. Elastic wave propagation analysis using explicit time integration method may be apt for solving larger domain problems such as three-dimensional elastic wave problem more efficiently.

• Structural Engineering and Mechanics
• /
• v.16 no.2
• /
• pp.153-176
• /
• 2003

From the equation of motion of a "bare" non-uniform beam (without any concentrated elements), an eigenfunction in term of four unknown integration constants can be obtained. When the last eigenfunction is substituted into the three compatible equations, one force-equilibrium equation, one governing equation for each attaching point of the concentrated element, and the boundary equations for the two ends of the beam, a matrix equation of the form [B]{C} = {0} is obtained. The solution of |B| = 0 (where ${\mid}{\cdot}{\mid}$ denotes a determinant) will give the "exact" natural frequencies of the "constrained" beam (carrying any number of point masses or/and concentrated springs) and the substitution of each corresponding values of {C} into the associated eigenfunction for each attaching point will determine the corresponding mode shapes. Since the order of [B] is 4n + 4, where n is the total number of point masses and concentrated springs, the "explicit" mathematical expression for the existing approach becomes lengthily intractable if n > 2. The "numerical assembly method"(NAM) introduced in this paper aims at improving the last drawback of the existing approach. The "exact"solutions in this paper refer to the numerical results obtained from the "continuum" models for the classical analytical approaches rather than from the "discretized" ones for the conventional finite element methods.

• Journal of computational fluids engineering
• /
• v.19 no.1
• /
• pp.27-33
• /
• 2014

A new and efficient implicit scheme is proposed to obtain a steady-state solution in time integration and the comparison of characteristics with the approximation ways for the implicit method to solve the incompressible Navier-Stokes equations is provided. The conservative, finite-volume cell-vertex upwind scheme and artificial compressibility method using dual time stepping for time accuracy is applied in this paper. The numerical results obtained indicate that the direct application of Jacobian matrix to the Lower and upper sweeps of implicit LU-SGS leads to better performance as well as convergence regardless of CFL number and true time step than explicit scheme and approximation of Jacobian matrix. The flow simulation around box in uniform flow with unstructured meshes is demonstrated to check the validity of the current formulation.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.4
• /
• pp.465-474
• /
• 2012

In this study, a dynamic simulation model that considers many spherical particles and multibody dynamics (MBD) entities is developed. Plenteous spherical particles are solved using the Discrete Element Method (DEM) technique and simulated on a GPU board in a PC. A fast algorithm is used to calculate the Hertzian contact forces between many spherical particles, and NVIDIA CUDA is used to increase the calculation speed. The explicit integration method is applied to solve the many spheres. MBD entities are simulated by recursive formulation. Constraints are reduced by recursive formulation, and the implicit generalized alpha method is applied to solve the dynamic model. A new algorithm is developed to simulate the DEM and MBD models simultaneously. As a numerical example, a truck car model and gear model are developed. The results show that the proposed algorithm using a general-purpose GPU in a PC has many advantages.

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

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.10 no.4
• /
• pp.141-150
• /
• 1990

Numerical experiments of sballow water equations are performed under various boundary conditions by finite element method to simulate the circulation in estuaries and coastal areas. Galerkin method is employed to discretize spatial domain, and for time integration, finite difference method (Crank-Nicolson scheme) is used. This method is tested in five problems, in which first four cases have analytic solutions. The computed values are well in agreement with the analytic solutions in four experiments and the result of the last 2-dimensional ease is resonable. Implicit and two step Lax-Wendroff schemes in time domain are compared, and the results when using four node bilinear and triangular elements are presented. Consequently it takes very long time for complex problems requiring many elements to integrate all the time steps using the implicit schemes. And the explicit scheme requires careful consideration in selecting the time step and the grid size to obtain the desired accuracy.