• Journal of the Earthquake Engineering Society of Korea
• /
• v.9 no.6 s.46
• /
• pp.67-73
• /
• 2005

Precise analysis of soil-structure interaction requires a proper description of soil profile. However, such approach becomes generally nearly unpractical for soil exhibiting material discontinuity and complex geometry since meshes should match that material discontinuity line. To overcome these difficulties, a different numerical integration method is adopted in this paper, which enables to integrate easily over an element with material discontinuity without regenerating mesh fellowing the discontinuity line. As a result the mesh is highly structured, loading to very regular silliness matrix. The influence of the shape of soil profile on the response is examined and it is seen that the proposed soil-structure analysis method can be easily used on soil-structure interaction problems with complicated soil profile and produce reliable results regardless of material discontinuities.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.13 no.3
• /
• pp.21-28
• /
• 2009

The precise analysis of soil-pile-structure interaction requires a proper description of soil layer, pile, and structure. In commonly used finite element simulations, mesh boundaries should match the material discontinuity line. However, in practice, the geometry of soil profiles and piles may be so complex that mesh alignment becomes a wasteful and difficult task. To overcome these difficulties, a different integration method is adopted in this paper, which enables easy integration over a regular element with material discontinuity regardless of the location of the discontinuity line. By applying this integration method, the mesh can be generated rapidly and in a highly structured manner, leading to a very regular stiffness matrix. The influence of the shape of the soil profile and piles on the response is examined, and the validity of the proposed soil-pile structure interaction analysis method is demonstrated through several examples. It is seen that the proposed analysis method can be easily used on soil-pile-structure interaction problems with complex interfaces between materials to produce reliable results regardless of the material discontinuity line.

• Journal of the Korea Concrete Institute
• /
• v.11 no.3
• /
• pp.181-192
• /
• 1999

The progressive failure following strain localization in concrete can be analyzed effectively using finite element modeling of fracture process zone of concrete with a finite element embedded discontinuity. In this study, a finite element with embedded discontinuous line is utilized for the analysis of progressive failure in concrete. The finite element with embedded discontinuity is a kind of discrete crack element, but the difficulties in discrete crack approach such as remeshing or adding new nodes along with crack growth can be avoided. Using a discontinuous shape function for this element, the displacement discontinuity is embedded within an element and its constitutive equation is modeled from the modeling of fracture process zone. The element stiffness matrix is derived and its dual mapping technique for numerical integration is employed. Then, a finite element analysis program with employed algorithms is developed and failure analysis results using developed finite element program are verified through the comparison with experimental data and other analysis results.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1998.04a
• /
• pp.407-412
• /
• 1998

In this work, a new time integration method is proposed using the generalized derivative concept to simulate the dynamic phenomena having sudden constraint occurring in dynamic contact/impact problems. By the adoption of the generalized derivative concept and jump assumption, discontinuity can be incorporated in time integration and as a result, the algorithm does not need any other special consideration of jumps in dynamic field variables due to sudden constraint like dynamic contact-release conditions. To observe the characteristics of the proposed time integration method, the stability and convergence analyses are carried out. In numerical tests, several dynamic contact/impact problems are analyzed by straightforward application of the proposed time integration method with the exterior penalty method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.3
• /
• pp.321-327
• /
• 2007

This study presents a new gridless finite difference method for solving elastic crack problems. The method constructs the Taylor expansion based on the MLS(Moving Least Squares) method and effectively calculates the approximation and its derivatives without differentiation process. Since no connectivity between nodes is required, the modeling of discontinuity embedded in the domain is very convenient and discontinuity effect due to crack is naturally implemented in the construction of difference equations. Direct discretization of the governing partial differential equations makes solution process faster than other numerical schemes using numerical integration. Numerical results for mode I and II crack problems demonstrates that the proposed method accurately and efficiently evaluates the stress intensity factors.

• The Transactions of the Korean Institute of Electrical Engineers
• /
• v.38 no.2
• /
• pp.100-105
• /
• 1989

This paper presents an adaptive finite element method for magnetostatic force computation using Maxwell's stress tensor. Mesh refinements are performed automatically by interelement magnetic field intensity discontinuity errors and element force errors. In initial mesh, the computed forces for different integration paths give great differences, but converge to a certain value as mesh division is performed by the adaptive scheme, We obtained good agreement between analytic solutions and numerical values in typical examples.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.04a
• /
• pp.60-67
• /
• 1998

This paper deals with the development of a variable-node element and its application to the adaptive h-version mesh refinement-recovery for the incompressible viscous flow analysis. The element which has variable mid-side nodes can be used in generating the transition zone between the refined and unrefined elements and efficiently used for construction of a refined mesh without generating distorted elements. A modified Gaussian quadrature is needed to evaluate the element matrices due to the discontinuity of derivatives of the shape functions used for the element. The penalty function method which can reduce the number of independent variables is adopted for the purpose of computational efficiency and the selective reduced integration is carried out for the convection and pressure terms to preserve the stability of solution. For the economical analysis of transient problems, not only the mesh refinement but also the mesh recovery is needed. The numerical examples show that the optimal mesh for the finite element analysis of a wind around the structures can be obtained automatically by the proposed scheme.

• Structural Engineering and Mechanics
• /
• v.48 no.1
• /
• pp.125-150
• /
• 2013

This paper intends to present an application of isogeometric analysis in crack problems. An isogeometric formula is developed based on NURBS basis functions - enriched and adopted via X-FEM enrichment functions. The proposed method which is represented by the combination of the two above-mentioned methods, first by using NURBS functions models the geometry exactly and then by defining level set function on domain, identifies available discontinuity in elements. Additional DOFs are allocated to elements containing the crack and X-FEM enrichment functions enrich approximate solution. Moreover, a subelement refinement technique is used to improve the accuracy of integration by the Gauss quadrature rule. Finally, several numerical examples are illustrated to demonstrate the effectiveness, robustness and accuracy of the proposed method during calculation of crack parameters.

• Proceedings of the KIEE Conference
• /
• 1988.11a
• /
• pp.24-27
• /
• 1988

This paper presents an adaptive finite element method for magnetostatic force computation using Maxwell's stress tensor Mesh refinements are performed automatically by interelement flux density discontinuity errors and element force errors. In initial mesh, the computed forces for different Integration paths give great differences. but converge to a certain value as mesh division is performed by the adaptive scheme, We obtained good agreement between analytic solutions and numerical values In typical examples.

• Wind and Structures
• /
• v.1 no.1
• /
• pp.111-125
• /
• 1998

This paper deals with the development of variable-node element and its application to the adaptive h-version mesh refinement-recovery for the incompressible viscous flow analysis. The element which has variable mid-side nodes can be used in generating the transition zone between the refined and unrefined element and efficiently used for the construction of a refined mesh without generating distorted elements. A modified Guassian quadrature is needed to evaluate the element matrices due to the discontinuity of derivatives of the shape functions used for the element. The penalty function method which can reduce the number of the independent variables is adopted for the purpose of computational efficiency and the selective reduced integration is carried out for the convection and pressure terms to preserve the stability of solution. For the economical analysis of transient problems in which the locations to be refined are changed in accordance with the dynamic distribution of velocity gradient, not only the mesh refinement but also the mesh recovery is needed. The numerical examples show that the optimal mesh for the finite element analysis of a wind around the structures can be obtained automatically by the proposed scheme.