• Journal of Mechanical Science and Technology
• /
• 제18권3호
• /
• pp.460-470
• /
• 2004

A computational analysis of engineering problems with moving domain or/and boundary according to either Lagrangian or Eulerian approach may encounter inherent numerical difficulties, the extreme mesh distortion in the former and the material boundary indistinctness in the latter. In order to overcome such defects in classical numerical approaches, the ALE(arbitrary Lagrangian Eulerian) method is widely being adopted in which the finite element mesh moves with arbitrary velocity. This paper is concerned with the ALE finite element formulation, aiming at the dynamic response analysis of baffled fuel-storage container in yawing motion, for which the coupled time integration scheme, the remeshing and smoothing algorithm and the mesh velocity determination are addressed. Numerical simulation illustrating theoretical works is also presented.

• Steel and Composite Structures
• /
• 제20권2호
• /
• pp.295-316
• /
• 2016

This paper presents the results of experimental investigation, numerical calculation and theoretical analysis on axial compression ratio limit values for steel reinforced concrete (SRC) special shaped columns. 17 specimens were firstly intensively carried out to investigate the hysteretic behavior of SRC special shaped columns subjected to a constant axial load and cyclic reversed loads. Two theories were used to calculate the limits of axial compression ratio for all the specimens, including the balanced failure theory and superposition theory. It was found that the results of balanced failure theory by numerical integration method cannot conform the reality of test results, while the calculation results by employing the superposition theory can agree well with the test results. On the basis of superposition theory, the design limit values of axial compression ratio under different seismic grades were proposed for SRC special shaped columns.

• 한국지반공학회논문집
• /
• 제22권9호
• /
• pp.45-59
• /
• 2006

탄소성 구성방정식은 주로 미분방적식(rate equation)으로 이루어져 있기 때문에 유한요소법 등을 이용한 지반구조물 해석시 미분방정식들에 대한 수치적분을 수행할 수 있는 방법이 필요하다. 구조물의 거동을 해석할시 미분방정식들을 위한 적분방법은 해석결과의 정확성과 유한요소법 모델링의 안전성에 큰 영향을 미치고 있다. 본 논문에서는 최근에 개발되어 사용되고 있는 흙에 관한 구성모델인 "Two-surface soil plasticity model (Manzari and Dafalias 1997)"을 Implicit return-mapping 수치적분방법을 이용하여 실행하는 과정을 제시한다. 본 연구에서 사용된 수치적분방법은 Closest-Point-Projection Method(CPPM) 방법으로 탄성 예측자-소성 교정자(elastic predictor-plastic corrector) 개념을 Implicit Backward Euler방법으로 체계화 시킨 알고리듬이다. 본 연구에서 수행한 "Two-surface soil plasticity model"은 조립토의 비선형거동을 해석하며, Bounding surface 개념 및 비선형 등방경화와 이동경화법칙을 사용하는 모델이다. 본 연구는 CPPM 방법이 정확하고 안정되며 유용한 수치적분을 수행할 수 있는 알고리듬이라는 것을 제시한다. 또한, CPPM 알고리듬은 구성방정식의 해를 반복적으로 해석하는 동안 "Consistent tangent operator $d{\sigma}/d{\varepsilon}$"를 제공하므로, 비선형 유한요소 해석이 2차(quadratic convergence rate)의 수렴 조건을 만족하는데 기여한다는 것을 보여준다.

• Coupled systems mechanics
• /
• 제8권1호
• /
• pp.71-93
• /
• 2019

Finite element simulations of solid mechanics problems often involve the use of Non-Confirming Meshes (NCM) to increase accuracy in capturing nonlinear behavior, including damage and plasticity, in part of a solid domain without an undue increase in computational costs. In the presence of material nonlinearity and plasticity, higher-order variables are often needed to capture nonlinear behavior and material history on non-conforming interfaces. The most popular formulations for coupling non-conforming meshes are dual methods that involve the interpolation of a traction field on the interface. These methods are subject to the Ladyzhenskaya-Babuska-Brezzi (LBB) stability condition, and are therefore limited in their implementation with the higher-order elements needed to capture nonlinear material behavior. Alternatively, the enriched discontinuous Galerkin approach (EDGA) (Haikal and Hjelmstad 2010) is a primal method that provides higher order kinematic fields on the interface, and in which interface tractions are computed from local finite element estimates, therefore facilitating its implementation with nonlinear material models. The inclusion of higher-order interface variables, however, presents the issue of preserving material history at integration points when a increase in integration order is needed. In this study, the enriched discontinuous Galerkin approach (EDGA) is extended to the case of small-deformation plasticity. An interface-driven Gauss-Kronrod integration rule is proposed to enable adaptive enrichment on the interface while preserving history-dependent material data at existing integration points. The method is implemented using classical J2 plasticity theory as well as the pressure-dependent Drucker-Prager material model. We show that an efficient treatment of interface variables can improve algorithmic performance and provide a consistent approach for coupling non-conforming meshes in inelasticity.

• Kyungpook Mathematical Journal
• /
• 제61권2호
• /
• pp.383-393
• /
• 2021

In this paper, numerical techniques are presented for solving initial value problems of fractional differential equations with variable coefficients. The method is derived by applying a Taylor vector approximation. Moreover, the operational matrix of fractional integration of a Taylor vector is provided in order to transform the continuous equations into a system of algebraic equations. Furthermore, numerical examples demonstrate that this method is applicable and accurate.

• Structural Engineering and Mechanics
• /
• 제30권3호
• /
• pp.279-296
• /
• 2008

This paper presents an exact integration for the hypersingular boundary integral equation of two-dimensional elastostatics. The boundary is discretized by straight segments and the physical variables are approximated by discontinuous quadratic elements. The integral for the hypersingular boundary integral equation analysis is given in a closed form. It is proven that using the exact integration for discontinuous boundary element, the singular integral in the Cauchy Principal Value and the hypersingular integral in the Hadamard Finite Part can be obtained straightforward without special treatment. Two numerical examples are implemented to verify the correctness of the derived exact integration.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 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.

• Journal of Advanced Marine Engineering and Technology
• /
• 제34권7호
• /
• pp.943-950
• /
• 2010

An efficient meshfree method based on a stabilized conforming nodal integration method is developed for elastoplastic contact analysis of metal forming processes. In this approach, strain smoothing stabilization is introduced to eliminate spatial instability in Galerkin meshfree methods when the weak form is integrated by a nodal integration. The gradient matrix associated with strain smoothing satisfies the integration constraint for linear exactness in the Galerkin approximation. Strain smoothing formulation and numerical procedures for path-dependent problems are introduced. Applications of metal forming analysis are presented, from which the computational efficiency has been improved significantly without loss of accuracy.

• 한국정밀공학회:학술대회논문집
• /
• 한국정밀공학회 2012년도 춘계학술대회 논문집
• /
• pp.521-522
• /
• 2012

• 한국소음진동공학회논문집
• /
• 제24권1호
• /
• pp.45-53
• /
• 2014

For the time integration solution of the impulsive dynamic contact problem of high-speed rotating disks formulated by the finite element technique, the velocity and acceleration contact constraints as well as the displacement contact constraint are imposed for the numerical stability without spurious oscillations. The solution of the present technique is checked by the numerical simulation using the concentric high-speed rotating disks with the clearance and impulsive loading. It is shown that the almost steady state solution agrees with the corresponding analytical solution of the elasticity and that the differentiated constraints are crucial for the numerical stability of such high-speed contact problems of the disks under impulsive loading.