• Computational Structural Engineering
• /
• v.3 no.3
• /
• pp.53-54
• /
• 1990

응력 변형율의 관계가 시간에 대한 미분의 형태로 나타나는 비선형 탄소성 혹은 점탄소성 재질을 갖는 구조물이나 지만의 거동 문제를 유한요소법 등의 방법을 이용하여 해결하려고 하는 경우 주어진 외력에 의한 새로운 응력이나 응력 강화 현상을 표현하는 여러 재료 상수값들을 구하기 위해서는 적분을 요하게 되며 일반적으로 수치해석적 방법에 의해 수행된다. 이러한 수치해석적 적분방법은 보다 정확한 결과를 얻기 위하여 알고리즘 자체의 정확성과 안정성이 요구된다. 정확성은 수치해석적 적분방법이 적용될 수 있는 step size에 관계없이 거의 동일한 결과치를 얻을 수 있느냐 하는 것을 말하고 안정성은 큰 step size에서도 수렴된 결과치를 얻을 수 있느냐 하는 것을 의미한다. 그 뿐만 아니라 비교적 복잡하고도 그 대상영역이 큰 문제를 해석하고자 할 때는 수렴속도 또한 빠른 해석방법이 바람직하게 된다. 따라서 본 기사에서는 여러가지 가능한 수치해석 적분 방법을 소개하고 그들의 장단점을 논하고자 한다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.2
• /
• pp.129-140
• /
• 1999

추계론적 해석은 구조계 내의 해석인수에 존재하는 공간적 또는 시간적 임의성이 구조계 반응에 미치는 영향에 대한 고찰을 목적으로 한다. 확률장은 구족계 내에서 특정한 확률분포를 가지는 것으로 가정된다. 구조계 반응에 대한 이들 확률장의 영향 평가를 위하여 통계학적 추계론적 해석과 비통계학적 추계론적 해석이 사용되고 있다. 본 연구에서는 비통계학적 추계론적 해석방법 중의 하나인 가중적분법을 제안하였다. 특히 구조계의 공간적 임의성이 큰 특성을 가지고 있는 반무한영역에 대한 적용 예를 제시하고자 한다. 반무한영역의 모델링에는 무한요소를 사용하였다. 제안된 방법에 의한 해석 결과는 통계학적 방법인 몬테카를로 방법에 의한 결과와 비교되었다. 제안된 가중적분법은 자기상관함수를 사용하여 확률장을 고려하므로 무한영역의 고려에 따른 해석의 모호성을 제거할 수 있다. 제안방법과 몬테카를로 방법에 의한 결과는 상호 잘 일치하였으며 공분산 및 표준편차는 무한요소의 적용에 의하여 매우 개선된 결과를 나타내었다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.21 no.4
• /
• pp.375-382
• /
• 2008

The particular integral formulation for two(2D) and three(3D) dimensional inelastic transient dynamic stress analysis is presented. The elastostatic equation is used for the complementary solution. Using the concept of global shape function, the particular integrals for displacement and traction rates are obtained to approximate acceleration of the inhomogeneous equation. The Houbolt time integration scheme is used for the time-marching process. The Newton-Raphson algorithm for plastic multiplier is used to solve the system equation. Numerical results of four example problems are given to demonstrate the validity and accuracy of the present formulation.

• Proceedings of the Korean Society of Coastal and Ocean Engineers Conference
• /
• 1991.07a
• /
• pp.18-24
• /
• 1991

항내 정온도 해석은 일반적으로 유한차분법, 유한요소법 및 경계적분 방정식법 등의 엄밀해법과 근사 경계적분법, 고산의 방법 및 파향선법 등의 근사해법으로 구분된다. 엄밀해법은 지배방정식을 이산화 이외의 근사를 사용하지 않고 푸는 수치계산 방법으로 임의형상에의 적용성과 엄밀성이 뛰어나나 대상으로 하는 파의 파장이 짧고 항의 규모가 큰 경우에는 계산용량이 증대되여 실용적이지 못하다.(중략)

• Journal of the Korean Magnetics Society
• /
• v.1 no.2
• /
• pp.79-84
• /
• 1991

A variational approach employing localized functional is presented to solve alternating magnetic field problems with open boundary. The functional used in the approach consists of the domain integral of finite element region only and the boundary integral of the interfacial boundary between the finite and infinite element regions. The boundary integral is obtained by transforming the infinite domain integral for the infinite element region into the interfacial boundary integral. The proposed algorithm is then applied to a simple two-dimensional problem where the analytic solutions are available. It is shown that the algorithm makes it possible to yield good agreements between the numerical and analytic solutions. and that it requires less computer storage memory and computation time than the conventional finite element method due to the reduction of the computing region.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2009.04a
• /
• pp.63-66
• /
• 2009

본 논문에서는 트림영역이 있는 NURBS 평면을 등기하 해석할 수 있는 방법을 제시한다. 기존 등기하 해석법으로 트림 NURBS 곡면을 해석하기 위해서는 해석 도메인이 여러 개의 사각형 패치로 분할되어있어야 한다. 그러나 본 연구에서 제안한 방법은 CAD에서 제공하는 트림곡선의 정보를 해석에 직접 사용할 수 있기 때문에 CAD 모델을 별도로 재구성해야하는 번거로움이 없다. NURBS 곡선 투영법을 이용하여 트림되는 요소를 찾고, 트림된 요소는 쿼드트리 분할법과 NEFEM에서 사용된 적분방법을 동시에 고려하면 어떤 경우의 트림 요소라도 적분이 가능하다. 다양한 수치 예제를 통하여 제안한 해석 방법을 검증하고, 기존의 등기하해석법으로 해석하기 어려운 다수의 트림영역이 존재하는 NURBS 평면을 해석하여 본 방법의 유용성을 검토한다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.3
• /
• pp.151-157
• /
• 2021

This paper is devoted to a convergence study of the nonlocal integral operator in peridynamics. The implicit formulation can be an efficient approach to obtain the static/quasi-static solution of crack propagation problems. Implicit methods require constly large-matrix operations. Therefore, convergence is important for improving computational efficiency. When the radial influence function is utilized in the nonlocal integral equation, the fractional Laplacian integral equation is obtained. It has been mathematically proved that the condition number of the system matrix is affected by the order of the radial influence function and nonlocal horizon size. We formulate the static crack problem with peridynamics and utilize Newton-Raphson methods with a preconditioned conjugate gradient scheme to solve this nonlinear stationary system. The convergence behavior and the computational time for solving the implicit algebraic system have been studied with respect to the order of the radial influence function and nonlocal horizon size.

• Computational Structural Engineering
• /
• v.8 no.4
• /
• pp.137-148
• /
• 1995

A new path independent contour integral formulus for the distinct calculation of mode I stress intensity factors in two dimensional linear elastic fracture mechanics problems is presented. This method is based on p-convergence concepts and can be easily appended to existing finite element computer codes. In this study, the stress state at crack tip has been investigated and the path independence of J-integral values has been tested with respect to different contours expressed by normalized distance apart from the crack tip. Numerical results by p-convergence for the problems such as centrally cracked panels, single and double edged cracks in rectangular panels have been compared with those by the conventional h-convergence. The comparison demonstrates the accuracy and stability of the proposed method.

• Proceedings of the Korean Institute of Information and Commucation Sciences Conference
• /
• 2007.10a
• /
• pp.219-222
• /
• 2007

The paper deals with the rigorous analysis of three layered media structures. The dyadic Green's function for three layer medium is derived. The Green's functions belonging to the kernel of the integral equation are expressed as Sommerfeld integrals, in which surface wave effects are automatically included. We propose this integral representation as the most appropriate in the spatial domain analysis of slive structure. Also, we used extraction method for the convergence of this integral function. Finally, some numerical results are presented. These computed value show good agreement with proposed this method.

• Journal of Advanced Navigation Technology
• /
• v.7 no.2
• /
• pp.180-190
• /
• 2003

By using accurate closed-form Green's functions obtained from real-axis integration method, the full-wave analysis of CPW discontinuities are performed in space domain. In solving MPIE(Mixed Potential Integral Equation), Galerkin's scheme is employed with the linear basis functions on the triangular elements in air-dielectric boundary. In the singular integral arising when the observation point and source point coincides, the surface integral is transformed into the line integral and the integral is evaluated by regular integration. By using the Green's function from the real-axis integration method, the discontinuities are characterized accurately.