• Journal of the Earthquake Engineering Society of Korea
• /
• v.7 no.5
• /
• pp.11-17
• /
• 2003

It is well known that the double integration of measured acceleration records is one of the most difficult signal processing, particularly in the measurements on civil engineering structures, The measured accelerations of civil engineering structures are usually non-stationary and contain non-gaussian low-frequency noises, which can be significant causes of numerical instabilities in double Integration, For the de-noising of this kind of signals, wavelet transform can be very effective because of its inherent processing features for non-stationary signals, In this paper, the de-noising algorithm for the double integration is proposed using the modified wavelet transform, which is extended version of ordinary wavelet transform to process non-gaussian and low-frequency noises, using the median filter concept, The example studies show that the integration can be improved by the proposed method.

• Composites Research
• /
• v.23 no.4
• /
• pp.7-13
• /
• 2010

A volume integral equation method (VIEM) is introduced for the solution of elastostatic problems in an unbounded isotropic elastic solids containing interacting multiple anisotropic inclusions subject to remote uniaxial tension. The method is applied to two-dimensional problems involving long parallel cylindrical inclusions. A detailed analysis of stress field at the interface between the matrix and the central inclusion is carried out for square and hexagonal packing of the inclusions. Effects of the number of anisotropic inclusions and various fiber volume fractions on the stress field at the interface between the matrix and the central inclusion are also investigated in detail. The accuracy of the method is validated by solving the single inclusion problem for which solutions are available in the literature.

• 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.

• Geophysics and Geophysical Exploration
• /
• v.26 no.1
• /
• pp.1-7
• /
• 2023

This study derives the expressions of vector gravity and gravity gradient tensor due to an elliptical cylinder. The vector gravity for an arbitrary three-dimensional (3D) body is obtained by differentiating the gravitational potential, including the triple integral, according to the shape of the body in each axis direction. The vector gravity of the 3D body with axial symmetry is integrated along the axial direction and reduced to a double integral. The complex Green's theorem using complex conjugates subsequently converts the double integral into a one-dimensional (1D) closed-line integral. Finally, the vector gravity due to the elliptical cylinder is derived using 1D numerical integration by parameterizing a boundary of the elliptical cross-section as a closed line. Similarly, the gravity gradient tensor due to the elliptical cylinder is second-order differentiated from the gravitational potential, including the triple integral, and integrated along the vertical axis direction reducing it to a double integral. Consequently, all the components of the gravity gradient tensor due to an elliptical cylinder are derived using complex Green's theorem as used in the case of vector gravity.

• Journal of Ocean Engineering and Technology
• /
• v.12 no.1
• /
• pp.85-98
• /
• 1998

본 논문의 수치해법은 경계치문제를 풀기 위하여 코시이론(Cauchy's theorem)을 사용하였다. 경계치문제는 완전한 물체표면조건과 자유표면조건을 만족시키는 초기치문제로 귀결된다. 현 수치해법에서 무한영역은 수치계산 영역인 비선형 영역과 선형 자유표면조건을 만족하는 선형영역으로 나누어진다. 선형영역의 해는 과도 그린(Green)함수를 사용하여 정합조건을 부과함으로써, 수치계산은 비선형 영역에서만 수행된다. 본 논문에서 저자는 수치계산 영역에서 코시이론을 사용하여 적분방정식을 도출하였고, 무한영역의 해는 정합면에서 과도 그린함수를 사용하여 표현하였다. 본 수치계산에서 자유표면에 요소 재분배법을 적용함으로써 쇄파현상에 대해서도 안정적인 수치해석을 할 수 있었다. 본 논문에서 개발된 수치방법을 적용한 문제는 다음과 같다. 첫째는 자유표면에서 실린더가 강제동요하는 경우에 자유표면형상과 힘을 계산하여 이전의 실험치 및 계산치와 비교하였다. 두번째로는 실린더가 자유수면하에서 일정한 속도로 항주하는 경우에는 조파저항과 양력을 계산하여 고차 스펙트럴법과 비교하였다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.5
• /
• pp.525-533
• /
• 2010

Direct time integration method such as Newmark method is numerically performed under the assumption that continuous load function such as constant amplitude load can be treated as a discontinuous load fuction. It is due that the load can be treated as a constant value at the given time period regardless of variation of load at the time increment interval. It means the numerical results should be accompanied by the error due to approximation of load fuction. In contrast, the load function is calculated by convolution integral for the given time interval at finite element equation based on Gurtin's variation equation. Therefore. precise numerical results can be obtained by Gurtin's method because of convolution integral for the continuous load fuction curve even at the variation of load function in the given time interval. In this study, we prove that Gurtin's method can be more suitable than Newmark method in the problem of constant amplitude loading, using the numerical results for the free end of the one-dimensional rod. This study also shows that Gurtin's method is more effective in constant amplitude loading than in constant loading. The accuracy and the validity are verified by comparison between the results of in-house FORTRAN code and ADINA, a commercial software supporting Newmark method.

• Journal of the Korean GEO-environmental Society
• /
• v.15 no.3
• /
• pp.41-48
• /
• 2014

Experimental study of sedimentation and self-weight consolidation has been primary research area in dredged soil. However, good quality of the dredged soil and minimum water pollution caused by the pumping of reclaimed soil require intensive study of the flow characteristics of dredged material due to dumping. In this study, continuity and the equilibrium equations for mass flow assuming single phase was derived to simulate mass flow in dredged containment area. To optimize computation and modeling time for three dimensional geometry and boundary conditions, depth integration is applied to governing equations to consider three dimensional topography of the site. Petrov-Galerkin formulation is applied in spatial discretization of governing equations. Generalized trapezoidal rule is used for time integration, and Newton iteration process approximated the solution. DG and CDG technique were used for weighting matrix in discontinuous test function in dredged flow analysis, and numerical stability was evaluated by performed a square slump simulation. A comparative analysis for numerical methods showed that DG method applied to SU / PG formulation gives minimal pseudo oscillation and reliable numerical results.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.16 no.4
• /
• pp.37-52
• /
• 2012

In this study, a new 4-node general shell element with 6 DOFs per node is presented. Drilling rotational degrees of freedom are introduced by the variational principle with an independent rotation field. In formulation of the element, substitute transverse shear strain fields are used to avoid shear locking, while four nonconforming modes are applied in the in-plane displacement fields as a remedy for membrane locking. In addition, a direct modification method for nonconforming modes is employed in the numerical implementation of nonconforming modes to represent constant strain states. A 9-points integration rule is adopted for volume integration in the computation of the element stiffness matrix. With the combined use of these techniques, the developed shell element has no spurious zero energy modes, and can represent a constant strain state. Several numerical tests are carried out to evaluate the performance of the new element developed. The test results show that the behavior of the elements is satisfactory.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.2
• /
• pp.113-121
• /
• 2005

In order to resolve a common numerical integration inaccuracy of meshfree methods, we introduce an improved natural clement method called Petrov-Galerkin natural element method(PG-NEM). While Laplace basis function is being taken for the trial shape function, the test shape function in the present method is differently defined such that its support becomes a union of Delaunay triangles. This approach eliminates the inconsistency of tile support of integrand function with the regular integration domain, and which preserves both simplicity and accuracy in the numerical integration. In this paper, the validity of the PG-NEM is verified through the representative benchmark problems in 2-d linear elasticity. For the comparison, we also analyze the problems using the conventional Bubnov-Galerkin natural element method(BG-NEM) and constant strain finite clement method(CS-FEM). From the patch test and assessment on convergence rate, we can confirm the superiority of the proposed meshfree method.

• Journal of Korea Game Society
• /
• v.18 no.5
• /
• pp.123-132
• /
• 2018

Various methods have been proposed to efficiently animate the motion of soft objects in realtime. In order to maintain the topology between the elements of the objects, it is required to employ constraint forces, which limit the size of the time steps for the numerical integration and reduce the efficiency. To tackle this, an implicit method with larger steps was proposed. However, the method is, in essence, a linear system with a large matrix, of which solution requires heavy computations. Several approximate methods have been proposed, but the approximation is obtained with an increased damping and the loss of accuracy. In this paper, new integration method based on harmonic oscillation with better stability was proposed, and it was further stabilized with the hybridization with approximate implicit method. GPU parallelism can be easily implemented for the method, and large-scale soft objects can be simulated in realtime.