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

• Structural Engineering and Mechanics
• /
• v.85 no.4
• /
• pp.485-500
• /
• 2023

The non-linear static analysis of reinforced concrete (RC) structures using the three-dimensional (3D) finite element method is a time-consuming and challenging task. Moreover, this type of analysis encounters numerical problems such as the lack of convergence of results in the stages of growth and propagation of cracks in the structure. The time integration analysis along with the mass scaling (MS) technique is usually used to overcome these limitations. Despite the use of this method in the 3D finite element analysis of RC structures, a comprehensive study has not been conducted so far to assess the effects of the MS method on the accuracy of results. This study aims to evaluate the accuracy of the MS method in the non-linear quasi-static finite element analysis of RC structures. To this aim, different types of RC structures were simulated using the finite element approach based on the implicit time integration method and the mass scaling technique. The influences of effective parameters of the MS method (i.e., the allowable values of increase in the mass of the RC structure, the relationship between the duration of the applied load and fundamental vibration period of the RC structure, and the pattern of applied loads) on the accuracy of the simulated results were investigated. The accuracy of numerical simulation results has been evaluated through comparison with existing experimental data. The results of this study show that the achievement of accurate structural responses in the implicit time integration analyses using the MS method involves the appropriate selection of the effective parameters of the MS method.

• Journal of computational fluids engineering
• /
• v.12 no.3
• /
• pp.29-40
• /
• 2007

An implicit discontinuous Galerkin method for the two-dimensional Euler equations was developed on unstructured triangular meshes. The method can achieve high-order spatial accuracy by using hierachical basis functions based on Legendre polynomials. Numerical tests were conducted to estimate the convergence order of numerical solutions to the Ringleb flow and the supersonic vortex flow for which analytic solutions are available. Also, the flows around a 2-D circular cylinder and an NACA0012 airfoil were numerically simulated. The numerical results showed that the implicit discontinuous Galerkin methods couples with a high-order representation of curved solid boundaries can be an efficient method to obtain very accurate numerical solutions on unstructured meshes.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.30-40
• /
• 2007

The implicit discontinuous Galerkin method for the two-dimensional Euler equations was developed on unstructured triangular meshes, which can achieve higher-order accuracy by wing hierachical basis functions based on Legendre polynomials. Numerical tests were conducted to estimate the convergence order of numerical solutions to the Ringleb flow and the supersonic vortex flow for which analytic solutions are available. And, the flows around a circle and a NACA0012 airfoil was also numerically simulated. Numerical results show that the implicit discontinuous Galerkin methods with higher-order representation of curved solid boundaries can be an efficient higher-order method to obtain very accurate numerical solutions on unstructured meshes.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.25 no.4
• /
• pp.262-295
• /
• 2021

This paper introduces mono-implicit general linear methods, a special class of general linear methods, which are implicit in the output solution for the numerical integration of differential algebraic equations. We show how L-stable inherent Runge-Kutta members can be derived. The procedures for implementation have been discussed. The numerical test on the problem considered shows that the methods have improved accuracy when compared to RADAU IIA and the results from MATLAB ode15s, which have been taken as our reference solution.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.25 no.4
• /
• pp.224-261
• /
• 2021

The interest in implicit-explicit (IMEX) integration methods has emerged as an alternative for dealing in a computationally cost-effective way with stiff ordinary differential equations arising from practical modeling problems. In this paper, we introduce implicit-explicit second derivative linear multi-step methods (IMEX SDLMM) with error control. The proposed IMEX SDLMM is based on second derivative backward differentiation formulas (SDBDF) and recursive SDBDF. The IMEX second derivative schemes are constructed with order p ranging from p = 1 to 8. The methods are numerically validated on well-known stiff equations.

• Land and Housing Review
• /
• v.1 no.1
• /
• pp.59-65
• /
• 2010

The paper presents a comparison between a semi-implicit time integration linear finite element implementation and fully-implicit nonlinear Newton-Raphson finite element implementation of a triphasic small strain mixture formulation of an elastic partially saturated porous medium. The pore air phase pressure pa is assumed atmospheric, i.e., $p_a$ = 0, although the formulation and implementation are general to handle increase in pore air pressure as a result of loading, if needed. The solid skeleton phase is assumed linear isotropic elastic and partially saturated 'consolidation' in the presence of surface infiltration and traction is simulated. The verification of the implementation against an analytical solution for partially saturated pore water flow (no deformation) and comparison between the two implementations is presented and the important of the porosity-dependent nature of the partially saturated permeability is assessed on comparison with a commercial code for the partially saturated flow with deformation. As a result, the response of partially saturated permeability subjected to the porosity influences on the saturation of a soil, and the different behaviors of the partially saturated soil between staggered and monolithic coupled programs is worth of attention because the negative pore water pressure in the partially saturated soil depends on the difference.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.6 s.177
• /
• pp.1499-1512
• /
• 2000

The dynamic explicit finite element method and the static implicit finite element method are applied effectively to analyze total auto-body panel stamping processes, which include the forming stage , the trimming stage and the spring-back stage.\The explicit time integration method has better merits in the forming stage including highly complicated three-dimensional contact conditions. On the contrary, the implicit time integration method is better for analyzing spring-back since the complicated contact conditions are removed and the computing time to get the final static state is short. In this work, brief descriptions of the formulation and the factor study are presented. Further, the simulated results for the total auto-body panel stamping processes are shown and discussed. The formability and the weld line movement in stamping with Tailor Welded Blanks were investigated through QTR-OTR-FRT.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.12
• /
• pp.1966-1973
• /
• 2001

The Jacobian matrix of the equations of motion of a system using velocity transformation technique is derived via variation methods to apply the implicit integration algorithm, DASSL. The concept of generalized coordinate partitioning is used to parameterize the constraint set with independent generalized coordinates. DASSL is applied to determine independent generalized coordinates and velocities. Dependent generalized coordinates, velocities, accelerations and Lagrange multipliers are explicitly retained in the formulation to satisfy all of the governing kinematic and dynamic equations. The derived Jacobian matrix of a system is proved to be valid and accurate both analytically and through solution of numerical examples.

• Journal of the Korea Computer Graphics Society
• /
• v.10 no.2
• /
• pp.17-26
• /
• 2004

본 연구는 적은 감쇠(damping) 효과를 가지면서도 안정적으로 의복의 움직임을 생성할 수 있는 시뮬레이션 방법을 제안한다. 최근 의복 시뮬레이션 분야에서는 뛰어난 수치적 안정성을 지닌 암시적 적분법(implicit integration)이 널리 이용되어 왔다. 그러나, 암시적 적분법은 선형화 과정에서 과도한 감쇠효과가 발생하고, 이 감쇠효과는 내부 변형 없는 회전이동 시에도 과도하게 나타나므로, 옷감의 회전 이동이나 주름의 자연스런 움직임을 크게 방해한다. 본 연구는 불필요한 감쇠효과를 줄이고 안정적으로 의복의 움직임을 표현할 수 있는 새로운 암시적 적분법을 제안한다. 제안하는 방법용 선형화 과정에서 발생하는 내부 감쇠력을 회전이용이 아닌 순수한 내부 변형에 대해서만 발생하도록 모델링하며, 이렇게 계산된 내부 감쇠력을 역학식에 안정적으로 반영한다. 실험 결과는 제안하는 방법이 기존의 방법들에 비해 보다 다양한 재질의 옷감을 과도한 감쇠 효과 없어도 안정적으로 실시간에 생성할 수 있음을 보여주었다.