• Structural Engineering and Mechanics
• /
• v.75 no.3
• /
• pp.311-325
• /
• 2020

In applying the spectral stochastic finite element methods to the frequency response analysis, the conventional methods are known to give unstable and inaccurate results near the natural frequencies. To address this issue, a new sensitivity based stabilized formulation for stochastic frequency response analysis is proposed in this paper. The main difference over the conventional spectral methods is that the polynomials of random variables are applied to both numerator and denominator in approximating the harmonic response solution. In order to reflect the resonance behavior of the structure, the denominator polynomials is constructed by utilizing the natural frequency sensitivity and the random mode superposition. The numerator is approximated by applying a polynomial chaos expansion, and its coefficients are obtained through the Galerkin or the spectral projection method. Through various numerical studies, it is seen that the proposed method improves accuracy, especially in the vicinities of structural natural frequencies compared to conventional spectral methods.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.13 no.4
• /
• pp.267-280
• /
• 2009

In this paper we investigate a simple two-level additive Schwarz preconditioner for the P1 symmetric interior penalty Galerkin method of the Poisson equation on rectangular meshes. The construction is based on the decomposition of the global space of piecewise linear polynomials into the sum of local subspaces, each of which corresponds to an element of the underlying mesh, and the global coarse subspace consisting of piecewise constants. This preconditioner is a direct combination of the block Jacobi iteration and the cell-centered finite difference method, and thus very easy to implement. Explicit upper and lower bounds for the maximum and minimum eigenvalues of the preconditioned matrix system are derived and confirmed by some numerical experiments.

• Water Engineering Research
• /
• v.4 no.1
• /
• pp.45-58
• /
• 2003

Two-dimensional depth-averaged advection-dispersion equation was simulated using FEM. In the straight rectangular channel, the advection-dispersion processes are simulated so that these results can be compared with analyti-cal solutions for the transverse line injection and the point injection. In the straight domain the standard Galerkin method with the linear basis function is found to be inadequate to the advection-dispersion analysis compared to the upwind finite element scheme. The experimental data in the S-curved channel were compared with the result by the numerical model using SUPG(Streamline upwind Petrov-Galerkin) method.

• The Transactions of the Korean Institute of Electrical Engineers
• /
• v.36 no.4
• /
• pp.252-258
• /
• 1987

Time periodic finite element solutions for sinusoidally excited electromagnetic field problems in moving media are presented. Solutions by the Galerkin method contain spurious oscillations when grid Peclet number is more than one. To suppress these oscillations an upwind finite element method using two different time periodic test functions is introduced. One is multiplied to second and first-order space derivative terma and the other to the time derivative term. Test functions are obtained from trial functions by adding or subtracting quadratic bias functions with appropriate scaling factors. Phase differences are considered between trial functions and bias functions. For simple interpretations of the phase differences, complex scaling factors are used. The proposed method is developed to give nodally exact solutions for uniform grid spacing in one dimensional problems. Based on the one dimensional results, a two dimensional upwinding scheme is also derived.

• 연구논문집
• /
• s.28
• /
• pp.123-135
• /
• 1998

A shape function for element free Galerkin method is carved from Shepard interpolant of singular weight and consistency condition. Thus present shape function is an interpolation and has no singularities. The shape function is applied to cantilever bending problems and gives good results in comparison with beam theory. Finally it is shown that the coupling with finite element method is made easily without any additional treaties.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.3
• /
• pp.211-236
• /
• 2019

This paper is devoted to the study and investigation of the Runge-Kutta discontinuous Galerkin method for a system of differential equations consisting of two hyperbolic conservation laws. The numerical coupling flux which is used at a given interface (x = 0) is the upwind flux. Moreover, in the linear case, we derive optimal convergence rates in the $L_2$-norm, showing an error estimate of order ${\mathcal{O}}(h^{k+1})$ in domains where the exact solution is smooth; here h is the mesh width and k is the degree of the (orthogonal Legendre) polynomial functions spanning the finite element subspace. The underlying temporal discretization scheme in time is the third-order total variation diminishing Runge-Kutta scheme. We justify the advantages of the Runge-Kutta discontinuous Galerkin method in a series of numerical examples.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2000.04b
• /
• pp.324-331
• /
• 2000

In this study, a new algorithm analyzing dynamic crack propagation problem by the coupling technique of Meshfree Method and Finite Element Method is proposed. The coupling procedure of two methods is presented with a short description of Meshfree Method especially, Element-free Galerkin (EFG) method. The elastodynamic fracture theory is presented, and a numerical implementation procedure for dynamic fracture analysis by Meshfree Method is also discussed. A couple of dynamic crack propagation problems illustrate the performance of the propsed technique. The accuracy of numerical solutions by the developed algorithm are compared with those of analytical solutions and experimental ones.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.1
• /
• pp.87-94
• /
• 2003

A geometric non-linear finite element formulation of CRTS reflector subjected to displacement loads, corresponding to the successional assembly steps of the reflector, is presented in order to determine the initial static equilibrium state based on the displacement incremental method. Parametric analyses of the influence of cables and mechanical properties of the reflector on the shape error between reference and equilibrium surfaces have been studied. These results of the present study are compared with the others using Galerkin mothod and NASS 98 program to demonstrate the feasibility.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.4 no.3
• /
• pp.53-61
• /
• 1984

A nonlinear formulation of a beam finite element(NB6) on the total Lagrangian mode for the geometrically nonlinear analysis of two-dimensional elastic framed structures is presented. The NB6 beam element has been degenerated from the three-dimensional continuum by introducing the deep beam assumptions and consists of three reference nodes and three relative nodes. The element characteristics are derived by discretizing the beam equations of motion using the Galerkin weighted residual method and are reduced-integrated repeatedly for each loading step by the Newton-Raphson iteration techpique. Several numerical examples are given to demonstrate the accuracy and versatility of the proposed nonlinear NB6 beam element.

• Advances in nano research
• /
• v.7 no.2
• /
• pp.99-108
• /
• 2019

Higher-order theories are very important to investigate the mechanical properties and behaviors of nanoscale structures. In this study, a free vibration behavior of SiNW resting on elastic foundation is investigated via Eringen's nonlocal elasticity theory. Silicon Nanowire (SiNW) is modeled as simply supported both ends and clamped-free Euler-Bernoulli beam. Pasternak two-parameter elastic foundation model is used as foundation. Finite element formulation is obtained nonlocal Euler-Bernoulli beam theory. First, shape function of the Euler-Bernoulli beam is gained and then Galerkin weighted residual method is applied to the governing equations to obtain the stiffness and mass matrices including the foundation parameters and small scale parameter. Frequency values of SiNW is examined according to foundation and small scale parameters and the results are given by tables and graphs. The effects of small scale parameter, boundary conditions, foundation parameters on frequencies are investigated.