• Journal of the Korea institute for structural maintenance and inspection
• /
• v.20 no.5
• /
• pp.93-101
• /
• 2016

In this study, the governing equations of motion for multi-cracked nonuniform nanobeam based on nonlocal elasticity theory and embedded in an elastic medium were derived. DTM(differential transformation method) was applied to vibration analysis of multi-cracked nonuniform nanobeam based on nonlocal elasticity theory and embedded in an elastic medium. The non-dimensional natural frequencies of this nanobeam were obtained for eoe, crack stiffness and elastic medium stiffness with various boundary conditions. The results obtained by this method was compared with previous works and showed the close agreement between two methods. The important conclusions obtained by this study are as follows : 1. As the length of nanobeam is shorter, the effect of scale coefficient is greater. 2. The locations of crack change non-dimensional natural frequency, In the case of fixed-fixed ends, the non-dimensional natural frequency is the biggest in the first crack location of 0.6L of nanobeam length, and the smallest in both ends. In the case of fixed-free ends, the closer the location of first crack go tho the free end, the bigger the non-dimensional natural frequency. 3. As the stiffness of crack is greater, the non-dimensional natural frequency is smaller, And the effect of crack stiffness is similar on both fixed-free ends and fixed-fixed ends. 4. The bigger the stiffness of elastic medium, the greater the non - dimensional natural frequency.

• Journal of the Korean Geotechnical Society
• /
• v.22 no.9
• /
• pp.45-59
• /
• 2006

The successful performance of any numerical geotechnical simulation depends on the accuracy and efficiency of the numerical implementation of constitutive model used to simulate the stress-strain (constitutive) response of the soil. The corner stone of the numerical implementation of constitutive models is the numerical integration of the incremental form of soil-plasticity constitutive equations over a discrete sequence of time steps. In this paper a well known two-surface soil plasticity model is implemented using a generalized implicit return mapping algorithm to arbitrary convex yield surfaces referred to as the Closest-Point-Projection method (CPPM). The two-surface model describes the nonlinear behavior of coarse-grained materials by incorporating a bounding surface concept together with isotropic and kinematic hardening as well as fabric formulation to account for the effect of fabric formation on the unloading response. In the course of investigating the performance of the CPPM integration method, it is proven that the algorithm is an accurate, robust, and efficient integration technique useful in finite element contexts. It is also shown that the algorithm produces a consistent tangent operator $\frac{d\sigma}{d\varepsilon}$ during the iterative process with quadratic convergence rate of the global iteration process.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.1
• /
• pp.29-38
• /
• 2013

This paper deals with geometrical non-linear analyses of the tapered cantilever column subjected to the sub-tangential follower force at the free end. Cross-sections of the column whose flexural rigidities are functionally varied with the axial coordinate. The differential equations governing the elastica of such column are derived on the basis of the large deformation theory. These differential equations have three unknown parameters of the vertical and horizontal deflections and rotation at the free end. These differential equations are numerically solved by the iteration technique for obtaining three unknowns and elastica of the deformed column. For validating theories developed herein, laboratory scaled experiments are conducted.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.37 no.3
• /
• pp.187-195
• /
• 2024

This paper presents two-dimensional boundary value problems of the stress-based gradient elasticity within the smoothed finite element method (S-FEM) framework. Gradient elasticity is introduced to address the limitations of classical elasticity, particularly its struggle to capture size-dependent mechanical behavior at the micro/nano scale. The Ru-Aifantis theorem is employed to overcome the challenges of high-order differential equations in gradient elasticity. This theorem effectively splits the original equation into two solvable second-order differential equations, enabling its incorporation into the S-FEM framework. The present method utilizes a staggered scheme to solve the boundary value problems. This approach efficiently separates the calculation of the local displacement field (obtained over each smoothing domain) from the non-local stress field (computed element-wise). A series of numerical tests are conducted to investigate the influence of the internal length scale, a key parameter in gradient elasticity. The results demonstrate the effectiveness of the proposed approach in smoothing stress concentrations typically observed at crack tips and dislocation lines.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2010.05a
• /
• pp.213-217
• /
• 2010

댐붕괴파 (dam-break flow)나 지진해일에 의해 발생하는 undular bore와 충격파 (shock) 현상을 동수압 및 분산효과를 고려하여 수치모의를 수행하였다. 완전비선형 Boussinesq-type equations 모형을 이용하여, 동수압 및 분산 효과를 고려하였다. 방정식은 4차 정확도의 유한체적법을 이용하여 해석하였고, 시간적으로도 4차정확도의 기법을 이용하여 고차미분항에 대한 수치분산을 억제하였다. 다양한 경우의 1차원과 2차원 공간에서의 수치모의를 수행하고 검증을 수행하였다. 그 결과, 완전비선형 Boussinesq-type equations 모형은 천수방정식 (shallow water equations) 기반의 모형에서 재현이 불가능한 undular bore 등을 재현 하는 등, 전반적으로 천수방정식 기반의 모형 보다 물리적으로도 타당하고 정량적으로도 실험결과와 잘 일치하는 경향을 보였다. 즉, 댐붕괴파나 지진해일 등에 의한 범람 모의에 있어 동수압과 분산 효과의 중요성이 공학적으로도 매우 중요한 고려사항 임이 나타났다.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.142-149
• /
• 1996

이차원 노즐을 통하여 저밀도 환경으로 팽창하는 희박류의 분석에 있어서 불연속좌표법과 결합된 유한차분법(finite-difference method coupled with the discrete-ordinate method, FDDO)과 직접모사법(direct-simulation Monte-Carlo method, DSMC)이 비교되었다. FDDO를 이용한 분석에서는 충돌적분모델을 도입하여 간단해진 볼츠만식(Boltzmann equation)이 불연속좌표법을 이용하여 물리적 공간에서는 연속이나 분자속도 공간에서는 불연속좌표로 표시되는 편미분방정식군으로 변환되어 유한차분법에의하여 수치해석 되었다. 직접모사법에서는 분자모델로 가변강구모델(variable hard sphere model, VHS)이, 충돌샘플링모델로는 비시계수법(no time counter method, NTC)이 채택되었다. 전혀 다른 두 가지 방법에 의한 노즐 내부에서의 유체흐름 해석결과는 매우 잘 일치하였으며, 노즐 외부의 plume 영역에서는 FDDO에 의한 해석결과가 직접모사법에 의한 해석결과에 비하여 약간 느린 팽창을 보였다.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.4
• /
• pp.48-59
• /
• 1996

The numerical step in the unsteady viscous flow analysis can be divided in the space analysis step satisfying continuity equation and the time marching step. In this study the spectral method is applied to solve the pressure Poisson equation in the space analysis step. If the highest order differential term of the pressure Poisson equation is transformed by Fourier series, pressure arid its first derivatives can be expressed by the integral form of Fourier series. So Gibb's phenomena can be eliminated and the spectral method can be applied to non-periodic problems. The numerical analysis of unsteady viscous flow around 2-dimensional circular cylinder and wing is carried out and compared for verification.

• Proceedings of the KIEE Conference
• /
• 2003.07d
• /
• pp.2594-2596
• /
• 2003

유한요소법(FEM)은 많은 공학문제를 해결하는 가장 중요한 방법 중 하나로 인식되고 있다. 본 논문에서는 자동제어 및 신호처리 문제해결에 효율적이며 강력한 수치해석 패키지인 CEMTool환경에서 유한요소법을 이용하여 일반적인 편미분방정식 Solver 구현에 관한사항을 논의하고자 한다. 기본적으로 영역정보 및 노드수 등의 정보를 입력받아 각 노드의 정보를 출력하는 Mesh함수를 구현하며, 생성된 요소방정식들을 조립하는 Assemble함수를 작성한 뒤, Boundary함수를 통해 경계조건을 적용시킨 후 선형행렬 방정식을 풀어 전체노드의 값을 찾아내는 Solve함수를 구현하는 과정을 알아본다. 구현된 FEM Solver의 전체적인 구조를 통해 구현시 고려해야 할 사항을 논의하며 기본적인 편미분방정식의 예제를 통해 FEM PDE Solver의 동작과정을 검증할 것이다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.4
• /
• pp.215-222
• /
• 2019

The heat conduction problem with discontinuous material coefficients generally consists of the conservative equation, boundary condition, and interface condition, which should be additionally satisfied in the solution procedure. This feature often makes the development of new numerical schemes difficult as it induces a layered singularity in the solution fields; thus, a special approximation is required to capture the singular behavior. In addition to the approximation, the construction of a total system of equations is challenging. In this study, a wedge function is devised for enriching the approximation, and the interface condition itself is embedded in the moving least squares(MLS) derivative approximation to consistently satisfy the interface condition. The heat conduction problem is then discretized in a strong form using the developed derivative approximation, which is named as the interface immersed MLS difference method. This method is able to efficiently provide a numerical solution for such interface problems avoiding both numerical quadrature as well as extra difference equations related to the interface condition enforcement. Numerical experiments proved that the developed numerical method was highly accurate and computationally efficient at solving the heat conduction problem with interfacial jump as well as the problem with a geometrically induced interfacial singularity.

• Journal of the Society of Disaster Information
• /
• v.10 no.1
• /
• pp.123-140
• /
• 2014

This paper proposes a simple and effective method for determining the dynamic characteristics of revolution shells. This is a weighted residual method in which the collocation points are taken at the roots of orthogonal polynomial. In this paper the collocation method is employed to replace a partical differential eqations by a system of ordinary differential equations in time, and the resulting equations are solved by two different numerical methods of time integration : an implicit method and an explicit method. The proposed approach is formulated in some detail. The versatility and accuracy are illustrated through several numerical examples. The method appears to be relatively easy to set up and gives satisfactory results.