• Coupled systems mechanics
• /
• 제12권4호
• /
• pp.365-389
• /
• 2023

Dynamic analysis of a typical concrete gravity dam-reservoir system is formulated by FE-(FE-TE) approach (i.e., Finite Element-(Finite Element-Truncation Element)). In this technique, dam and reservoir are discretized by plane solid and fluid finite elements. Moreover, the H-W (i.e., Hagstrom-Warburton) high-order condition imposed at the reservoir truncation boundary. This task is formulated by employing a truncation element at that boundary. It is emphasized that reservoir far-field is excluded from the discretized model. The formulation is initially reviewed which was originally proposed in a previous study. Thereafter, the response of Pine Flat dam-reservoir system is studied due to horizontal and vertical ground motions for two types of reservoir bottom conditions of full reflective and absorptive. It should be emphasized that study is carried out under high order of H-W condition applied on the truncation boundary. The initial part of study is focused on the time harmonic analysis. In this part, it is possible to compare the transfer functions against corresponding responses obtained by FE-(FE-HE) approach (referred to as exact method). Subsequently, the transient analysis is carried out. In that part, it is only possible to compare the results for low and high normalized reservoir length cases. Therefore, the sensitivity of results is controlled due to normalized reservoir length values.

• Journal of Mechanical Science and Technology
• /
• 제15권5호
• /
• pp.671-680
• /
• 2001

Boundary condition is one of the major factors to influence the numerical stability and solution accuracy in numerical analysis. One of the most important physical boundary conditions in the flowfield analysis is the wall boundary condition imposed on the body surface. To solve a two-dimensional Euler equation, totally four numerical wall boundary conditions should be prescribed. Two of them are supplied by the flow tangency condition. The other two conditions, therefore, should be prepared additionally in a suitable way. In this paper, four different sets of wall boundary conditions are proposed and then applied to solve high-speed flowfields around a quarter circle geometry. A two-dimensional compressible Euler solver is prepared based on the finite volume method. This solver hires three different upwind schemes; Steger-Warmings flux vector splitting, Roes flux difference splitting, and Lious advection upstream splitting method. It is found that the way to specify the additional numerical wall boundary conditions strongly affects the overall stability and accuracy of the upwind schemes in high-speed flow calculation. The optimal wall boundary conditions should be also chosen very carefully depending on the numerical schemes used to solve the problem.

• Coupled systems mechanics
• /
• 제1권2호
• /
• pp.183-204
• /
• 2012

A method for the analysis of wave scattering in 3-D elastic full space is developed by means of the coupled boundary-volume integral equation, which takes into account the effects of both the boundary of inclusions and the uctuation of the wave field. The wavenumber domain formulation is used to construct the Krylov subspace by means of FFT. In order to achieve the wavenumber domain formulation, the boundary-volume integral equation is transformed into the volume integral equation. The formulation is also focused on this transform and its numerical implementation. Several numerical results clarify the accuracy and effectiveness of the present method for scattering analysis.

• 한국농공학회논문집
• /
• 제46권5호
• /
• pp.61-68
• /
• 2004

Methodologies of the finite element and boundary element are combined to achieve an efficient and accurate analysis model of frame structure containing shear wall. This model analyzes the frame by employing the finite element method and the shear wall by boundary element method. This study is applicable to a specific situation, where the boundary element is surrounded by finite elements. By employing FE dominant method in which boundary stiffness matrix is transformed into finite element stiffness matrix, boundary element and finite element method are combined to analyze frame structure with walls.

• Structural Engineering and Mechanics
• /
• 제12권1호
• /
• pp.101-118
• /
• 2001

The reliable prediction of elastic vibrations of wetted complex structures, as ships, tanks, offshore structures, propulsion components etc. represent a theoretical and numerical demanding task due to fluid-structure interaction. The paper presented is addressed to the vibration analysis by a combined FE-BE-procedure based on the added mass concept utilizing a direct boundary integral formulation of the potential fluid problem in interior and exterior domains. The discretization is realized by boundary element collocation method using conventional as well as infinite boundary element formulation with analytical integration scheme. Particular attention is devoted to modelling of interior problems with both several separate or communicating fluid domains as well as thin-walled structures wetted on both sides. To deal with this specific kind of interaction problems so-called "virtual" boundary elements in areas of cut outs are placed to satisfy the kinematical conditions in partial connected fluid domains existing in realistic tank systems. Numerical results of various theoretical and practical examples demonstrate the performance of the BE-methodology presented.

• 대한기계학회논문집B
• /
• 제25권6호
• /
• pp.829-836
• /
• 2001

A dependable boundary reconstruction technique is proposed. The finite element method is used for the analysis of the direct heat conduction problem to realize the deformable grid system. An appropriate strategy for grid update is suggested. A complete sensitivity analysis is performed to obtain the derivatives required for restoration of the optimal boundary. With the result of the sensitivity analysis, the unknown boundary is sought using the sequential quadratic programming. The method is applied to reconstruction of boundaries with sinusoidal, step, and cavity form. The overall performance of the proposed method is examined by comparison between the estimated the exact boundaries.

• 한국전산유체공학회지
• /
• 제20권4호
• /
• pp.44-50
• /
• 2015

The computational efficiency of flow simulations with Multi-resolution analysis (MRA) was enhanced via the boundary treatment of the computational domain. In MRA, an adaptive dataset to a solution is constructed through data decomposition with interpolating polynomial and thresholding. During the decomposition process, the basis points of interpolation should exceed the boundary of the computational domain. In order to resolve this problem, the weight coefficients of interpolating polynomial were adjusted near the boundaries. By this boundary treatment, the computational efficiency of MRA was enhanced while the numerical accuracy of a solution was unchanged. This modified MRA was applied to two-dimensional steady Euler equations and the enhancement of computational efficiency and the maintenance of numerical accuracy were assessed.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 2007년도 정기 학술대회 논문집
• /
• pp.259-264
• /
• 2007

In this paper, to analyze the initial valued non-homogeneous elastic half space by the scaled boundary analysis, the infinite element approach was introduced. The free surface of the initial valued non-homogeneous elastic half space was mode1ed as a circumferential direction of boundary scaled boundary coordinate. The infinite element was used to represent the infinite length of the free surface. The initial value of material property(elastic modulus) was considered by the combination of the position of the sealing center and the power function of the radial direction. By use of the mapping type infinite element, the consistent e1ements formulation could be available. The performance and the feasibility of proposed approach are examined by two numerical examples.

• Structural Engineering and Mechanics
• /
• 제39권5호
• /
• pp.669-682
• /
• 2011

A Chebyshev spectral method (CSM) for the dynamic analysis of non-uniform Timoshenko beams under various boundary conditions and concentrated masses at their ends is proposed. The matrix-based Chebyshev spectral approach was used to construct the spectral differentiation matrix of the governing differential operator and its boundary conditions. A matrix condensation approach is crucially presented to impose boundary conditions involving the homogeneous Cauchy conditions and boundary conditions containing eigenvalues. By taking advantage of the standard powerful algorithms for solving matrix eigenvalue and generalized eigenvalue problems that are embodied in the MATLAB commands, chebfun and eigs, the modal parameters of non-uniform Timoshenko beams under various boundary conditions can be obtained from the eigensolutions of the corresponding linear differential operators. Some numerical examples are presented to compare the results herein with those obtained elsewhere, and to illustrate the accuracy and effectiveness of this method.

• 한국항해항만학회지
• /
• 제36권7호
• /
• pp.561-569
• /
• 2012

In this paper, the boundary element method is applied to solve the diffraction of waves by multiple vertical cylinders under the assumption of linear wave theory. A numerical analysis by boundary element method is based on Green's theorem and introduced to an integral equation for the fluid velocity potential around the cylinders. The numerical results obtained in this study are compared with the experimental data and the results of the theory using multiple scattering techniques. The comparisons show strong agreement. This numerical analysis method developed by using boundary element method could be used broadly for the design of various offshore structures to be constructed in coastal zones in the future.