• Proceedings of the International Microelectronics And Packaging Society Conference
• /
• 2003.11a
• /
• pp.247-250
• /
• 2003

In this paper, the boundary element analysis of viscoelastic strain energy release ,ate G(t) for the cracked linear viscoelastic materials is attempted. This study proposes the G(t) equation and the calculating method of G(t) by time-domain boundary element analysis for the viscoelastic solids. The G(t) is defined as the derivative of the viscoelastic potential energy $\Pi(t)$ with respect to crack length a. Two example problems are presented to show the applicability of the proposed method to the analysis of the cracked linear viscoelastic solids. Numerical results of example problems show the accuracy and effectiveness of the proposed method.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.22 no.3
• /
• pp.9-18
• /
• 1985

The numerical calculation for solving boundary-value problem related to potential flows with a free surface is carried out by application of the localized finite element method. Only forced motion of 2-D body in infinitely deep fluid is considered, although this schemes is equally applicable to any first order time-harmonic problems of similar nature. The infinite domain of the fluid is separated into the inner flow field and the outer flow field with common inter-surface boundary. The finite element method is applied to obtain the solution in the inner flow field and the Green functions are utilized to represent the solution in the outer flow field. At the inter-surface boundary, the continuity of the value of potential and the normal derivative of the potential(i.e. matching condition) is conserved. The present method has better computational efficiency than the previous LFEM and the integral equation method of Frank. This enhanced computational efficiency is presumably due to the fact that the present method gives a symmetric coefficient matrix and requires less computational time in calculating the influence coefficient matrix of Green function than the integral equation method. And the irregular frequency desen't exist because the uniqueness of the solution is assured by the such that the exact free surface condition is satisfied on the boundary of the localized finite element region(i.e. inner region). As an example of the above method, the hydrodynamic forces for the circular cylinder and the rectangular cylinders are calculated. In the computed results, the small number of singularity distribution segments($3{\sim}6$) give good result relative to Ursell's and Vugts'.

• Journal of Ocean Engineering and Technology
• /
• v.7 no.2
• /
• pp.9-18
• /
• 1993

Numerical wave tank is a robust tool by which the nonlinear interactions between the body and the free-surface can be treated in time-domain. In this paper, a two-dimensional numerical wave tank based on the Spectral/Boundary-Element Method is developed, and applied successfully to the study of nonlinear wave diffraction around a submerged circular cylinder. Particularly, it is shown that the high-order wave components of significant wave height are developed in the lee-side of the cylinder and that these waves result in a negative drift force on the circular cylider.

• Computational Structural Engineering
• /
• v.10 no.1
• /
• pp.129-134
• /
• 1997

Interfacial stress singularity in a unidirectional two-dimensional laminate model consisting of an elastic fiber and a viscoelastic matrix has been investigated using the time-domain boundary element method. First, the interfacial singular stresses between the fiber and the matrix of a unidirectional laminate subjected to a uniform transverse tensile strain have been investigated near the free surface, but without any defect or any edge crack. Such a stress singularity might lead to fiber-matrix debonding or interfacial edge cracks. Then, the overall stress intensity factor for the case of a small interfacial edge crack of length a has been computed.

• Journal of Ocean Engineering and Technology
• /
• v.38 no.2
• /
• pp.53-63
• /
• 2024

Submerged breakwaters can be installed underneath floating structures to reduce the external wave loads acting on the structure. The objective of this study was to establish an optimization analysis framework to determine the corresponding shape and position of the submerged breakwater that can minimize or maximize the external forces acting on the floating structure. A two-dimensional frequency-domain boundary element method (FD-BEM) based on the linear potential theory was developed to perform the hydrodynamic analysis. A metaheuristic algorithm, the advanced particle swarm optimization, was newly coupled to the FD-BEM to perform the optimization analysis. The optimization analysis process was performed by calling FD-BEM for each generation, performing a numerical analysis of the design variables of each particle, and updating the design variables using the collected results. The results of the optimization analysis showed that the height of the submerged breakwater has a significant effect on the surface piercing body and that there is a specific area and position with an optimal value. In this study, the optimal values of the shape and position of a single submerged breakwater were determined and analyzed so that the external force acting on a surface piercing body was minimum or maximum.

• 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.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.6
• /
• pp.447-457
• /
• 2018

In this study, a two-dimensional fully nonlinear transient wave numerical tank was developed using a desingularized indirect boundary integral equation method. The desingularized indirect boundary integral equation method is simpler and faster than the conventional boundary element method because special treatment is not required to compute the boundary integral. Numerical simulations were carried out in the time domain using the fourth order Runge-Kutta method. A mixed Eulerian-Lagrangian approach was adapted to reconstruct the free surface at each time step. A numerical damping zone was used to minimize the reflective wave in the downstream region. The interpolating method of a Gaussian radial basis function-type artificial neural network was used to calculate the gradient of the free surface elevation without element connectivity. The desingularized indirect boundary integral equation using an isolated point source and radial basis function has no need for information about the element connectivity and is a meshless method that is numerically more flexible. In order to validate the accuracy of the numerical wave tank based on the desingularized indirect boundary integral equation method and meshless technique, several numerical simulations were carried out. First, a comparison with numerical results according to the type of desingularized source was carried out and confirmed that continuous line sources can be replaced by simply isolated sources. In addition, a propagation simulation of a $2^{nd}$-order Stokes wave was carried out and compared with an analytical solution. Finally, simulations of propagating waves in shallow water and propagating waves over a submerged bar were also carried and compared with published data.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.18 no.11
• /
• pp.1157-1169
• /
• 2008

This paper presents spectral element formulation which approximates Lamb wave propagation by PZT transducers bonded on a thin plate. A two layer beam model under 2-D plane strain condition is introduced to simulate high-frequency dynamic responses induced by a piezoelectric (PZT) layer rigidly bonded on a base plate. Mindlin-Herrmann and Timoshenko beam theories are employed to represent the first symmetric and anti-symmetric Lamb wave modes on a base plate, respectively. The Euler-Bernoulli beam theory and 1-D linear piezoelectricity are used to model the electro-mechanical behavior of a PZT layer. The equations of motions of a two layer beam model are derived through Hamilton's principle. The necessary boundary conditions associated with the electro-mechanical properties of a PZT layer are formulated in the context of dual functions of a PZT layer as an actuator and a sensor. General spectral shape functions of response field and the associated boundary conditions are obtained through equations of motions converted into frequency domain. Detailed spectrum element formulation for composing the dynamic stiffness matrix of a two layer beam model is presented as well. The validity of the proposed spectral element is demonstrated through numerical examples.

• Journal of Ocean Engineering and Technology
• /
• v.15 no.2
• /
• pp.11-21
• /
• 2001

The focus of this paper is on the numerical investigation of obliquely incident wav interactions with a system composed of fully submerged and floating dual buoy/vertical-flexible-membrane breakwaters placed in parallel with spacing between two systems. The fully submerged two systems allow surface and bottom gaps to enable wave transmission over and under the system. The problem is formulated based on the two-dimensional multi-domain hydro-elastic linear wave-body interaction theory. The hydrodynamic interaction of oblique incident waves with the combination of the rigid and flexible bodies was solved by the distribution of the simple sources (modified Bessel function of the second kind) that satisfy the Helmholz governing equation in fluid domains. A boundary element program for three fluid domains based on a discrete membrane dynamic model and simple source distribution method is developed. Using this developed computer program, the performance of various dual systems varying buoy radiuses and drafts, membrane lengths, gaps, spacing, mooring-lines stiffness, mooring types, water depth, and wave characteristics is thoroughly examined. It is found that the fully submerged and floating dual buoy/membrane breakwaters can, if it is properly tuned to the coming waves, have good performances in reflecting the obliquely incident waves over a wide range of wave frequency and headings.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.58 no.10
• /
• pp.1936-1941
• /
• 2009

This paper deals with a new methodology for topology optimization in which the topology of the design domain may change during the shape optimization process. To achieve this, the concept of the topological gradient is introduced to compute the sensitivity of an objective function when a small hole is drilled in the domain. Based on shape and topological sensitivity values, the shape and topology of the design domain may be simultaneously changed during design iterations if necessary. To verify the advantages and also to facilitate understanding of the method itself, two electrostatic design problems have been tested by using 2D finite element analysis: the first is the inverse problem of a simple dielectric model and the second is the rotor design of a MEMS actuator.