• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1988년도 가을 학술발표회 논문집
• /
• 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.

• 전산구조공학
• /
• 제2권1호
• /
• pp.55-64
• /
• 1989

본 논문에서는 경계요소법과 비선형 유한요소법의 각 장점을 이용하여 반무한 영역을 가진 구조체의 해석방법을 논하였다. 여기서, 반무한 경계요소는 Melan의 반무한 평면에 대한 해로부터 구성하였다. 비선형 유한요소는 지하구조물에서 주로 접할 수 있는 탄소성 재료의 비균질성 또는 불규칙성을 모형화하기 위하여 사용하였다. 본 조합방법의 검증을 위하여 얕은 터널에 일정한 내압이 작용하는 경우를 택하여, 비선형 유한요소법과 조합방법의 결과를 비교하였다. 비교결과, 개발된 조합방법이 다른 해석방법에 비해 충분한 정확도를 가짐을 알 수 있었다.

• Steel and Composite Structures
• /
• 제49권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
• /
• 제15권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.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1994년도 봄 학술발표회 논문집
• /
• 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.

• 한국전산유체공학회지
• /
• 제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.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2004년도 춘계 학술대회논문집
• /
• 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
• /
• 제39권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.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1996년도 가을 학술발표회 논문집
• /
• 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.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1998년도 가을 학술발표회 논문집
• /
• 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.