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

• Proceedings of the Korean Institute of Electrical and Electronic Material Engineers Conference
• /
• 1995.11a
• /
• pp.48-51
• /
• 1995

We have developed the three dimensional mesh generator for three dimensional process simulation using the FEM(Finite Element Method). Tetrahedron element construct the presented three dimensional mesh, which is suitable for the simulation of three dimensional behavior of the LOCOS. The simulation of thermal oxidation is one of the problem in scale downed semiconductor processes. As three dimensional simulators use the huge size of the memory, we use the efficient method that generates the new nodes inside the growing oxide and removes the nodes nearby the SiO2/Si interface in silicon. The resented three dimensional mesh generator was designed to be used in various process simulations, for instance thermal oxidation, silicidation, nitridation, ion implantation, diffusion, and so on.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.6
• /
• pp.1004-1014
• /
• 2001

The fully automatic algorithm from initial finite element mesh generation to remeshing in two dimensional geometry is introduced using bubble packing method (BPM) for finite element analysis. BPM determines the node placement by force-balancing configuration of bubbles and the triangular meshes are made by Delaunay triangulation with advancing front concept. In BPM, we suggest two node-search algorithms and the adaptive/recursive bubble controls to search the optimal nodal position. To use the automatically generated mesh information in FEA, the new enhanced bandwidth minimization scheme with high efficiency in CPU time is developed. In the remeshing stage, the mesh refinement is incorporated by the control of bubble size using two parameters. And Superconvergent Patch Recovery (SPR) technique is used for error estimation. To verify the capability of this algorithm, we consider two elasticity problems, one is the bending problem of short cantilever beam and the tension problem of infinite plate with hole. The numerical results indicate that the algorithm by BPM is able to refine the mesh based on a posteriori error and control the mesh size easily by two parameters.

• Journal of Mechanical Science and Technology
• /
• v.20 no.1
• /
• pp.110-124
• /
• 2006

The multi scale wavelet-Galerkin method implemented in an adaptive manner has an advantage of obtaining accurate solutions with a substantially reduced number of interpolation points. The method is becoming popular, but its numerical efficiency still needs improvement. The objectives of this investigation are to present a new numerical scheme to improve the performance of the multi scale adaptive wavelet-Galerkin method and to give detailed implementation procedure. Specifically, the subdomain technique suitable for multiscale methods is developed and implemented. When the standard wavelet-Galerkin method is implemented without domain subdivision, the interaction between very long scale wavelets and very short scale wavelets leads to a poorly-sparse system matrix, which considerably worsens numerical efficiency for large-sized problems. The performance of the developed strategy is checked in terms of numerical costs such as the CPU time and memory size. Since the detailed implementation procedure including preprocessing and stiffness matrix construction is given, researchers having experiences in standard finite element implementation may be able to extend the multi scale method further or utilize some features of the multiscale method in their own applications.

• Computational Structural Engineering
• /
• v.1 no.2
• /
• pp.121-130
• /
• 1988

Adaptive reinements yield a large sparse system of equations. In order to solve such a system, the core storage requirement is an important consideration. Accordingly, an iterative method which minimizes the core storage and provides a high rate of convergence is called for. In this paper the conjugate gradient algorithms with various preconditionings including the incomplete Cholesky decomposition are examined.

• Journal of Korean Society for Atmospheric Environment
• /
• v.9 no.E
• /
• pp.347-357
• /
• 1993

The Lagrangian grid-free numerical method, the discrete vortex method, was applied to solve the Navier-Stokes euqations. This method avoids the introduction of numerical viscosity swamping the real physical viscosity at high Reynolds number, unlike Eulerian method, e.g. finite difference and element methods. The boundary integral equation method for the potential flow solution was included to make the discrete vortex method more feasible for complex geometries. The fast adaptive multipole expansion method was incorporated to reduce the computational time from $O(N^2)$ to O(N) for the computations of vortex-vortex interactions. The test problems were air flow around one circular cylinder and two circular cylinders in tandem with various gaps. The numerical results were in excellent gareement with the experimental and other computational results. The applicabilty of the method was discussed with the indoor and the outdoor air pollution problems, especially the contaminant transport in the recirculation regions.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.21 no.2
• /
• pp.283-294
• /
• 2023

APro, developed in KAERI for the process-based total system performance assessment (TSPA) of deep geological disposal systems, performs finite element method (FEM)-based multiphysics analysis. In the FEM-based analysis, the mesh element quality influences the numerical solution accuracy, memory requirement, and computation time. Therefore, an appropriate mesh structure should be constructed before the mesh stability analysis to achieve an accurate and efficient process-based TSPA. A generic reference case of DECOVALEX-2023 Task F, which has been proposed for simulating stationary groundwater flow and time-dependent conservative transport of two tracers, was used in this study for mesh stability analysis. The relative differences in tracer concentration varying mesh structures were determined by comparing with the results for the finest mesh structure. For calculation efficiency, the memory requirements and computation time were compared. Based on the mesh stability analysis, an approach based on adaptive mesh refinement was developed to resolve the error in the early stage of the simulation time-period. It was observed that the relative difference in the tracer concentration significantly decreased with high calculation efficiency.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.19 no.4 s.74
• /
• pp.429-440
• /
• 2006

This paper comprises two specific objectives. The first is to examine the applicability of ordinary kriging interpolation(OK) to the p-adaptivity of the finite element method that is based on variogram modeling. The second objective Is to present the adaptive procedure by the hierarchical p-refinement in conjunction with a posteriori error estimator using the modified S.P.R. (superconvergent patch recovery) method. The ordinary kriging method that is one of weighted interpolation techniques is applied to obtain the estimated exact solution from the stress data at the Gauss points. The weight factor is determined by experimental and theoretical variograms for interpolation of stress data apart from the conventional interpolation methods that use an equal weight factor. In the p-refinement, the analytical domain has to be refined automatically to obtain an acceptable level of accuracy by increasing the p-level non-uniformly or selectively. To verify the performance of the modified S.P.R. method, the new error estimator based on limit value has been proposed. The validity of the proposed approach has been tested with the help of some benchmark problems of linear elastic fracture mechanics such as a centrally cracked panel, a single edged crack, and a double edged crack.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2008.04a
• /
• pp.114-119
• /
• 2008

A refined three-dimensional finite element interaction model between the high-speed train and railway bride deck has been developed in the present study. Analytical predictions of vertical deflections for a railway bridge are compared with in-situ test results and a good agreement is achieved. Then, input variables employed in the analytical comparisons are selected as random variables for the limit state functions. followed by risk assessment. For this purpose, a linear adaptive weighted response surface method has been developed and applied. A typical railway bridge has been selected and the limit state functions are employed from UIC and Korean specifications in the comparative studies. The results reveal that Korean specifications give significantly risky reliability indices in comparison with UIC specifications. It is thus encouraged from the above that the present linear adaptive weighted response surface method can be an alternative for the fast estimation of nonlinear structural systems.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.3
• /
• pp.503-512
• /
• 1990

In this study, a penalty method effective for the case that has no snap-back phenomenon, is proposed and an adaptive method which choose the penalty method or Riks' type method, is suggested for the case of snap-back problem. And for the case that loads are applied to one or more points of a structure, the Riks' method is applied in general, but under certain condition choice of new incremental load parameter is suggested to accelerate the convergence rate. Finally, for the case that displacements of a structure are controlled at one or more points Riks' type method is proposed. In this case, the proposed method is applied in general but under certain condition it is recommended to choose other incremental displacement parameter to eliminate probable divergence. Five examples are analysed and compared with the result of published literature.