• 제목/요약/키워드: Numerical integration of discontinuity

검색결과 11건 처리시간 0.022초

불규칙한 다층 물성을 가지는 지반의 지진 응답 해석 (Earthquake Response Analysis considering Irregular Soil Layers)

  • 박장호;박재균
    • 한국지진공학회논문집
    • /
    • 제9권6호
    • /
    • pp.67-73
    • /
    • 2005
  • 엄밀한 지반-구조물 상호작용 해석을 위해서는 지반에 대한 묘사가 적절히 이루어져야 한다. 그러나 물성이 불연속적이거나 지층의 모양이 복잡할 때에는 유한요소망을 구성하는 것이 어려운데, 그 이유는 유한요소망이 지층의 모양을 정확히 나타내어야하기 때문이다. 이러한 어려움을 극복하기 위하여 본 연구에서는 지층의 불연속선을 따라 요소를 재배열하지 않고 엄밀하게 수치적분이 가능한 방법을 채용하였다. 그 결과로 정렬된 요소망을 그대로 이용하게 되어 강성행렬의 성질이 좋아지며, 해석 또한 신뢰성이 높아지게 되었다. 이 해석 방법을 이용하여 지진응답에 미지는 지층의 영향을 조사하였으며, 결과적으로 복잡하고 불연속적인 물성을 가지는 지층으로 이루어진 지반과 구조물의 상호작용을 해석하는 데에 쉽게 이용 될 수 있음을 보였다.

복잡한 지반층을 고려한 지반-말뚝-구조물의 상호작용 동해석 (Dynamic Analysis of Soil-Pile-Structure Interaction Considering a Complex Soil Profile)

  • 박장호;박재균
    • 한국지진공학회논문집
    • /
    • 제13권3호
    • /
    • pp.21-28
    • /
    • 2009
  • 지반-말뚝-구조의 상호작용을 정밀하게 해석하기 위해서는 토층, 말뚝 그리고 구조물의 적절한 묘사가 필요하다. 일반적으로 사용하는 유한요소해석의 경우에는 지반이나 구조물의 물성이 바뀌는 경계를 따라서 요소의 경계가 정해지게 된다. 그러나 실제로는 토층 단면과 말뚝의 형상이 매우 복잡하여 요소의 배열이 매우 어려운 작업이 될 수 있다. 이 어려움을 해결하기 위하여, 이 논문에서는 불연속선의 위치에 관계없이 규칙적인 요소를 사용하여 쉽게 적분을 가능하게 하는 다른 적분 방법을 채택하였다. 이 방법을 적용함으로 써 요소는 매우 빠르고 규칙적인 강성 매트릭스를 만든다. 구조물 응답에 대한 토층과 말뚝의 영향을 조사하였고, 예를 통하여 본 방법의 유효성을 보였다. 탄성 말뚝의 사용으로 20% 대의 가속도 감소 효과를 얻었고 지반 층의 모양에 따라 그 영향이 변하는 것을 확인하였다.

콘크리트의 변형률 국소화 및 진행성 파괴에 관한 연구 (Study on Strain Localization and Progressive Failure of Concrete)

  • 송하원;김형운;우승민
    • 콘크리트학회논문집
    • /
    • 제11권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.

충격해석을 위한 새로운 불연속 시간적분법의 개발 (Development of a New Discontinuous Time Integration Method for Transient Analysis of Impact Phenomena)

  • 조진연;김승조
    • 한국소음진동공학회:학술대회논문집
    • /
    • 한국소음진동공학회 1998년도 춘계학술대회논문집; 용평리조트 타워콘도, 21-22 May 1998
    • /
    • 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.

  • PDF

탄성균열해석을 위한 그리드 없는 유한차분법 (A Gridless Finite Difference Method for Elastic Crack Analysis)

  • 윤영철;김동조;이상호
    • 한국전산구조공학회논문집
    • /
    • 제20권3호
    • /
    • pp.321-327
    • /
    • 2007
  • 본 연구는 탄성균열문제를 신속하고 정확하게 해석할 수 있는 새로운 개념의 그리드(grid) 없는 유한차분법을 제시한다. 이동최소제곱법을 이용한 Taylor 전개식 구성을 통해 직접적인 미분계산 없이 근사함수와 그 미분을 손쉽게 계산한다. 그리드로 인한 절점 간의 종속성이 없어 해석영역 내의 불연속면 모델링이 용이하여 차분식 구성시 균열로 인한 불연속 효과를 고려하는 과정도 자연스럽다. 유한차분법에 근간을 두고 있어 지배 미분방정식을 직접 이산화하기 때문에 수치적분이 필요한 수치기법에 비해 계산속도도 빠르다. 모드 I과 모드 II 균열문제 해석을 통해 본 해석기법이 정확하고 효율적으로 응력확대계수를 계산할 수 있음을 보였다.

정자력 계산을 위한 적응 유한 요소법 (An Adaptive Finite Element Method for Magnetostatic Force Computations)

  • 박용규;박일한;정형석;정현교;이기식;한송엽
    • 대한전기학회논문지
    • /
    • 제38권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.

  • PDF

요소 세분화 및 재결합을 이용한 바람의 적응적 유한요소 해석 (Adaptive finite element wind analysis with mesh refinement and recovery)

  • 최창근;유원진;이은진
    • 한국전산구조공학회:학술대회논문집
    • /
    • 한국전산구조공학회 1998년도 봄 학술발표회 논문집
    • /
    • 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.

  • PDF

Combination of isogeometric analysis and extended finite element in linear crack analysis

  • Shojaee, S.;Ghelichi, M.;Izadpanah, E.
    • Structural Engineering and Mechanics
    • /
    • 제48권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.

정자력 계산을 위한 적응 유한 요소법 (An Adaptive Finite Element Method for Magnetostatic Force Computations)

  • 박용규;정현교;이기식;한송엽
    • 대한전기학회:학술대회논문집
    • /
    • 대한전기학회 1988년도 추계학술대회 논문집 학회본부
    • /
    • 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.

  • PDF

Adaptive finite element wind analysis with mesh refinement and recovery

  • Choi, Chang-Koon;Yu, Won-Jin
    • Wind and Structures
    • /
    • 제1권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.