• Proceedings of the Korea Concrete Institute Conference
• /
• 1998.10a
• /
• pp.450-455
• /
• 1998

In this paper, finite element with embedded discontinuous line is introduced in order to avoid the difficulties of adding new nodal points along with crack growth in discrete crack model. With the discontinuous element using discontinuous shape function, stiffness matrix of finite element is derived and dual mapping technique for numerical integration is employed. Using the finite element program made with employed algorithms, algorithm is verified and fracture analysis of simple concrete beam is performed.

• Computers and Concrete
• /
• v.23 no.4
• /
• pp.235-243
• /
• 2019

The discontinuous Galerkin method (DGM) has become widely used as it possesses several qualities, such as a natural ability to dealing with discontinuities. DGM has its major success related to fluid mechanics. Its major importance is the ability to deal with discontinuities and still provide high order of approximation. That is an important advantage when simulating cracking propagation. No remeshing is necessary during the propagation, since the crack path follows the interface of elements. However, DGM comes with the drawback of an increased number of degrees of freedom when compared to the classical continuous finite element method. Thus, it seems a natural approach to combine them in the same simulation obtaining the advantages of both methods. This paper proposes the application of the combined continuous-discontinuous Galerkin method (CDGM) to crack propagation. An important engineering problem is the simulation of crack propagation in concrete structures. The problem is characterized by discontinuities that evolve throughout the domain. Crack propagation is simulated using CDGM. Discontinuous elements are placed in regions with discontinuities and continuous elements elsewhere. The cohesive zone model describes the fracture process zone where softening effects are expressed by cohesive zones in the interface of elements. Two numerical examples demonstrate the capacities of CDGM. In the first example, a plain concrete beam is submitted to a three-point bending test. Numerical results are compared to experimental data from the literature. The second example deals with a full-scale ground slab, comparing the CDGM results to numerical and experimental data from the literature.

• Proceedings of the Korea Concrete Institute Conference
• /
• 1996.04a
• /
• pp.247-252
• /
• 1996

In the concrete structures, cracks occur in various causes and the cracks seriously affect the functions of structures. The analysis techniques of progressive crack in the concrete have been improved with the advance of numerical techniques. The discrete crack model used in finite element program for the analysis of progressive failure is very effective, but it can not be easily implemented into numerical procedures because of difficult handing of nodal points in finite element meshes for crack growth. This paper introduces one of the techniques which skips the difficulty. In this paper, the modeling of progressive failure using finite element formulation is explained for the analysis of concrete fracture. The discontinuous element using the discontinuous shape function and the dual mapping technique in the numerical integration are implemented into finite element code for this purpose. It is shown that developed finite element program can predict the quasi-brittle behavior of concrete including ultimate load. The comparisons of the analysis results with other data are also shown.

• Computational Structural Engineering
• /
• v.6 no.4
• /
• pp.83-88
• /
• 1993

A time-discontinuous Galerkin method based upon using a finite element formulation in time has evolved. This method, working from the differential equation viewpoint, is different from those which have been generally used. They admit discontinuities with respect to the time variable at each time step. In particular, the elements can be chosen arbitrarily at each time step with no connection with the elements corresponding to the previous step. Interpolation functions and weighting functions are taken to be discontinuous across inter-element boundaries. These methods lead to a unconditional stable higher-order accurate ordinary differential equation solver.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.5 no.6
• /
• pp.141-147
• /
• 1997

For the effective finite element analysis of the structures including material interfaces or contact surfaces, interface elements are proposed. Most of early works in this problem require not only iterative computation but also complex formulation because of the kinematic nonlinearities caused from the discontinuous behavior and the stress concentration phenomena. The proposed elements, however, are consistently formulated using relative displacements and tractions between top and bottom regular finite elements. The effectiveness of these elements are shown by solving various numerical sample problems including an leaf spring and comparing with results of general finite element analysis. As a result, more stable solutions are conveniently obtaines using interface elements than regular finite elements.

• Nuclear Engineering and Technology
• /
• v.52 no.6
• /
• pp.1137-1147
• /
• 2020

The discrete ordinates method (S_N) is one of the major shielding calculation method, which is suitable for solving deep-penetration transport problems. Our objective is to explore the available quadrature sets and to improve the accuracy in shielding problems involving strong anisotropy. The linear discontinuous finite-element (LDFE) quadrature sets based on the icosahedron (in short, ICLDFE quadrature sets) are developed by defining projected points on the surfaces of the icosahedron. Weights are then introduced in the integration of the discontinuous finite-element basis functions in the relevant angular regions. The multivariate secant method is used to optimize the discrete directions and their corresponding weights. The numerical integration of polynomials in the direction cosines and the Kobayashi benchmark are used to analyze and verify the properties of these new quadrature sets. Results show that the ICLDFE quadrature sets can exactly integrate the zero-order and first-order of the spherical harmonic functions over one-twentieth of the spherical surface. As for the Kobayashi benchmark problem, the maximum relative error between the fifth-order ICLDFE quadrature sets and references is only -0.55%. The ICLDFE quadrature sets provide better integration precision of the spherical harmonic functions in local discrete angle domains and higher accuracy for simple shielding problems.

• Structural Engineering and Mechanics
• /
• v.11 no.6
• /
• pp.591-604
• /
• 2001

In this paper, a finite element with embedded displacement discontinuity which eliminates the need for remeshing of elements in the discrete crack approach is applied for the progressive fracture analysis of concrete structures. A finite element formulation is implemented with the extension of the principle of virtual work to a continuum which contains internal displacement discontinuity. By introducing a discontinuous displacement shape function into the finite element formulation, the displacement discontinuity is obtained within an element. By applying either a nonlinear or an idealized linear softening curve representing the fracture process zone (FPZ) of concrete as a constitutive equation to the displacement discontinuity, progressive fracture analysis of concrete structures is performed. In this analysis, localized progressive fracture simultaneous with crack closure in concrete structures under mixed mode loading is simulated by adopting the unloading path in the softening curve. Several examples demonstrate the capability of the analytical technique for the progressive fracture analysis of concrete structures.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.6 no.1
• /
• pp.1-15
• /
• 2002

Advection-dominated transport problems possess difficulties in the design of numerical methods for solving them. Because of the hyperbolic nature of advective transport, many characteristic numerical methods have been developed such as the classical characteristic method, the Eulerian-Lagrangian method, the transport diffusion method, the modified method of characteristics, the operator splitting method, the Eulerian-Lagrangian localized adjoint method, the characteristic mixed method, and the Eulerian-Lagrangian mixed discontinuous method. In this paper relationships among these characteristic methods are examined. In particular, we show that these sometimes diverse methods can be given a unified formulation. This paper focuses on characteristic finite element methods. Similar examination can be presented for characteristic finite difference methods.

• Computers and Concrete
• /
• v.5 no.4
• /
• pp.401-416
• /
• 2008

The paper presents a methodology to model three-dimensional reinforced concrete members by means of embedded discontinuity elements based on the Continuum Strong Discontinuous Approach (CSDA). Mixture theory concepts are used to model reinforced concrete as a 3D composite material constituted of concrete with long fibers (rebars) bundles oriented in different directions embedded in it. The effects of the rebars are modeled by phenomenological constitutive models devised to reproduce the axial non-linear behavior, as well as the bond-slip and dowel action. The paper presents the constitutive models assumed for the components and the compatibility conditions chosen to constitute the composite. Numerical analyses of existing experimental reinforced concrete members are presented, illustrating the applicability of the proposed methodology.

• Journal of the Korean Mathematical Society
• /
• v.39 no.3
• /
• pp.363-376
• /
• 2002

Stability result is obtained for the approximation of the stationary Stokes problem with nonconforming elements proposed by Douglas et al [1] for the velocity and discontinuous piecewise constants for the pressure on qudrilateral elements. Optimal order $H^1$and $L^2$error estimates are derived.