• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1988.10a
• /
• pp.49-54
• /
• 1988

A procedure which may be useful in dealing with problems of half plane is considered. Boundary elements are combined with finite elements to facilitate their merits. Boundary elements for semi-infinite region are composed using the Melan's solution for half plane. Finite elements are used to model irregurarity or the nonhomogeneity of materials, which is usual in underground structures. In order to verify the procedure, a shallow tunnel under internal pressure is analysed using the finite element method, the boundary element method, and the combined method. It is shown that the developed Procedure is accurate enough compared with other method.

• Computational Structural Engineering
• /
• v.2 no.1
• /
• pp.55-64
• /
• 1989

A procedure which may be useful in dealing with problems of half plane is considered. Boundary elements are combined with nonlinear finite elements to facilitate their merits. Boundary elements for semi-infinite region are composed using the Melan's solution for half plane. Nonlinear finite elements are used to model irregularity or nonhomogeneity of elasto-plastic materials, which is usual in underground structures. In order to verify the procedure, a shallow tunnel under internal pressure is analysed using the nonlinear finite element method and combined method. It is shown that the developed procedure is accurate enough compared with other method.

• Steel and Composite Structures
• /
• v.49 no.6
• /
• pp.667-684
• /
• 2023

In the context of finite element method, a computational simulation is presented to study and analyze the dynamic behavior of regularly perforated functionally graded rotating beam for the first time. To investigate the effect of perforation configurations, both regular circular and squared perforation patterns are studied. To explore impacts of graded material distributions, both axial and transverse gradation profiles are considered. The material characteristics of graded materials are assumed to be smoothly and continuously varied through the axial or the thickness direction according the nonlinear power gradation law. A computational finite elements procedure is presented. The accuracy of the numerical procedure is verified and compared. Resonant frequencies, axial displacements as well as internal stress distributions throughout the perforated graded rotating cantilever beam are studied. Effects of material distributions, perforation patterns, as well as the rotating beam speed are investigated. Obtained results proved that the graded material distribution has remarkable effects on the dynamic performance. Additionally, circular perforation pattern produces more softening effect compared with squared perforation configuration thus larger values of axial displacements and maximum principal stresses are detected. Moreover, squared perforation provides smaller values of nondimensional frequency parameters at most of vibration modes compared with circular pattern.

• Structural Engineering and Mechanics
• /
• v.15 no.5
• /
• pp.563-578
• /
• 2003

Delaunay triangulation is combined with an adaptive finite element method for analysis of two-dimensional crack propagation problems. The content includes detailed descriptions of the proposed procedure which consists of the Delaunay triangulation algorithm and an adaptive remeshing technique. The adaptive remeshing technique generates small elements around the crack tips and large elements in the other regions. Three examples for predicting the stress intensity factors of a center cracked plate, a compact tension specimen, a single edge cracked plate under mixed-mode loading, and an example for simulating crack growth behavior in a single edge cracked plate with holes, are used to evaluate the effectiveness of the procedure. These examples demonstrate that the proposed procedure can improve solution accuracy as well as reduce total number of unknowns and computational time.

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

An automated procedure which allows adaptation of spatial and time discretization simultaneously in finite element analysis of linear elastodynamic problems is presented. For dynamic problems having responses dominated by high frequency modes, such as those with impact, explosive, traveling and earthquake loads high gradient stress regions change their locations from time to time. And the time step size may need to vary in order to deal with whole process ranging from transient phase to steady state phase. As the sizes of elements in space vary in different regions, the procedure also permits different time stepping. In such a way, the best performance attainable by the finite element method can be achieved. In this study, we estimate both of the kinetic energy error and stran energy error induced by spatial and time discretization in a consistent manner. Numerical examples are used to demonstrate the performance of the procedure.

• Journal of computational fluids engineering
• /
• v.9 no.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.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.91-98
• /
• 2004

In two dimensional incompressible flows, a preconditioning technique called Hierarchical Iterative Procedure(HIP) has been implemented on a stabilized finite element formulation. The stabilization has been peformed by a modified residual method proposed by Illinca et. al.[3]. The stabilization which is necessary to escape from the LBB constraint renders an equal order formulation. In this paper, we increased the order of interpolation whithin an element up to cubic. The conjugate gradient squared(CGS) method is used for the outer iteration, and the HIP for the preconditioning for the incompressible Navier-Stokes equation. The hierarchical elements has 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 showed the present procedure to be stable, very efficient and useful in flow analyses in conjunction with hierarchical elements.

• Structural Engineering and Mechanics
• /
• v.39 no.6
• /
• pp.767-793
• /
• 2011

An efficient edge-based smoothed finite element method (ES-FEM) has been recently developed for solving solid mechanics problems. The ES-FEM uses triangular elements that can be generated easily for complicated domains. In this paper, the complexity study of the ES-FEM based on triangular elements is conducted in detail, which confirms the ES-FEM produces higher computational efficiency compared to the FEM. Therefore, the ES-FEM offers an excellent platform for adaptive analysis, and this paper presents an efficient adaptive procedure based on the ES-FEM. A smoothing domain based energy (SDE) error estimate is first devised making use of the features of the ES-FEM. The present error estimate differs from the conventional approaches and evaluates error based on smoothing domains used in the ES-FEM. A local refinement technique based on the Delaunay algorithm is then implemented to achieve high efficiency in the mesh refinement. In this refinement technique, each node is assigned a scaling factor to control the local nodal density, and refinement of the neighborhood of a node is accomplished simply by adjusting its scaling factor. Intensive numerical studies, including an actual engineering problem of an automobile part, show that the proposed adaptive procedure is effective and efficient in producing solutions of desired accuracy.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1996.10a
• /
• pp.55-61
• /
• 1996

In the present study, a numerical formulation procedure fer free vibration analysis of thin-walled horizontally curved multi-girder bridges is presented. The presented finite element procedure consists of curved and straight beam elements including warping degree of freedom. The homogeneous solutions of curved beam equations were used for shape functions in numerical formulation to achieve good convergence. In the straight beam element, the third order hermite polynomials were used fer shape functions. The Gupta method was used to solve the eigenvalue problem efficiently. The developed numerical procedure was applied to investigate the characteristics of free vibration of horizontally curved multi-girder bridges with varing subtended angle.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.10a
• /
• pp.405-412
• /
• 1998

According to the relative stiffness between the half-space and plate or loading condition, some parts of the plate can be separated from the half-space. The finite element procedure to determine the contact area by considering the distribution of contact pressure between plate and the elastic half-space is developed. The vertical surface displacements of the elastic half-space can be obtained through the integrations of the Boussinesq's solution for a point load. The rectangular plate on the elastic half-space is modeled by the 8-node rectangular and 6-node triangular elements and the Mindlin plate theory is used in oder to consider the transverse shear effect. In this study, the contact area may be determined approximately by the analysis with rectangular elements. From this results, the mesh pattern is modified by using triangular and rectangular elements. The contact area can be determined by the new mesh pattern with a relatively sufficient accuracy.