• Journal of the Korean Society of Safety
• /
• v.19 no.1
• /
• pp.38-43
• /
• 2004

Many analysis methods, including finite element method, have been suggested and used for assessing the integrity of cracked structures. In the paper, in order to analyze arbitrary three dimensional cracks, the finite element alternating method is extended. The crack is modeled by the symmetric Galerkin boundary element method as a distribution of displacement discontinuities, which is formulated as singularity-reduced integral equations. And the finite element method is used to calculate the stress values for the uncracked body only. Applied the proposed method to several example problems for planner cracks in finite bodies, the accuracy and efficiency of the method were demonstrated.

• Structural Engineering and Mechanics
• /
• v.20 no.2
• /
• pp.225-239
• /
• 2005

Applications of the Trefftz boundary element method (BEM) to anti-plane electroelastic problems are presented in this paper. Both direct and indirect methods with domain decomposition are discussed in details. Each crack is treated as semi-infinite thin slit defined in a subregion, for which a particular solution of the anti-plane problem, satisfying exactly the crack-face condition, is derived. The stress intensity factors defined at each crack tip can be directly computed from the coefficients of the particular solution. The performance of the proposed formulation is assessed by two examples and comparison is made with results obtained by other approaches. The Trefftz boundary element approach is demonstrated to be suitable for the analysis of the anti-plane problem of piezoelectric materials.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.2
• /
• pp.553-563
• /
• 1996

An algorithm for automatic mesh generation of two-dimensional arbitrary planar domain is proposed by using Delaunay triangulation algorithm. An efficient algorithm is proposed for the construction of Delaunay triangulation algorithm over convex planar domain. From the definition of boundary, boundary nodes are first defined and then interior nodes are generated ensuring the Delaunay property. These interior nodes and the boundary nodes are then linked up together to produce a valid triangular mesh for any finite element analysis. Through the various example, it is found that high-quality triangular element meshes are obtained by Delaunay algorithm, showing the robustness of the current method. The proposed mesh generation scheme has been extended to automatic remeshing, which is applicable to FE analysis including large deformation and large distortion of elements.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.6 no.4
• /
• pp.45-51
• /
• 2003

The reflection and transmission of submerged breakwater with trapezoidal type are computed numerically using boundary element method. The analysis method is based on the wave pressure function with the contlnuit? in the analytical region including fluid and porous structures. Wane motion within the porous structures is simulated by introducing the linear dissipation coefficient and added mass coefficient. The results indicate that transmission and reflection coefficient are determined due to the change of slope of submerged breakwater with trapezoidal type.

• Computational Structural Engineering
• /
• v.8 no.2
• /
• pp.147-151
• /
• 1995

The time-domain boundary element analysis has been attempted for investigating the displacement and stress in an elastic-viscoelastic compound structure. The subdomain technique has been employed and the whole domain has been divided into two subdomains, which are respectively a homogeneous elastic zone and a homogeneous viscoelastic zone. The boundary element equation has been formulated using the equilibrium and continuity conditions at the common interface. The numerical results of example problem has been presented.

• Transactions of Materials Processing
• /
• v.11 no.6
• /
• pp.503-512
• /
• 2002

The drawbead is used to control material flow into the die during the binder wrap process and the stamping process in the sheet metal forming process. Since the dimension of drawbead is relatively small in comparison with the typical dimensions, it is difficult to include drawbeads in finite element analysis of the sheet metal forming process. It is because the mesh system has to be fine enough to describe the drawbead and the computation time is drastically increased. In this paper, simulation of drawbead forming has been carried out to obtain the equivalent boundary conditions in the binder wrap process and the stamping process. In order to investigate the effect of various die geometries, parameter studies are performed with the variation of parameters such as the blank length, the drawbead depth, the drawbead radius, the inclination of die and the friction coefficient.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2004.05a
• /
• pp.88-92
• /
• 2004

Mufflers have been important elements in engineering practice, which reduce various kinds of noises because of its general capability of application. Many kinds of methods :Ire applied to analyze their characteristics and to expect their performances. Some of the methods, conventionally, are based on the plane wave assumption for its simplicity. The shortcoming of this approach is its limitation of analysis ranges and parameters related with analysis and response. This research employs the boundary element method f3r the analysis of mufflers, which considers 3 dimensional scattering effects. This method can be used more appropriately fur tile analysis of mufflers because its analysis ability for higher frequency range which can be decided by the element size of the model. And, experimental analysis using a reciprocity theorem is conducted to verify the analysis results.

• Structural Engineering and Mechanics
• /
• v.54 no.3
• /
• pp.393-417
• /
• 2015

In this paper an 8-node quadrilateral assumed stress hybrid Mindlin plate element with $39{\beta}$ is presented. The formulation is based on complementary energy principle. The proposed element is free of shear locking and is capable of passing all the patch tests, especially the non-zero constant shear enhanced patch test. To accomplish this purpose, special attention is devoted to selecting boundary displacement interpolation and stress approximation in domain. The arbitrary order Timoshenko beam function is successfully used to derive the boundary displacement interpolation. According to the equilibrium equations, an appropriate stress approximation is rationally derived. Particularly, in order to improve element's accuracy, the assumed stress field is derived by employing $39{\beta}$ rather than conventional $21{\beta}$. The resulting element can be adopted to analyze both moderately thick and thin plates, and the convergence for the very thin case can be ensured theoretically. Excellent element performance is demonstrated by a wide of experimental evaluations.

• Structural Engineering and Mechanics
• /
• v.42 no.4
• /
• pp.471-488
• /
• 2012

An efficient one dimensional finite element model has been presented for the dynamic analysis of composite laminated beams, using the efficient layerwise zigzag theory. To meet the convergence requirements for the weak integral formulation, cubic Hermite interpolation is used for the transverse displacement ($w_0$), and linear interpolation is used for the axial displacement ($u_0$) and shear rotation (${\psi}_0$). Each node of an element has four degrees of freedom. The expressions of variationally consistent inertia, stiffness matrices and the load vector are derived in closed form using exact integration. The formulation is validated by comparing the results with the 2D-FE results for composite symmetric and sandwich beams with various end conditions. The employed finite element model is free of shear locking. The present zigzag finite element results for natural frequencies, mode shapes of cantilever and clamped-clamped beams are obtained with a one-dimensional finite element codes developed in MATLAB. These 1D-FE results for cantilever and clamped beams are compared with the 2D-FE results obtained using ABAQUS to show the accuracy of the developed MATLAB code, for zigzag theory for these boundary conditions. This comparison establishes the accuracy of zigzag finite element analysis for dynamic response under given boundary conditions.

• The Proceeding of the Korean Institute of Electromagnetic Engineering and Science
• /
• v.2 no.3
• /
• pp.26-31
• /
• 1991

The scattering property of TMz illuminated a elliptic dielectric cylinders with arbitrary cross section are analyzed by the boundary element techniques. The boundary element equations are for- mulated via Maxwell's equations, weighted residual of Green's theorem, and the boundary conditions. The unknown surface fields on the boundaries are then calculated by the boundary element integral equations. Once the surface fields are found, the scattered fields in far-zone and scattering widths (SW) are readily determined. To show the validity and usefulness of this formulation, computations are compared with those obtained using analytical method and one layer circular cylinder. As exten- sion to arbitrary cross-sectioned cylinders, plane wave scattering from a elliptic dielectric cylinders are numerically analyzed. A general computer program has been developed using the quadratic ele- ments(Higher order borndary elements) and the Gaussian quadrature.