• Journal of Wetlands Research
• /
• v.13 no.3
• /
• pp.395-409
• /
• 2011

The objective of this study is to develop an accurate and robust two-dimensional finite element method for turbulence simulation in open channels. The model is based on Streamline Upwind/Petrov-Galerkin finite element method and Boussinesq's eddy viscosity theory. The method developed in the study is depth-averaged mixing length model which assumes anisotropic and local equilibrium state of turbulence. The model calibration and validation were performed by comparing with analytical solutions and observed data. Several numerical simulations were carried out, which examined the performance of the turbulence model for the purpose of sensitivity analysis. The uniform channels that appear horizontal flow and vertical flow were carried out. The model was also applied to the Han river was in for the applicability test. The results were compared with the observed data. The suggested model displayed reasonable flow distribution compare to the observed data in natural river flow. As a result of this study, the two-dimensional finite element model provides a reliable results for flow distribution based on the turbulence simulation in open channels.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.6 no.1
• /
• pp.61-68
• /
• 1986

A set of the shape functions which perfectly satisfy the homogeneous Euler Equations has been proposed for deep beam problems. A finite element beam model using the proposed shape functions has been derived by the Galerkin weighted residual method and used to analyze the numerical examples without reduced shear integration, to show the accuracy and efficiency of the proposed shape functions. The result shows that the finite element model using the proposed shape functions gives very accurate solutions for both static and free vibration analyses. The concept of the proposed shape functions is thought to be applied for the finite element analysis of the elasto-static problems.

• Structural Engineering and Mechanics
• /
• v.85 no.5
• /
• pp.593-605
• /
• 2023

This work presents a comparison between analytical and finite element analysis for bending of porous sandwich functionally graded material (FGM) plates. The plate is rectangular and simply supported under static sinusoidal loading. Material properties of FGM are assumed to vary continuously across the face sheets thickness according to a power-law function in terms of the volume fractions of the constituents while the core is homogeneous. Four types of porosity are considered. A refined higher-order shear with normal deformation theory is used. The number of unknowns in this theory is five, as against six or more in other shear and normal deformation theories. This theory assumes the nonlinear variation of transverse shear stresses and satisfies its nullity in the top and bottom surfaces of the plate without the use of a shear correction factor. The governing equations of equilibrium are derived from the virtual work principle. The Navier approach is used to solve equilibrium equations. The constitutive law of the porous FGM sandwich plate is implemented for a 3D finite element through a subroutine in FORTRAN (UMAT) in Abaqus software. Results show good agreement between the finite element model and the analytical method for some results, but the analytical method keeps giving symmetric results even with the thickness stretching effect and load applied to the top surface of the sandwich.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.2 no.1
• /
• pp.111-129
• /
• 1998

An integral equation representation of cracks was presented, which differs from well-known "dislocation layer" representation. In this new representation, an integral equation representation of cracks was developed and coupled to the direct boundary-element method for treatment of cracks in plane finite bodies. The method was developed for in-plane(modes I and II) loadings only. In this paper, the method is formulated and applied to mode III problems involving smooth or kinked cracks in finite region. The results are compared to exact solutions where available and the method is shown to be very accurate despite of its simplicity.

• Journal of the Korean Mathematical Society
• /
• v.59 no.3
• /
• pp.519-548
• /
• 2022

In this paper, a stabilized-penalized collocated finite volume (SPCFV) scheme is developed and studied for the stationary generalized Navier-Stokes equations with mixed Dirichlet-traction boundary conditions modelling an incompressible biological fluid flow. This method is based on the lowest order approximation (piecewise constants) for both velocity and pressure unknowns. The stabilization-penalization is performed by adding discrete pressure terms to the approximate formulation. These simultaneously involve discrete jump pressures through the interior volume-boundaries and discrete pressures of volumes on the domain boundary. Stability, existence and uniqueness of discrete solutions are established. Moreover, a convergence analysis of the nonlinear solver is also provided. Numerical results from model tests are performed to demonstrate the stability, optimal convergence in the usual L² and discrete H¹ norms as well as robustness of the proposed scheme with respect to the choice of the given traction vector.

• Structural Engineering and Mechanics
• /
• v.8 no.2
• /
• pp.207-231
• /
• 1999

A versatile 4-node shell element which is useful for the analysis of arbitrary shell structures is presented. The element is developed by flat shell approach, i.e., by combining a membrane element with a Mindlin plate element. The proposed element has six degrees of freedom per node and permits an easy connection to other types of finite elements. In the plate bending part, an improved Mindlin plate has been established by the combined use of the addition of non-conforming displacement modes (N) and the substitute shear strain fields (S). In the membrane part, the nonconforming displacement modes are also added to the displacement fields to improve the behavior of membrane element with drilling degrees of freedom and the modified numerical integration (M) is used to overcome the membrane locking problem. Thus the element is designated as NMS-4F. The rigid link correction technique is adopted to consider the effect of out-of-plane warping. The shell element proposed herein passes the patch tests, does not show any spurious mechanism and does not produce shear and membrane locking phenomena. It is shown that the element produces reliable solutions even for the distorted meshes through the analysis of benchmark problems.

• Proceeding of KASS Symposium
• /
• 2008.05a
• /
• pp.95-100
• /
• 2008

A 4-node assumed strain finite element based on higher order shear deformation theory is developed to investigate the behaviours of symmetric and unsymmetric laminated composite plates. The present element is based on Reddy's higher order shear deformation theory so that it can consider the parabolic distribution of shear deformation through plate thickness direction. In particular, assumed strain method is adopted to alleviate the shear locking phenomena inherited plate elements based on higher order shear deformation theory. The present finite element has seven degrees of freedom per node and denoted as HSA4. Numerical examples are carried out for symmetric and unsymmetric laminated composite plate with various thickness values. Numerical results are compared with reference solutions produced by other higher order shear deformation theories.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.6
• /
• pp.1592-1602
• /
• 1990

Mixed mode surface crack in finite-width plate subjected to uniform tension and bending has been analyzed in 3-D problem by using boundary element method. The calculations were carried out for the surface crack angles(.a/pha.) of 0.deg., 15.deg., 30.deg., 45.deg., 60.deg., and 75.deg., and for the aspect ratio(a/c) of 0.2, 0.4, 0.6 and 1.0 to get stress intensity factors at the boundary points of the surface crack. For the aspect ratio of 1.0 and the surface crack angles, finite element method was used to check the results in this study. Comparison of the results from both methods showed good agreement.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1994.04a
• /
• pp.145-152
• /
• 1994

The dynamic interaction between tunnel structures and their surrounding soil medium due to impulse loading is investigated by a hybrid IEM/FEM methodology. A dynamic infinite element is developed for the efficient descretization of the far-field region of the unbounded soil medium. The shape functions of the infinite element are constructed based on the far-field solutions which are obtained by solving the 2-D elastic wave problems. Also they are devised to obtain a reasonable result over all frequency range. Numerical analysis is carried out to examine the response of the tunnel subjected to simple rectangular impulse. It is indicated that the results by the present method are in good accord with those by the boundary and finite element coupling method.

• Journal of Mechanical Science and Technology
• /
• v.17 no.1
• /
• pp.64-76
• /
• 2003

In this work, a new penalty formulation is proposed for the analysis of Mindlin-Reissner plates by using the element-free Galerkin method. A penalized weak form for the Mindlin-Reissner Plates is constructed through the exterior penalty method to enforce the essential boundary conditions of rotations as well as transverse displacements. In the numerical examples, some typical problems of Mindlin-Reissner plates are analyzed, and parametric studies on the order of integration and the size of influence domain are also carried out. The effect of the types of background cells on the accuracy of numerical solutions is observed and a proper type of background cell for obtaining optimal accuracy is suggested. Further, optimal order of integration and basis order of Moving Least Squares approximation are suggested to efficiently handle the irregularly distributed nodes through the triangular type of background cells. From the numerical tests, it is identified that unlike the finite element method, the proposed element-free Galerkin method with penalty technique gives highly accurate solution without shear locking in dealing with Mindlin-Reissner plates.