• Structural Engineering and Mechanics
• /
• v.7 no.6
• /
• pp.561-573
• /
• 1999

A $C^*$-convergence algorithm for finite element analysis has been proposed by Bigdeli and Kelly (1997) and elements for the first three levels applied to planar elasticity have been defined. The fourth level element for the new family is described in this paper and the rate of convergence for the $C^*$-convergence algorithm is investigated numerically. The new family adds derivatives of displacements as nodal variables and the number of nodes and elements can therefore be kept constant during refinement. A problem exists on interfaces where the derivatives are required to be discontinuous. This problem is addressed for curved boundaries and a procedure is suggested to resolve the excessive interelement continuity which occurs.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1996.10a
• /
• pp.210-217
• /
• 1996

Under strong shock loads including earthquake or blast, structures may start to crack in stress concentrated members. The continuous behavior of the structure changes to the discontinuous. In this study, numerical method analyzing continuous and discontinuous behavior of a structure is developed using a modified distinct element method. Equations of motion of each distinct element are integrated using the central difference method, one of the finite difference methods. Interactions between he elements are considered by an element and pore spring. The forces acting in the center of an element include contact stress transferred by element spring; tensile stress by pore spring; and external traction such as earthquake or blast load. To verify the proposed method, the behavior of the cantilever beam subject to the quasi-static concentrated force at the end is investigated. The failure behavior of the simply supported beam subject to the strong shock at the center is studied. The proposed method can predict the failure behavior of the structure due to the shock loading and the post-failure discontinuous behavior of the structure.

• Nuclear Engineering and Technology
• /
• v.24 no.3
• /
• pp.319-328
• /
• 1992

The Galerkin formulation of the finite element method is applied to the integral law of the first-order form of the one-group neutron transport equation in one-dimensional spherical geometry. Piecewise linear or quadratic Lagrange polynomials are utilized in the integral law for the angular flux to establish a set of linear algebraic equations. Numerical analyses are performed for the scalar flux distribution in a heterogeneous sphere as well as for the criticality problem in a uniform sphere. For the criticality problems in the uniform sphere, the results of the finite element method, with the use of continuous finite elements in space and angle, are compared with the exact solutions. In the heterogeneous problem, the scalar flux distribution obtained by using discontinuous angular and spatical finite elements is in good agreement with that from the ANISN code calculation.

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

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.198-198
• /
• 2011

The numerical simulation of dam break problem suffers from several challenges in terms of accuracy, stability, and versatility of the simulation algorithm since the water flow is generally discontinuous and presents abrupt variations. Thus, to obtain stable and accurate solutions, flow models for this purpose require numerical schemes provided with shock-capturing properties, and with the ability to work with flexible two-dimensional meshes. In this context, SU/PG method(Hughes and Brooks, 1979) is excellent candidate for the solution of the dam break problem. The weak formulation of the equations and the discontinuous polynomial basis lead to an accurate representation of bore waves(shocks). Furthermore, the discretization of the domain in finite elements is extremely effective in modeling complex geometries. In this study, a finite element model based on the SU/PG scheme is developed to solve shallow water equations and the model is applied to dam break problem. It is found that the present model accurately captures the bore wave that propagates downstream while spreading laterally and the depression wave that moves upstream. Furthermore, the propagation and formation of water surface profile compared favorably with those obtained by the previously published results.

• Nuclear Engineering and Technology
• /
• v.46 no.6
• /
• pp.783-796
• /
• 2014

The coarse mesh finite difference (CMFD) method is applied to the discontinuous finite element method based discrete ordinate calculation for source convergence acceleration. The three-dimensional (3-D) DFEM-Sn code FEDONA is developed for general geometry applications as a framework for the CMFD implementation. Detailed methods for applying the CMFD acceleration are established, such as the method to acquire the coarse mesh flux and current by combining unstructured tetrahedron elements to rectangular coarse mesh geometry, and the alternating calculation method to exchange the updated flux information between the CMFD and DFEM-Sn. The partial current based CMFD (p-CMFD) is also implemented for comparison of the acceleration performance. The modified p-CMFD method is proposed to correct the weakness of the original p-CMFD formulation. The performance of CMFD acceleration is examined first for simple two-dimensional multigroup problems to investigate the effect of the problem and coarse mesh sizes. It is shown that smaller coarse meshes are more effective in the CMFD acceleration and the modified p-CMFD has similar effectiveness as the standard CMFD. The effectiveness of CMFD acceleration is then assessed for three-dimensional benchmark problems such as the IAEA (International Atomic Energy Agency) and C5G7MOX problems. It is demonstrated that a sufficiently converged solution is obtained within 7 outer iterations which would require 175 iterations with the normal DFEM-Sn calculations for the IAEA problem. It is claimed that the CMFD accelerated DFEM-Sn method can be effectively used in the practical eigenvalue calculations involving general geometries.

• Journal of the Korean Mathematical Society
• /
• v.57 no.5
• /
• pp.1239-1266
• /
• 2020

In this paper, some novel discrete formulations for stabilizing the mixed finite element method Q1-Q0 (bilinear velocity and constant pressure approximations) are introduced and discussed for the generalized Stokes problem. These are based on stabilizing discontinuous pressure approximations via local jump operators. The developing idea consists in a reduction of terms in the local jump formulation, introduced earlier, in such a way that stability and convergence properties are preserved. The computer implementation aspects and numerical evaluation of these stabilized discrete formulations are also considered. For illustrating the numerical performance of the proposed approaches and comparing the three versions of the local jump methods alongside with the global jump setting, some obtained results for two test generalized Stokes problems are presented. Numerical tests confirm the stability and accuracy characteristics of the resulting approximations.

• Journal of the Korean Society of Safety
• /
• v.6 no.4
• /
• pp.45-52
• /
• 1991

The purpose of a fatigue crack arrest desing is to prvent a fatigue fracture of machine and structure resulted from unstable crack growth. In all cases of load transfer to second elements such as stringers, doublers or flangers, crack arrest is possible; arrest occuring when the fatigue crack reaches the second element. In the present work, the possibility of crack arrest and the design criterion of fatigue crack arrest in the variable thickness plates are examined numericaiiy by using fatigue crack arrest thresthod $\Delta$K$_{th}$of SA-508 reactor vessel steel and stress intensity factor which was obtained in the previous work as a result of 3-dimensional finite element analysis for CT type variable thickness plates having discontinuous interface.e.

• Structural Engineering and Mechanics
• /
• v.63 no.3
• /
• pp.329-345
• /
• 2017

This paper devoted to study dynamic interaction between crack and inclusion or void by developing the eXtended Finite Element Methods (XFEM). A novel XFEM approximation is presented for these structures containing multi discontinuities (void, inclusion, and crack). The level set methods are used so that elements that include a crack segment, the boundary of a void, or the boundary of an inclusion are not required to conform to discontinuous edges. The investigation covers the effects of a single circular or elliptical void / stiff inclusion, and multi stiff inclusions on the crack propagation path under dynamic loads. Both the void and the inclusion have a significant effect on the dynamic crack propagation path. The crack initially curves towards into the void, then, the crack moves round the void and propagates away the void. If a large void lies in front of crack tip, the crack may propagate into the void. If an enough small void lies in front of crack tip, the void may have a slight or no influence on the crack propagation path. For a stiff inclusion, the crack initially propagates away the inclusion, then, after the crack moves round the inclusion, it starts to propagate along its original path. As ${\delta}$ (the ratio of the elastic modulus of the inclusion to that of the matrix) increases, a larger curvature of the crack path deflection can be observed. However, as ${\delta}$ increases from 2 to 10, the curvature has an evident increase. By comparison, the curvature has a slight increase, as ${\delta}$ increases from 10 to 1000.