• Journal of the Korean Society for Precision Engineering
• /
• v.16 no.4 s.97
• /
• pp.197-204
• /
• 1999

The boundary element method was used for studying singularities of an interface crack with contact zones. The iterative procedure is applied to estimate the contact zone size. Because the contact zone size was extremely small in a tension field, a large number of Gaussian points were used for numerical integration of the Kernels. Stress extrapolation method and J-integral were used ofr determining stress intensity factors. When the interface crack was assumed to have opened tips, oscillatory singularities appear near the tips of the interface crack. But the interface crack with contact zone which Comninou suggested had no oscillatory behavior. The contact zone size under shear loading was much larger than that under tensile. The stress intensity factors computed by stress extrapolation method were close to those of Comninou's solution. And the stress intensity factor evaluated by J-integral was similar to that by stress extrapolation method.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.6
• /
• pp.1445-1451
• /
• 1994

To calculate the wave pressure acting on coastal structures under the design wave condition, it is often necessary to numerically reproduce the big standing wave profiles close to wave breaking condition. For this, the governing equation and all nonlinear terms occurring in boundary conditions should be effectively considered in the numerical wave profile. In particular, the velocity square term in the free surface boundary condition is very important. A boundary element method is applied here to calculate the standing wave profile with the velocity square term fully treated by Newton iterative method. In order to check the validity of the method, the numerical wave profiles are compared to ones calculated by the perturbation method, the Fourier approximation method and the hydraulic experiment.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2006.04a
• /
• pp.291-298
• /
• 2006

A parallelized topology design optimization method is developed on a distributed memory system. The parallelization is based on a domain decomposition method and a boundary communication scheme. For the finite element analysis of structural responses and design sensitivities, the PCG method based on a Krylov iterative scheme is employed. Also a parallelized optimization method of optimality criteria is used to solve large-scale topology optimization problems. Through several numerical examples, the developed method shows efficient and acceptable topology optimization results for the large-scale problems.

• Ocean Systems Engineering
• /
• v.6 no.3
• /
• pp.245-264
• /
• 2016

The iterative boundary element method (IBEM) developed originally before for cavitating two-dimensional (2-D) and three-dimensional (3-D) hydrofoils moving under free surface is modified and applied to the case of 2-D (two-dimensional) airfoils and 3-D (three-dimensional) wings over water. The calculation of the steady-state flow characteristics of an inviscid, incompressible fluid past 2-D airfoils and 3-D wings above free water surface is of practical importance for air-assisted marine vehicles such as some racing boats including catamarans with hydrofoils and WIG (Wing-In-Ground) effect crafts. In the present paper, the effects of free surface both on 2-D airfoils and 3-D wings moving steadily over free water surface are investigated in detail. The iterative numerical method (IBEM) based on the Green's theorem allows separating the airfoil or wing problems and the free surface problem. Both the 2-D airfoil surface (or 3-D wing surface) and the free surface are modeled with constant strength dipole and constant strength source panels. While the kinematic boundary condition is applied on the airfoil surface or on the wing surface, the linearized kinematic-dynamic combined condition is applied on the free surface. The source strengths on the free surface are expressed in terms of perturbation potential by applying the linearized free surface conditions. No radiation condition is enforced for downstream boundary in 2-D airfoil and 3-D wing cases and transverse boundaries in only 3-D wing case. The method is first applied to 2-D NACA0004 airfoil with angle of attack of four degrees to validate the method. The effects of height of 2-D airfoil from free surface and Froude number on lift and drag coefficients are investigated. The method is also applied to NACA0015 airfoil for another validation with experiments in case of ground effect. The lift coefficient with different clearance values are compared with those of experiments. The numerical method is then applied to NACA0012 airfoil with the angle of attack of five degrees and the effects of Froude number and clearance on the lift and drag coefficients are discussed. The method is lastly applied to a rectangular 3-D wing and the effects of Froude number on wing performance have been investigated. The numerical results for wing moving under free surface have also been compared with those of the same wing moving above free surface. It has been found that the free surface can affect the wing performance significantly.

• Structural Engineering and Mechanics
• /
• v.52 no.1
• /
• pp.115-136
• /
• 2014

In this study, an extended Kantorovich method, employing multi-term displacement functions, is applied to analyze the vibration problem of symmetrically laminated plates with arbitrary boundary conditions. The vibration behaviors of laminated plates are determined based on the variational principle of total energy minimization and the iterative Kantorovich method. The out-of-plane displacement is represented in the form of a series of a sum of products of functions in x and y directions. With a known function in the x or y directions, the formulation for the variation of total potential energy is transformed to a set of governing equations and a set of boundary conditions. The equations and boundary conditions are then numerically solved for the natural frequency and vibration mode shape. The solutions are verified with available solutions from the literature and solutions from the Ritz and finite element analysis. In most cases, the natural frequencies compare very well with the reference solutions. The vibration mode shapes are also very well modeled using the multi-term assumed displacement function in the terms of a power series. With the method used in this study, it is possible to solve the angle-ply plate problem, where the Kantorovich method with single-term displacement function is ineffective.

• Structural Engineering and Mechanics
• /
• v.16 no.3
• /
• pp.259-274
• /
• 2003

In this paper, a novel numerical solution technique, the differential cubature method is employed to study the buckling problems of thick plates with arbitrary quadrilateral planforms and non-uniform boundary constraints based on the first order shear deformation theory. By using this method, the governing differential equations at each discrete point are transformed into sets of linear homogeneous algebraic equations. Boundary conditions are implemented through discrete grid points by constraining displacements, bending moments and rotations of the plate. Detailed formulation and implementation of this method are presented. The buckling parameters are calculated through solving a standard eigenvalue problem by subspace iterative method. Convergence and comparison studies are carried out to verify the reliability and accuracy of the numerical solutions. The applicability, efficiency, and simplicity of the present method are demonstrated through solving several sample plate buckling problems with various mixed boundary constraints. It is shown that the differential cubature method yields comparable numerical solutions with 2.77-times less degrees of freedom than the differential quadrature element method and 2-times less degrees of freedom than the energy method. Due to the lack of published solutions for buckling of thick rectangular plates with mixed edge conditions, the present solutions may serve as benchmark values for further studies in the future.

• Structural Engineering and Mechanics
• /
• v.88 no.5
• /
• pp.439-449
• /
• 2023

Non-linear vibration characteristics of functionally graded CNT-reinforced composite (FG-CNTRC) cylindrical shell panel on elastic foundation have not been sufficiently examined. In this situation, this study aims at the profound numerical investigation of the non-linear vibration response of FG-CNTRC cylindrical panels on Winkler-Pasternak foundation by introducing an accurate and effective 2-D meshfree-based non-linear numerical method. The large-amplitude free vibration problem is formulated according to the first-order shear deformation theory (FSDT) with the von Karman non-linearity, and it is approximated by Laplace interpolation functions in 2-D natural element method (NEM) and a non-linear partial derivative operator H_NL. The complex and painstaking numerical derivation on the curved surface and the crucial shear locking are overcome by adopting the geometry transformation and the MITC3+ shell elements. The derived nonlinear modal equations are iteratively solved by introducing a three-step iterative solving technique which is combined with Lanczos transformation and Jacobi iteration. The developed non-linear numerical method is estimated through the benchmark test, and the effects of foundation stiffness, CNT volume fraction and functionally graded pattern, panel dimensions and boundary condition on the non-linear vibration of FG-CNTRC cylindrical panels on elastic foundation are parametrically investigated.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2009.04a
• /
• pp.534-535
• /
• 2009

This paper deals with BEM analysis of transient elastodynamic problems using domain decomposition method and particular integrals. The particular method is used to approximate the acceleration term in the governing equation. The domain decomposition method is examined to consider multi-region problems. The domain of the original problem is subdivided into sub-regions, which are modeled by the particular integral BEM. The iterative coupling employing Schwarz algorithm is used for the successive update of the interface boundary conditions until convergence is achieved. The numerical results, compared with those by ABAQUS, demonstrate the validity of the present formulation.

• Journal of the Korean Mathematical Society
• /
• v.34 no.4
• /
• pp.841-857
• /
• 1997

We present analysis and some computational methods for boundary optimal and feedback control problems for Navier-Stokes equations. We use one example to illustrate our methodology and ideas which are applicable to general control problems for Navier-Stokes equations. First, we discuss the existence of optimal solutions and derive an optimality system of equations from which an optimal solution may be computed. Then we present a gradient type iterative method. Finally, we present some numerical results.

• Computers and Concrete
• /
• v.11 no.5
• /
• pp.441-462
• /
• 2013

When using the conventional finite element method, a great number of grid nodes are necessary to describe the large and uneven temperature gradients in the concrete around cooling pipes when calculating the temperature field of mass concrete with cooling pipes. In this paper, the temperature gradient properties of the concrete around a pipe were studied. A new calculation method was developed based on these properties and an explicit iterative algorithm. With a small number of grid nodes, both the temperature distribution along the cooling pipe and the temperature field of the concrete around the water pipe can be correctly calculated with this new method. In conventional computing models, the cooling pipes are regarded as the third boundary condition when solving a model of concrete with plastic pipes, which is an approximate way. At the same time, the corresponding parameters have to be got by expensive experiments and inversion. But in the proposed method, the boundary condition is described strictly, and thus is more reliable and economical. And numerical examples were used to illustrate that this method is accurate, efficient and applicable to the actual engineering.