• Structural Engineering and Mechanics
• /
• 제28권2호
• /
• pp.129-152
• /
• 2008

This paper considers the solution of the stochastic differential equations (SDEs) with random operator and/or random excitation using the spectral SFEM. The random system parameters (involved in the operator) and the random excitations are modeled as second order stochastic processes defined only by their means and covariance functions. All random fields dealt with in this paper are continuous and do not have known explicit forms dependent on the spatial dimension. This fact makes the usage of the finite element (FE) analysis be difficult. Relying on the spectral properties of the covariance function, the Karhunen-Loeve expansion is used to represent these processes to overcome this difficulty. Then, a spectral approximation for the stochastic response (solution) of the SDE is obtained based on the implementation of the concept of generalized inverse defined by the Neumann expansion. This leads to an explicit expression for the solution process as a multivariate polynomial functional of a set of uncorrelated random variables that enables us to compute the statistical moments of the solution vector. To check the validity of this method, two applications are introduced which are, randomly loaded simply supported reinforced concrete beam and reinforced concrete cantilever beam with random bending rigidity. Finally, a more general application, randomly loaded simply supported reinforced concrete beam with random bending rigidity, is presented to illustrate the method.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2003년도 봄 학술발표회 논문집
• /
• pp.3-10
• /
• 2003

The use of frequency-dependent dynamic stiffness matrix (or spectral element matrix) in structural dynamics may provide very accurate solutions, while it reduces the number of degrees-of-freedom to improve the computational efficiency and cost problems. Thus, this paper develops a spectral element model for the thin plates moving with constant speed under uniform in-plane tension and gravity. The concept of Kantorovich method and the principle of virtual displacement is used in the frequency-domain to formulate the dynamic stiffness matrix. The present spectral element model is evaluated by comparing its solutions with the exact analytical solutions. The effects of moving speed, in-plane tension and gravity on the natural frequencies of the plate are numerically investigated.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제19권2호
• /
• pp.123-135
• /
• 2015

In this paper, a reduced-order modeling(ROM) of Burgers equations is studied based on pseudo-spectral collocation method. A ROM basis is obtained by the proper orthogonal decomposition(POD). Crank-Nicolson scheme is applied in time discretization and the pseudo-spectral element collocation method is adopted to solve linearlized equation based on the Newton method in spatial discretization. We deliver POD-based algorithm and present some numerical experiments to show the efficiency of our proposed method.

• Journal of applied mathematics & informatics
• /
• 제38권3_4호
• /
• pp.311-334
• /
• 2020

The solution of the elliptic partial differential equation has interface singularity at the points which are either the intersections of interfaces or the intersections of interfaces with the boundary of the domain. The singularities that arises in the elliptic interface problems are very complex. In this article we propose an exponentially accurate nonconforming spectral element method for these problems based on [7, 18]. A geometric mesh is used in the neighbourhood of the singularities and the auxiliary map of the form z = ln ξ is introduced to remove the singularities. The method is essentially a least-squares method and the solution can be obtained by solving the normal equations using the preconditioned conjugate gradient method (PCGM) without computing the mass and stiffness matrices. Numerical examples are presented to show the exponential accuracy of the method.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 1997년도 추계학술대회논문집; 한국과학기술회관; 6 Nov. 1997
• /
• pp.402-407
• /
• 1997

For cylindrical shells, the closed-form solutions are limited only to the cases with special boundary and/or loading conditions. Though the finite element method is certainly a powerful solution approach for the general structural dynamics problems, it is known to provide reliable solutions only in the low frequency region due to the inherent high sensitivities of structural and numerical modeling errors. Instead, the spectral element method has been proved to provide extremely accurate dynamic responses even in the high frequency region. Since the wave characteristics of a cylindrical shell becomes identical to that of a flat plate as the frequency increases, an equivalent plate model (EPM) representing the high-frequency dynamic characteristics of a cylindrical shell is introduced herein. The EPM-based spectral element analysis solutions are compared with the known analytical solutions for the corresponding cylindrical shell to confirm the validity of the present modeling approach.

• Geomechanics and Engineering
• /
• 제12권1호
• /
• pp.161-183
• /
• 2017

This paper investigates the damage identification of the concrete pile element through axial wave propagation technique using computational and experimental studies. Now-a-days, concrete pile foundations are often common in all engineering structures and their safety is significant for preventing the failure. Damage detection and estimation in a sub-structure is challenging as the visual picture of the sub-structure and its condition is not well known and the state of the structure or foundation can be inferred only through its static and dynamic response. The concept of wave propagation involves dynamic impedance and whenever a wave encounters a changing impedance (due to loss of stiffness), a reflecting wave is generated with the total strain energy forked as reflected as well as refracted portions. Among many frequency domain methods, the Spectral Finite Element method (SFEM) has been found suitable for analysis of wave propagation in real engineering structures as the formulation is based on dynamic equilibrium under harmonic steady state excitation. The feasibility of the axial wave propagation technique is studied through numerical simulations using Elementary rod theory and higher order Love rod theory under SFEM and ABAQUS dynamic explicit analysis with experimental validation exercise. Towards simulating the damage scenario in a pile element, dis-continuity (impedance mismatch) is induced by varying its cross-sectional area along its length. Both experimental and computational investigations are performed under pulse-echo and pitch-catch configuration methods. Analytical and experimental results are in good agreement.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제21권3호
• /
• pp.109-124
• /
• 2017

The numerical solution of the Stokes equation with discontinuous viscosity and singular force term is challenging, due to the discontinuity of pressure, non-smoothness of velocity, and coupled discontinuities along interface.In this paper, we give an efficient algorithm to solve this problem by employing spectral Legendre and Chebyshev approximations.First, we present the algorithm for a problem defined in rectangular domain with straight line interface. Then it is generalized to a domain with smooth curve boundary and interface by employing spectral element method. Numerical experiments demonstrate the accuracy and efficiency of our algorithm and its spectral convergence.

• 한국해양공학회:학술대회논문집
• /
• 한국해양공학회 2006년 창립20주년기념 정기학술대회 및 국제워크샵
• /
• pp.113-120
• /
• 2006

In this paper, generation mechanisms of ocean freak waves are briefly introduced in the context of wave-current interaction phenomena. As an accurate and efficient numerical tool, the spectral element method is presented with general features and specific treatment for the wave-current interaction problem. The present model of the fluid motion is based on the Navier-Stokes equations incorporating a velocity-pressure formulation. In order to deal with the free surface motion, an Arbitrary Lagrangian-Eulerian (ALE) description is adopted. As an intermediate stage of development, solution procedure and characteristic aspects of the present modeling and numerical method features are addressed in detail, and numerical results for wave-current interaction is left as further study.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 1997년도 추계학술대회논문집; 한국과학기술회관; 6 Nov. 1997
• /
• pp.160-165
• /
• 1997

Structural boundary condition is very important as a part of a structural system because it determines the dynamic characteristics of the structure. It is often to experience that experimental measurements of structural dynamic characteristics are somewhat different from the analytical predictions in which idealized boundary conditions are usually assumed. However, real structural boundary conditions are not so ideal; not perfectly clamped, for instance. Thus this paper introduces a new method to determine the non-ideal structural boundary conditions in the frequency domain. In this method, structural boundary conditions are modeled by both extensional (vertical) and torsional elastic springs. The effective springs are then determined from experimental FRFs (frequency response functions) by using the spectral element method (SEM). For a cantilevered beam experiments are conducted to determine the real boundary conditions in terms of effective springs. Dynamic characteristics (analytically predicted) based on identified boundary conditions are found to be much closer to experimental measurements when compared with those based on ideal boundary conditions.

• 호남수학학술지
• /
• 제37권3호
• /
• pp.299-315
• /
• 2015

In this paper, we study the first-order system least-squares (FOSLS) spectral method for parabolic partial differential equations. There were lots of least-squares approaches to solve elliptic partial differential equations using finite element approximation. Also, some approaches using spectral methods have been studied in recent. In order to solve the parabolic partial differential equations in parallel, we consider a parallel numerical method based on a hybrid method of the frequency-domain method and first-order system least-squares method. First, we transform the parabolic problem in the space-time domain to the elliptic problems in the space-frequency domain. Second, we solve each elliptic problem in parallel for some frequencies using the first-order system least-squares method. And then we take the discrete inverse Fourier transforms in order to obtain the approximate solution in the space-time domain. We will introduce such a hybrid method and then present a numerical experiment.