• Journal of Mechanical Science and Technology
• /
• v.18 no.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
• /
• v.20 no.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.

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

• Coupled systems mechanics
• /
• v.8 no.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.

• Structural Engineering and Mechanics
• /
• v.30 no.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.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1998.04a
• /
• 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
• /
• v.34 no.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.

• Kyungpook Mathematical Journal
• /
• v.61 no.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.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2012.05a
• /
• pp.521-522
• /
• 2012

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.24 no.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.