• Computational Structural Engineering
• /
• v.7 no.4
• /
• pp.85-101
• /
• 1994

A family of variational principles governing the dynamics of laminated plate has been derived using a variationally consistent shear deformable discrete laminated plate theory with particular reference to finite element procedures. The theoretical basis for the derivation is Sandhu's generalized procedure for the variational formulation of linear coupled boundary value problem. As the bilinear mapping to write the operator matrix of the field equations in self-adjoint form, convolution product was employed. Boundary conditions, initial conditions and probable internal discontinuity were explicitly included in the governing functionals. Some interesting extensions and specializations of the general variational principle were presented, which can provide many different finite element formulations for the problem.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.5
• /
• pp.1033-1042
• /
• 1994

Dynamic behavior of flexible rectangular liquid containers is analyzed by a two-dimensional coupled boundary element-finite element method. The irrotational motion of inviscid and incompressible ideal fluid is modeled by boundary elements and the motion of structure by finite elements. A singularity free integral formulation is employed for the implementation of boundary element method. Coupling is performed by using compatibility and equilibrium conditions along the interface between the fluid and structure. The fluid-structure interaction effects are reflected into the coupled equation of motion as added fluid mass matrix and sloshing stiffness matrix. By solving the eigen-problem for the coupled equation of motion, natural frequencies and mode shapes of coupled system are obtained. The free surface sloshing motion and hydrodynamic pressure developed in a flexible rectangular container due to horizontal and vertical ground motions are computed in time domain.

• Water for future
• /
• v.21 no.3
• /
• pp.299-307
• /
• 1988

A model is developed to predict the long-term beach evolution near the long groin considering the combined effects of variation of sea level, wave refraction and diffraction. A numerical solution for this problem is solved by considering the equation as a system subject to the boundary condition for longshore transport rate. One possible method is the centered Crank-Nicolson type implicit scheme. The results which ard obtained by applying this numerical model at Songdo beach, Pohang are as follows. Owing to the approximation used in the calculation of the refraction and diffraction coefficients, the discrepancy between the predicted and actual shoreline occurs to the interior of long groin. However, the shape of shoreline at the exterier of long groins agrees well.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.3
• /
• pp.185-192
• /
• 1993

A numerical model has been used for the prediction of wave agitations in a harbor which are induced by the intrusion and transformation of incident waves. Based on linear wave theory a mild-slope equation has been used. A partial absorbing boundary condition has been used on solid boundary. Functional has been derived following Chen and Mei(l974)'s technique based on Hybrid Element Method which uses finite discretisation in the inner region and analytical solution of Helmholtz equation in the outer region. Final simultaneous equation has been solved using the Gaussian Elimination Method. Helmholtz natural period and second peak period of seiche in Donghae Harbor coincide very well with the results from numerical calculation. Computed amplification factors show good agreement, especially when the reflection coefficient on solid boundary is 0.99, with those of measurements.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.3
• /
• pp.19-30
• /
• 1990

This note describes an application of Hamiton's principle to nonlinear free-surface flow problems. Two functionals are constructed based on classical Hamilton's principle with a modification due to the presence of a free surface. As an effort towards the development of an efficient numerical scheme for our problem, we present the following three test results: i) The bounding principles of the eigenvalues for the linear dispersion relation. ii) By assuming steady solitary waves, an approximate relation between the amplitudes and the speeds of solitary waves are derived from the two functionals constructed. Their numerical results are compared with those of Longuet-Higgins & Fenton(1974). iii) The shapes and charicteristics of solitary waves are computed from two sets of functionals by varying the number of total finite elements in the fluid domain.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.19 no.6
• /
• pp.607-613
• /
• 2009

A new formulation of NDIF method to the algebraic eigenvalue problem is introduced to efficiently extract natural frequencies of arbitrarily shaped plates with the simply supported boundary condition. NDIF method, which was developed by the authors for the free vibration analysis of arbitrarily shaped membranes and plates, has the feature that it yields highly accurate natural frequencies compared with other analytical methods or numerical methods(FEM and BEM). However, NDIF method has the weak point that it needs the inefficient procedure of searching natural frequencies by plotting the values of the determinant of a system matrix in the frequency range of interest. A new formulation of NDIF method developed in the paper doesn't require the above inefficient procedure and natural frequencies can be efficiently obtained by solving the typical algebraic eigenvalue problem. Finally, the validity of the proposed method is shown in several case studies, which indicate that natural frequencies by the proposed method are very accurate compared to other exact, analytical, or numerical methods.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.2
• /
• pp.86-97
• /
• 1993

A numerical method is developed for the nonlinear motion of two-dimensional wedges and axisymmetric-forced-heaving motion using Semi-Largrangian scheme under assumption of potential flows. In two-dimensional-problem Cauchy's integral theorem is applied to calculate the complex potential and its time derivative along boundary. In three-dimensional-problem Rankine ring sources are used in a Green's theorem boundary integral formulation to salve the field equation. The solution is stepped forward numerically in time by integrating the exact kinematic and dynamic free-surface boundary condition. Numerical computations are made for the entry of a wedge with a constant velocity and for the forced harmonic heaving motion from rest. The problem of the entry of wedge compared with the calculated results of Champan[4] and Kim[11]. By Fourier transform of forces in time domain, added mass coefficient, damping coefficient, second harmonic forces are obtained and compared with Yamashita's experiment[5].

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.3
• /
• pp.331-340
• /
• 2010

The basic theory and application of new adaptive finite element algorithm have been proposed in this study including the adaptive hp-refinement strategy, and the effective method for constructing hp-approximation. The hp-adaptive finite element concept needs the integrals of Legendre shape function, nonuniform p-distribution, and suitable constraint of continuity in conjunction with irregular node connection. The continuity of hp-adaptive mesh is an important problem at the common boundary of element interface. To solve this problem, the constraint of continuity has been enforced at the common boundary using the connectivity mapping matrix. The effective method for constructing of the proposed algorithm has been developed by using hierarchical nature of the integrals of Legendre shape function. To verify the proposed algorithm, the problem of simple cantilever beam has been solved by the conventional h-refinement and p-refinement as well as the proposed hp-refinement. The result obtained by hp-refinement approach shows more rapid convergence rate than those by h-refinement and p-refinement schemes. It it noted that the proposed algorithm may be implemented efficiently in practice.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.1 no.4
• /
• pp.59-68
• /
• 1997

In analysis of underground structures, the effects of artificial boundary conditions are considered as one of the major reasons for differences from experimental results. These phenomena can be overcome by using the boundary elements which satisfy the multi-layered half space conditions. The fundamental solutions of multi-layered half-space for boundary element method is formulated satisfying the transmission and reflection of waves at each layer interface and radiation conditions at bottom layer. The governing equations can be obtained from the displacements at each layer which are expressed in terms of harmonic functions. All types of waves can be included using the complete response from semi-infinite integrals with respect to horizontal wavenumbers using expansion of Fourier series and Hankel transformation. Two dimensional Green's functions are derived from cylindrical Navier equations and potentials performing infinite integration in y-direction. In this case, it is effective to transform into two dimensional problem using semi-analytical integration and sinusoidal Bessel function. Some verifications are given to show the accuracy and efficiency of the developed method, and numerical examples to demonstrate the dynamic behavior of underground with various properties.

• Journal of IKEEE
• /
• v.2 no.1 s.2
• /
• pp.160-165
• /
• 1998

This paper describes a method for field analysis inside the flat-type brushless DC motor using 2-D field simulator. Rigorous field analysis entail 3-D analysis. However, this analysis is not often appropriate for system designs because of the time and cost involved. For field analysis in this study, the 3-D problem is reduced to a 2-D boundary value problem by introducing a cylindrical cutting plane at the mean radius of the magnets. Independent of sizes and shapes of systems, the exact 2-D field results can be obtained with reasonable predictability.