• Journal of Mechanical Science and Technology
• /
• v.17 no.1
• /
• pp.64-76
• /
• 2003

In this work, a new penalty formulation is proposed for the analysis of Mindlin-Reissner plates by using the element-free Galerkin method. A penalized weak form for the Mindlin-Reissner Plates is constructed through the exterior penalty method to enforce the essential boundary conditions of rotations as well as transverse displacements. In the numerical examples, some typical problems of Mindlin-Reissner plates are analyzed, and parametric studies on the order of integration and the size of influence domain are also carried out. The effect of the types of background cells on the accuracy of numerical solutions is observed and a proper type of background cell for obtaining optimal accuracy is suggested. Further, optimal order of integration and basis order of Moving Least Squares approximation are suggested to efficiently handle the irregularly distributed nodes through the triangular type of background cells. From the numerical tests, it is identified that unlike the finite element method, the proposed element-free Galerkin method with penalty technique gives highly accurate solution without shear locking in dealing with Mindlin-Reissner plates.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.10a
• /
• pp.615-622
• /
• 2002

In this paper, the modification of double projection method for the adaptive analysis of Element-free Galerkin(EFG) method were proposed. As results of the double projection method, the smoothed error profile that is adequate for adaptive analysis was obtained by re-projection of error that means the differences of EFG stress and projected stress. However, it was found that the efficiency of double projection method is degraded as increase of the numerical integration order. Since, the iterative refinement to single step error estimation made the same effect as increasing of integration order, the application of the iterative refinement base on double projection method could be produced the inadequately refined analysis model. To overcome this defect, a modified scheme of double projection were proposed. In the numerical example, the results did not show degradation of double projection effect in iterative refinement and the efficiency of proposed scheme were proved.

• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.625-630
• /
• 2000

In this paper, the governing equation and the boundary conditions of deploying rods are derived by using Hamilton's principle. The Galerkin method using the comparison function of the instantaneous natural modes is adopted by which the governing equation is discretized. Based on the discretized equations, the time integration analysis is performed and the longitudinal vibrations for the deploying and the retrieving velocity are analyzed.

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

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.4
• /
• pp.485-500
• /
• 2000

In this paper the adaptive crack propagation analysis based on the estimated local and global error in the element-free Galerkin (EFG) method is presented. It is possible to keep consistency and accuracy of analysis in each propagation step by adaptive analysis. The adaptivity analysis in crack propagation is achieved by adding and removing the node along the background integration cell that are refined or recovered as estimated error. These errors are obtained by calculating the difference between the values of the projected stresses and original EFG stresses. To evaluate the performance of proposed adaptive procedure, the convergence behavior is investigated lot several examples. The results of these examples show the efficiency of proposed scheme in crack propagation analysis.

• Interaction and multiscale mechanics
• /
• v.1 no.3
• /
• pp.339-355
• /
• 2008

A Galerkin meshfree method is presented for analyzing shear deformable cylindrical panels. Based upon the analogy between the cylindrical panel and the curved beam a pure bending mode for cylindrical panel is rationally constructed. The meshfree approximation employed herein is characterized by an enhanced moving least square or reproducing kernel basis function that can exactly represent the pure bending mode and thus meets the requirement of Kirchhoff mode reproducing condition. The variational form is discretized using the efficient stabilized conforming nodal integration with a smoothed nodal gradient based curvature. The resulting meshfree formulation satisfies the integration constraint for bending exactness. Moreover, it is shown here that the smoothed gradient preserves several desired properties which are valid for the standard gradient obtained by direct differentiation, such as partition of nullity and reproduction of a constant strain field. The efficacy of the proposed approach is demonstrated by two benchmark cylindrical panel examples.

• Structural Engineering and Mechanics
• /
• v.59 no.5
• /
• pp.901-920
• /
• 2016

MeshFree methods have become popular owing to the ease with which high stress gradients can be identified and node density distribution can be reformulated to accomplish faster convergence. This paper presents a strategy for nodal density refinement with strain energy as basis in Element-Free Galerkin MeshFree technique. Two popular flat plate problems are considered for the demonstration of the proposed strategies. Issue of integration errors introduced during nodal density refinement have been addressed by suggesting integration cell refinement. High stress effects around two symmetrical semi-circular notches under in-plane axial load have been addressed in the first problem. The second considers crack propagation under mode I and mode II fracture loading by the way of introducing high stress intensity through line crack. The computational efficacy of the adaptive refinement strategies proposed has been highlighted.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.474-479
• /
• 2004

According to our previous study, it is confirmed that the Petrov-Galerkin natural element method (PGNEM) completely resolves the numerical integration inaccuracy in the conventional Bubnov-Galerkin natural element method (BG-NEM). This paper is an extension of PG-NEM to two-dimensional nonlinear dynamic problem. For the analysis, a constant average acceleration method and a linearized total Lagrangian formulation is introduced with the PG-NEM. At every time step, the grid points are updated and the shape functions are reproduced from the relocated nodal distribution. This process enables the PG-NEM to provide more accurate and robust approximations. The representative numerical experiments performed by the test Fortran program, and the numerical results confirmed that the PG-NEM effectively and accurately approximates the nonlinear dynamic problem.

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