• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.209-212
• /
• 2002

Since Berkhoff proposed the mild-slope equation in 1972, it has widely been used for calculation of shallow water wave transformation. Recently, it was extended to give an extended mild-slope equation, which includes the bottom slope squared term and bottom curvature term so as to be capable of modeling wave transformation on rapidly varying topography. These equations were derived by integrating the Laplace equation vertically. In the present study, we develop a finite element model to solve the Laplace equation directly while keeping the same computational efficiency as the mild-slope equation. This model assumes the vertical variation of wave potential as a cosine hyperbolic function as done in the derivation of the mild-slope equation, and the Galerkin method is used to discretize . The computational domain was discretized with proper finite elements, while the radiation condition at infinity was treated by introducing the concept of an infinite element. The upper boundary condition can be either free surface or a solid structure. The applicability of the developed model was verified through example analyses of two-dimensional wave reflection and transmission. .

• 전기의세계
• /
• 제29권4호
• /
• pp.247-255
• /
• 1980

In this study, an application of Finite Element Method which, in principle, based on variational calculus has been presented for the two-dimensional analysis of magnetic flux distribution in the shell type core of single phase transformer. The necessary stationarity condition of energy functional and boundary conditions were determined under the assumptions that the electromagnetic field considered is stationary and that the effect of eddy current is negligible. In the process of application the domain of magnetic field was divided into triangle subsectional elements and then the matrix equations were constructed for the respective triangular element and for those of all after the manipulation of minimization process to the vector potential of magnetic field at the each vertex of the element. Furthermore the numerical computation for the equations was guided by the Gaussian Elimination Methods. As the results obtained, it is found that the aspect of magnetic flux distribution inside the core as well as the leakage flux profile at the vicinity of the inner leg of the core is not much different from the well-known distribution profile of magnetic flux, however, the procedure shows to possess the merit of the uniquely deterministic nature for the flux distribution at the desired points.

• 한국지진공학회:학술대회논문집
• /
• 한국지진공학회 2002년도 춘계 학술발표회 논문집
• /
• pp.117-124
• /
• 2002

The major purpose of this study is to determine the dynamic behavior of soil-pile-structure interaction system considering the underground cavity. For the analysis, a numerical method fur ground response analysis using FE-BE coupling method is developed. The total system is divided into two parts so called far field and near field. The far field is modeled by boundary element formulation using the multi-layered dynamic fundamental solution that satisfied radiational condition of wave. And this is coupled with near field modeled by finite elements. For the verification of dynamic analysis in the frequency domain, both forced vibration analysis and free-field response analysis are performed. The behavior of soil non-linearity is considered using the equivalent linear approximation method. As a result, it is shown that the developed method can be an efficient numerical method to solve the seismic response analysis considering the underground cavity in 2D problem.

• 한국소음진동공학회논문집
• /
• 제23권11호
• /
• pp.973-981
• /
• 2013

The main role of an acoustic diffuser is to diffuse reflected sound field spatially. Since the pioneering work of Schroeder, there have been investigations to improve its performance by using shape/sizing optimization methods. In this paper, a gradient-based topology optimization algorithm is newly presented to find the optimal distribution of reflecting materials for maximizing diffuser performance. Time-harmonic acoustic analysis in a two-dimensional acoustic domain is carried out where the domain is discretized by finite elements. Perfectly matched layers are placed to surround the domain to simulate non-reflecting boundary conditions. Design variables are assigned to each element of which material properties are interpolated between those of air and those of a rigid body. An approach to extract the reflected field from the total acoustic field is employed. To validate the effectiveness of the proposed method, design problems are solved at different frequencies. The performance of the optimized diffusers obtained by the proposed method is compared against that of the conventional Schroeder diffusers.

• 한국정밀공학회:학술대회논문집
• /
• 한국정밀공학회 1996년도 추계학술대회 논문집
• /
• pp.548-554
• /
• 1996

A new method to locate modular fixtures using an optimization technique is proposed. The optimal fixture arrangement is derived to minimize the elastic deformation of a workpiece. That is, a fixture arrangement is regarded better if it minimized the elastic deformation of the workpiece while fixing a workpart of course. In this approach, the workpiece is projected into two dimensional domain to simplify the 3-dimensional fixture arrangement problem into 2-dimensional one. Thus the problem is reduced to find the optimal positions of one horizontal clamp and three locators which minimize the total deformation of the workpiece and the design variables are the location of the contact points between the boundary of the workpiece and the 4-fixels. The Genetic Algorithm is used for the optimization by mapping each design variables to a gene of a chromosome. The fitness value is the total strain energy of the workpiece calculated by the fin element analysis.

• 대한기계학회논문집A
• /
• 제26권4호
• /
• pp.775-786
• /
• 2002

A Mixed Volume and Boundary Integral Equation Method is applied for the effective analysis of elastic wave scattering problems and plane elastostatic problems in unbounded solids containing general anisotropic inclusions and voids or isotropic inclusions. It should be noted that this newly developed numerical method does not require the Green's function for anisotropic inclusions to solve this class of problems since only Green's function for the unbounded isotropic matrix is involved in their formulation for the analysis. This new method can also be applied to general two-dimensional elastodynamic and elastostatic problems with arbitrary shapes and number of anisotropic inclusions and voids or isotropic inclusions. In the formulation of this method, the continuity condition at each interface is automatically satisfied, and in contrast to finite element methods, where the full domain needs to be discretized, this method requires discretization of the inclusions only. Finally, this method takes full advantage of the pre- and post-processing capabilities developed in FEM and BIEM. Through the analysis of plane elastostatic problems in unbounded isotropic matrix with orthotropic inclusions and voids or isotropic inclusions, and the analysis of plane wave scattering problems in unbounded isotropic matrix with isotropic inclusions and voids, it will be established that this new method is very accurate and effective for solving plane wave scattering problems and plane elastic problems in unbounded solids containing general anisotropic inclusions and voids/cracks or isotropic inclusions.

• Coupled systems mechanics
• /
• 제3권3호
• /
• pp.267-289
• /
• 2014

A proper physical modeling of infilled building frame-foundation beam-soil mass interaction system is needed to predict more realistic and accurate structural behavior under static vertical loading. This is achieved via finite element method considering the superstructure, foundation and soil mass as a single integral compatible structural unit. The physical modelling is achieved via use of finite element method, which requires the use of variety of isoparametric elements with different degrees of freedom. The unbounded domain of the soil mass has been discretized with coupled finite-infinite elements to achieve computational economy. The nonlinearity of soil mass plays an important role in the redistribution of forces in the superstructure. The nonlinear behaviour of the soil mass is modeled using hyperbolic model. The incremental-iterative nonlinear solution algorithm has been adopted for carrying out the nonlinear elastic interaction analysis of a two-bay two-storey infilled building frame. The frame and the infill have been considered to behave in linear elastic manner, whereas the subsoil in nonlinear elastic manner. In this paper, the computational methodology adopted for nonlinear soil-structure interaction analysis of infilled frame-foundation-soil system has been presented.

• Korea-Australia Rheology Journal
• /
• 제14권4호
• /
• pp.161-174
• /
• 2002

A new technique for numerical calculation of viscoelastic flow based on the combination of Neural Net-works (NN) and Brownian Dynamics simulation or Stochastic Simulation Technique (SST) is presented in this paper. This method uses a "universal approximator" based on neural network methodology in combination with the kinetic theory of polymeric liquid in which the stress is computed from the molecular configuration rather than from closed form constitutive equations. Thus the new method obviates not only the need for a rheological constitutive equation to describe the fluid (as in the original Calculation Of Non-Newtonian Flows: Finite Elements St Stochastic Simulation Techniques (CONNFFESSIT) idea) but also any kind of finite element-type discretisation of the domain and its boundary for numerical solution of the governing PDE's. As an illustration of the method, the time development of the planar Couette flow is studied for two molecular kinetic models with finite extensibility, namely the Finitely Extensible Nonlinear Elastic (FENE) and FENE-Peterlin (FENE-P) models.P) models.

• Earthquakes and Structures
• /
• 제5권2호
• /
• pp.161-205
• /
• 2013

We study seismically induced, anti-plane strain wave motion in a non-homogeneous geological region containing tunnels. Two different scenarios are considered: (a) The first models two tunnels in a finite geological region embedded within a laterally inhomogeneous, layered geological profile containing a seismic source. For this case, labelled as the first boundary-value problem (BVP 1), an efficient hybrid technique comprising the finite difference method (FDM) and the boundary element method (BEM) is developed and applied. Since the later method is based on the frequency-dependent fundamental solution of elastodynamics, the hybrid technique is defined in the frequency domain. Then, an inverse fast Fourier transformation (FFT) is used to recover time histories; (b) The second models a finite region with two tunnels, is embedded in a homogeneous half-plane, and is subjected to incident, time-harmonic SH-waves. This case, labelled as the second boundary-value problem (BVP 2), considers complex soil properties such as anisotropy, continuous inhomogeneity and poroelasticity. The computational approach is now the BEM alone, since solution of the surrounding half plane by the FDM is unnecessary. In sum, the hybrid FDM-BEM technique is able to quantify dependence of the signals that develop at the free surface to the following key parameters: seismic source properties and heterogeneous structure of the wave path (the FDM component) and near-surface geological deposits containing discontinuities in the form of tunnels (the BEM component). Finally, the hybrid technique is used for evaluating the seismic wave field that develops within a key geological cross-section of the Metro construction project in Thessaloniki, Greece, which includes the important Roman-era historical monument of Rotunda dating from the 3rd century A.D.

• 전산구조공학
• /
• 제7권4호
• /
• pp.127-136
• /
• 1994

2차원과 3차원모델에 대해 진동전파구간에 설치된 방진구의 진동차단효과가 주파수영역의 경계요소법에 의해 반무한영역지반과 충진지반에서 검토되었다. 연구결과는, 방진구의 효율이 진동원과 관련된 방진구의 위치 및 깊이에 큰 영향을 받으며, 그들의 최적크기는 진동주파수에 대해 민감한 것을 나타낸다.