• 한국전산유체공학회지
• /
• 제9권3호
• /
• pp.57-65
• /
• 2004

The incompressible Navier-Stokes equations in two dimensions are stabilized by a modified residual method, and then discretized by hierarchical elements. The stabilization is necessary to escape from the Ladyzhenskaya-Babuska-Brezzi(LBB) constraint and hence to achieve an equal order formulation. To expedite a standard iterative method such as the conjugate gradient squared(CGS) method, a preconditioning technique called the Hierarchical Iterative Procedure(HIP) has been applied. In this paper, we increased the order of interpolation within an element up to cubic. The hierarchical elements have been used to achieve a higher order accuracy in fluid flow analyses, but a proper efficient iterative procedure for higher order finite element formulation has not been available so far The numerical results by the present HIP for the lid driven cavity flow and others showed the present procedure to be stable, very efficient and useful in flow analyses in conjunction with hierarchical elements.

• 대한수학회보
• /
• 제53권2호
• /
• pp.423-440
• /
• 2016

The bubble stabilization technique of Chebyshev-Legendre high-order element methods for one dimensional advection-diffusion equation is analyzed for the proposed scheme by Canuto and Puppo in [8]. We also analyze the finite element lower-order preconditioner for the proposed stabilized linear system. Further, the numerical results are provided to support the developed theories for the convergence and preconditioning.

• Steel and Composite Structures
• /
• 제47권3호
• /
• pp.365-374
• /
• 2023

This paper presents a new scheme for constructing locking-free finite elements in thick and thin plates, called Discrete Shear Gap element (DSG), using multiphase material topology optimization for triangular elements of Reissner-Mindlin plates. Besides, common methods are also presented in this article, such as quadrilateral element (Q4) and reduced integration method. Moreover, when the plate gets too thin, the transverse shear-locking problem arises. To avoid that phenomenon, the stabilized discrete shear gap technique is utilized in the DSG3 system stiffness matrix formulation. The accuracy and efficiency of DSG are demonstrated by the numerical examples, and many superior properties are presented, such as being a strong competitor to the common kind of Q4 elements in the static topology optimization and its computed results are confirmed against those derived from the three-node triangular element, and other existing solutions.

• Geomechanics and Engineering
• /
• 제9권4호
• /
• pp.415-426
• /
• 2015

Researches have been done to discover ways to strengthen peat soil deposits. In this model study, fibrous peat that is the most compressible types of peat has been reinforced with precast peat columns stabilized with ordinary Portland cement and polypropylene fibres. Rowe cell consolidation tests as well as plate load tests (PLTs) were conducted on various types of test samples to evaluate the strength and deformation of untreated peat and peat reinforced by various types of columns. PLTs were conducted in a specially designed and fabricated circular steel test tank. The compression index ($C_c$) and recompression index ($C_r$) of fibrous peat samples reduced considerably upon use of precast columns. Also, PLT results confirmed the results obtained from Rowe cell tests. Use of polypropylene fibres added to cement further decreased ($C_c$) and ($C_r$) and increased load bearing capacity of untreated peat. Finite element method (FEM) using Plaxis 3D was carried out to evaluate the stress distributions along various types of tested samples and also, to compare the deformations obtained from FEM analysis with the actual maximum deformations found from PLTs. FEM results indicate that most of the induced stresses are taken on the upper portion of tested samples and reach their maximum values below the loading plate. Also, a close agreement was found between actual deformation values obtained from PLTs and values resulted from FEM analysis for various types of tested samples.

• Geomechanics and Engineering
• /
• 제1권3호
• /
• pp.193-204
• /
• 2009

A 10.4-m high highway embankment retained behind mechanically stabilized earth (MSE) walls is under construction in the northeastern part of the Indian state of Bihar. The structure is constructed with compacted, micaceous, grey, silty sand, reinforced with polyester (PET) geogrids, and faced with reinforced cement concrete fascia panels. The connections between the fascia panels and the geogrids failed on several occasions during the monsoon seasons of 2007 and 2008 following episodes of heavy rainfall, when the embankment was still under construction. However, during these incidents the MSE embankment itself remained by and large stable and the collateral damages were minimal. The observational data during these incidents presented an opportunity to develop and calibrate a simple procedure for estimating rainfall induced pore water pressure development within MSE embankments constructed with backfill materials that do not allow unimpeded seepage. A simple analytical finite element model was developed for the purpose. The modeling results were found to agree with the observational and meteorological records from the site. These results also indicated that the threshold rainwater infiltration flux needed for the development of pore water pressure within an MSE embankment is a monotonically increasing function of the hydraulic conductivity of backfill. Specifically for the MSE embankment upon which this study is based, the analytical results indicated that the instabilities could have been avoided by having in place a chimney drain immediately behind the fascia panels.

• Geomechanics and Engineering
• /
• 제12권2호
• /
• pp.257-267
• /
• 2017

This paper presents laboratory investigation of stabilization of subgrade soil. One type of soil and three types of stabilizers i.e., hydrated lime, class F fly ash and polypropylene fibres are selected in the study. Atterberg limit, compaction, california bearing ratio (CBR), unconfined compressive strength and triaxial shear strength tests are conducted on unstabilized and stabilized soil for varying percentage of stabilizers to analyze the effect of stabilizers on the properties of soil. Vertical compressive strains at the top of unstabilized and stabilized subgrade soil were found out by elasto-plastic finite element analysis using commercial software ANSYS. Strategy for design of optimum pavement section was based on extension in service life (TBR) and reduction in layer thickness (LTR). Extension in service life of stabilized subgrade soil is 6.49, 4.37 and 3.26 times more due to lime, fly ash and fibre stabilization respectively. For a given service life of the pavement, there is considerable reduction in layer thicknesses due to stabilization. It helps in reduction in construction cost of pavement and saving in natural resources as well.

• Geomechanics and Engineering
• /
• 제36권3호
• /
• pp.277-293
• /
• 2024

Shored Mechanically Stabilized Earth (SMSE) walls are types of soil retaining structures that increase soil stability under static and dynamic loads. The damage caused by an earthquake can be determined by evaluating the probabilistic seismic response of SMSE walls. This study aimed to assess the seismic performance of SMSE walls and provide fragility curves for evaluating failure levels. The generated fragility curves can help to improve the seismic performance of these walls through assessing and controlling variables like backfill surface settlement, lateral deformation of facing, and permanent relocation of the wall. A parametric study was performed based on a non-linear elastoplastic constitutive model known as the hardening soil model with small-strain stiffness, HSsmall. The analyses were conducted using PLAXIS 2D, a Finite Element Method (FEM) program, under plane-strain conditions to study the effect of the number of geogrid layers and the axial stiffness of geogrids on the performance of SMSE walls. In this study, three areas of damage (minor, moderate, and severe) were observed and, in all cases, the wall has not completely entered the stage of destruction. For the base model (Model A), at the highest ground acceleration coefficient (1 g), in the moderate damage state, the fragility probability was 76%. These values were 62%, and 54%, respectively, by increasing the number of geogrids (Model B) and increasing the geogrid stiffness (Model C). Meanwhile, the fragility values were 99%, 98%, and 97%, respectively in the case of minor damage. Notably, the probability of complete destruction was zero percent in all models.

• 한국전산유체공학회지
• /
• 제12권1호
• /
• pp.35-42
• /
• 2007

This paper describes a recent development on the divergence free basis function based on a hermite stream function and verifies its validity by comparing results with those from a modified residual method known as one of stabilized finite element methods. It can be shown that a proper choice of degrees of freedom at a node with a proper arrangement of the hermite interpolation functions can yield solenoidal or divergent free interpolation functions for the velocities. The well-known cavity problem has been chosen for validity of the present algorithm. The comparisons from numerical results between the present and the modified residual showed the present method yields better results in both the velocity and the pressure within modest Reynolds numbers(Re = 1,000).

• 대한수학회지
• /
• 제50권5호
• /
• pp.1129-1163
• /
• 2013

Based on finite element discretization, two linearization approaches to the defect-correction method for the steady incompressible Navier-Stokes equations are discussed and investigated. By applying $m$ times of Newton and Picard iterations to solve an artificial viscosity stabilized nonlinear Navier-Stokes problem, respectively, and then correcting the solution by solving a linear problem, two linearized defect-correction algorithms are proposed and analyzed. Error estimates with respect to the mesh size $h$, the kinematic viscosity ${\nu}$, the stability factor ${\alpha}$ and the number of nonlinear iterations $m$ for the discrete solution are derived for the linearized one-step defect-correction algorithms. Efficient stopping criteria for the nonlinear iterations are derived. The influence of the linearizations on the accuracy of the approximate solutions are also investigated. Finally, numerical experiments on a problem with known analytical solution, the lid-driven cavity flow, and the flow over a backward-facing step are performed to verify the theoretical results and demonstrate the effectiveness of the proposed defect-correction algorithms.

• Steel and Composite Structures
• /
• 제21권2호
• /
• pp.429-442
• /
• 2016

The aim of this paper is to determine temperature and stress distributions in a ceramic based on Partially Stabilized Zirconia coated steel piston crown by using plasma spraying for improving performance of a marine diesel engine. Effects of coating constituent and thickness on temperature and stress distributions were investigated including comparisons with results from an uncoated piston by means of finite element method namely ANSYS. Temperature developed at the coated surface is significantly higher than that of the uncoated piston. The maximum stress components occur between bond coat and adjacent ceramic layer. Provided that coating thickness is constant as 0.5 mm, when numbers of layers increase, magnitude of the normal stress decrease about 34.1% on the base metal surface according to uncoated piston, but the base metal surface temperature of the steel piston increase about 13.1%.