• Journal of Adhesion and Interface
• /
• v.6 no.1
• /
• pp.1-6
• /
• 2005

The high accurate simulation of very thin glue layers based on the finite element method is still connected to many problems which result from the necessity to construct a complicated mesh of essentially different sizes of elements. This can lead to a loss of accuracy, unstable calculations and even loss of convergence. However, the implementation of special transmission elements along the glue ling and special edge-elements in the near-edge region would lead to a dramatic decrease of number of finite elements in the mesh and thus, prevent unsatisfactory phenomena in numerical analysis and extensive computation time. The theoretical basis for such special elements is the knowledge about appropriate transmission conditions and the edge effects near the free boundary of the adhesive layer. Therefore, recently proposed so-called non-classical transmission conditions and the behavior near the free edge are investigated in the context of the single-lap tensile-shear test of adhesive technology.

• Journal of the Korean Society for Precision Engineering
• /
• v.23 no.6 s.183
• /
• pp.128-135
• /
• 2006

Rational dynamic modeling and analysis method f3r complex structures are studied with special attention to slide way joints. For modeling of slide way joints, a general modeling technique is used by using the influence coefficients method which is applied to the conversion of detailed finite element model to the equivalent reduced joint model. The theoretical part of this method is illustrated and the method is applied to the structure with slide way joint. In this method, the non-linearity of the contact surfaces is considered within a proper range and the boundary effect of the joint model could be eliminated. The proposed method was applied to finite element modal analysis of a clamp jointed cantilever beam and slide way joints of the vertical type lathe. The method can also be used to other kinds of joint modeling. The results of these analysis were compared with those of Yoshimura models and rigid joint models, which demonstrated the practical applicability of the proposed method.

• Computational Structural Engineering
• /
• v.10 no.2
• /
• pp.213-223
• /
• 1997

Calculation of elastic wave fields has important applications in a variety of engineering fields including NDE (Non-destructive evaluation). Scattering problems have been investigated by numerous authors with different solution schemes. For simple geometries of the scatterers (e.g., cylinders or spheres), the analysis of steady-state elastic wave scattering has been carried out using analytical techniques. For arbitrary geometries and multiple inclusions, numerical methods have been developed. Special finite element methods, e.g., the infinite element method and a hybrid method called the Global-Local finite element method have also been developed for this purpose. Recently, the boundary integral equation method has been used successfully to solve scattering problems. In this paper, a volume integral equation method (VIEM) is proposed as a new numerical solution scheme for the solution of general elasto-dynamic problems in unbounded solids containing multiple inclusions and voids or cracks. A boundary integral equation method (BIEM) is also presented for elastic wave scattering problems. The relative advantage of the volume and boundary integral equation methods for solving scattering problems is discussed.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.237-255
• /
• 2013

In this paper, a meshfree-enriched finite element method (ME-FEM) is introduced for the large deformation analysis of nonlinear path-dependent problems involving contact. In linear ME-FEM, the element formulation is established by introducing a meshfree convex approximation into the linear triangular element in 2D and linear tetrahedron element in 3D along with an enriched meshfree node. In nonlinear formulation, the area-weighted smoothing scheme for deformation gradient is then developed in conjunction with the meshfree-enriched element interpolation functions to yield a discrete divergence-free property at the integration points, which is essential to enhance the stress calculation in the stage of plastic deformation. A modified variational formulation using the smoothed deformation gradient is developed for path-dependent material analysis. In the industrial metal forming problems, the mortar contact algorithm is implemented in the explicit formulation. Since the meshfree-enriched element shape functions are constructed using the meshfree convex approximation, they pose the desired Kronecker-delta property at the element edge thus requires no special treatments in the enforcement of essential boundary condition as well as the contact conditions. As a result, this approach can be easily incorporated into a conventional displacement-based finite element code. Two elasto-plastic problems are studied and the numerical results indicated that ME-FEM is capable of delivering a volumetric locking-free and pressure oscillation-free solutions for the large deformation problems in metal forming analysis.

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

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1997.10a
• /
• 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.

• Journal of the Korean Society for Precision Engineering
• /
• v.22 no.5 s.170
• /
• pp.205-213
• /
• 2005

This paper addresses a method for numerical simulation in the pressing process of hot glass. Updated Lagrangian finite element formulations are employed for the flow and energy equations to accommodate moving meshes. The model is assumed axi-symmetric and creep flow is assumed due to the high viscosity. Commercial software ANSYS is used to solve the coupled flow and energy equations. Moving contact points as well as free surface during the pressing are effectively calculated and updated by utilizing API functions of CAD software Unigraphics. The mesh distortion problem near the wall is overcome by automatic remeshing, and the temperatures of the new mesh are conveniently interpolated by using a unique function of ANSYS. The developed model is applied to the pressing process of TV glasses. In conclusion, the presented method shows that the pressing process accompanying moving boundary can be simulated by effectively combining general purpose software without resorting to special dedicated codes.

• Coupled systems mechanics
• /
• v.7 no.3
• /
• pp.353-371
• /
• 2018

The Wavenumber or more accurately Wavenumber-FD approach was initially introduced for two-dimensional dynamic analysis of concrete gravity dam-reservoir systems. The technique was formulated in the context of pure finite element programming in frequency domain. Later on, a variation of the method was proposed which was referred to as Wavenumber-TD approach suitable for time domain type of analysis. Recently, it is also shown that Wavenumber-FD approach may be applied for three-dimensional dynamic analysis of concrete arch dam-reservoir systems. In the present study, application of its variation (i.e., Wavenumber-TD approach) is investigated for three-dimensional problems. The method is initially described. Subsequently, the response of idealized Morrow Point arch dam-reservoir system is obtained by this method and its special cases (i.e., two other well-known absorbing conditions) for time harmonic excitation in stream direction. All results for various considered cases are compared against the exact response for models with different values of normalized reservoir length and reservoir base/sidewalls absorptive conditions.

• Journal of Magnetics
• /
• v.18 no.2
• /
• pp.95-104
• /
• 2013

This paper presents a general analytical method for predicting the magnetic fields of different Halbach magnet arrays with or without back iron mounted on slotless permanent magnet (PM) linear machines. By using Fourier decomposition, the magnetization components of four typical Halbach magnet arrays are determined. By applying special synthetic boundary conditions on the PM surfaces, the expressions of their magnetic field distributions are derived based on the magnetic scalar potential (MSP), which are simpler than those based on the magnetic vector potential (MVP). The correctness of the method is validated by finite element analysis. The harmonics of airgap flux density waveforms of these Halbach magnet arrays with or without back iron are also compared and optimized.

• Structural Engineering and Mechanics
• /
• v.68 no.6
• /
• pp.713-720
• /
• 2018

In this paper a new numerical method to determine the elastic curve of the simply supported beams of variable cross-section is demonstrated. In general case it needs to solve linear or small nonlinear second order differential equations with prescribed boundary conditions. For numerical solution the initial values of the slope and the deflection of the end cross-section of the beam is necessary. For obtaining the initial values a lively procedure is developed: it is a special application of the shooting method because boundary value problems can be transformed into initial value problems. As a result of these transformations the initial values of the differential equations are obtained with high accuracy. Procedure is applied for calculating of elastic curve of a simply supported beam of variable cross-section. Results of these numerical procedures, analytical solution of the linearized version and finite element method are compared. It is proved that the suggested procedure yields technically accurate results.