• Title/Summary/Keyword: computational finite elements procedure

Search Result 38, Processing Time 0.028 seconds

A Coupling of Finite Elements and Boundary Elements for Half Plane Problems (반무한영역 문제에 대한 유한요소와 경계요소의 조합)

  • 김문겸;임윤묵
    • 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.

  • PDF

A Nonlinear Analysis of Half Plane Problems Using Coupling of Finite Elements and Boundary Elements (유한요소와 경계요소의 조합에 의한 반무한 영역 문제의 비선형해석)

  • 김문겸;임윤묵
    • 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.

  • PDF

Vibrations and stress analysis of perforated functionally graded rotating beams

  • Alaa A. Abdelrahman;Hanaa E. Abd-El-Mottaleb;Mohamed G. Elblassy;Eman A. Elshamy
    • 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.

Adaptive finite elements by Delaunay triangulation for fracture analysis of cracks

  • Dechaumphai, Pramote;Phongthanapanich, Sutthisak;Bhandhubanyong, Paritud
    • 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.

An Adaptive Procedure in Finite Element Analysis of Elastodynamic Problems (적응적 방법을 이용한 동적 유한요소해석)

  • 최창근;정흥진
    • 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.

  • PDF

An Incompressible Flow Computation by a Hierarchical Iterative and a Modified Residual Method (계층적 반복과 수정 잔여치법에 의한 비압축성 유동 계산)

  • Kim J. W.
    • 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.

An Incompressible Flow Computation by a Hierarchical Iterative Preconditioning (계층적 반복의 예조건화에 의한 비압축성 유동 계산)

  • Kim J. W.;Jeong C. R.
    • 한국전산유체공학회:학술대회논문집
    • /
    • 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.

  • PDF

An edge-based smoothed finite element method for adaptive analysis

  • Chen, L.;Zhang, J.;Zeng, K.Y.;Jiao, P.G.
    • 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.

Free Vibration Analysis of Horizontally Curved Multi-Girder Bridges (수평곡선 격자형교의 자유진동해석)

  • 윤기용;강영종
    • 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.

  • PDF

A Study on the Determination of Contact Area of a Plate on Elastic Half-Space (탄성지반 위에 놓인 평판의 접촉영역 결정에 관한 연구)

  • 정진환;이외득;김동석
    • 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.

  • PDF