• Journal of applied mathematics & informatics
• /
• v.17 no.1_2_3
• /
• pp.73-91
• /
• 2005

The application of Green's function in calculation of flow characteristics around submerged and floating bodies due to a regular wave is presented. It is assumed that the fluid is homogeneous, inviscid and incompressible, the flow is irrotational and all body motions are small. Two methods based on the boundary integral equation method (BIEM) are applied to solve associated problems. The first is a low order panel method with triangular flat patches and uniform distribution of velocity potential on each panel. The second method is a high order panel method in which the kernels of the integral equations are modified to make it nonsingular and amenable to solution by the Gaussian quadrature formula. The calculations are performed on a submerged sphere and some floating spheroids of different aspect ratios. The excellent level of agreement with the analytical solutions shows that the second method is more accurate and reliable.

• Wind and Structures
• /
• v.28 no.6
• /
• pp.347-354
• /
• 2019

This study provides an approximate analytical solution to the fractional vibration problem of thin plate governing anomalous motion of plate subjected to in-plane loads. The method of variable separable is employed to transform the fractional partial differential equations under consideration into a fractional ordinary differential equation in temporal variable and a bi-harmonic plate equation in spatial variable. The technique of conformable fractional derivative is utilized to solve the resulting fractional differential equation and the approach of finite sine integral transform method is used to solve the accompanying bi-harmonic plate equation. The deflection field which measures the transverse displacement of the plate is expressed in terms of product of Bessel and trigonometric functions via the temporal and spatial variables respectively. The obtained solution reduces to the solution of the free vibration problem of thin plate in literature. This work shows that conformable fractional derivative is an efficient mathematical tool for tracking analytical solution of fractional partial differential equation governing anomalous vibration of thin plates.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.7
• /
• pp.759-768
• /
• 2012

This paper proposes a method to determine the elastic follow-up factors for the $C(t)$-integral under secondary stress. The rate of creep crack growth for transient creep is correlated with the $C(t)$-integral. Elastic follow-up behavior, which occurs in structures under secondary loading, prevents a relaxation of stress during transient creep. Thus, both the values of $C(t)$ and creep crack growth increase as increasing elastic follow-up. An estimation solution for $C(t)$ was proposed by Ainsworth and Dean based on the reference stress method. To predict the value of $C(t)$ using this solution, an independent method to determine the elastic follow-up factors for cracked bodies is needed. This paper proposed that the elastic follow-up factors for $C(t)$ can be determined by elastic-plastic analyses using the plastic-creep analogy. Finite element analyses were performed to verify this method.

• Structural Engineering and Mechanics
• /
• v.78 no.5
• /
• pp.585-597
• /
• 2021

This paper presents a comparative study of analytical method, finite element method (FEM) and Multilayer Perceptron (MLP) for analysis of a contact problem. The problem consists of a functionally graded (FG) layer resting on a half plane and pressed with distributed load from the top. Firstly, analytical solution of the problem is obtained by using theory of elasticity and integral transform techniques. The problem is reduced a system of integral equation in which the contact pressure are unknown functions. The numerical solution of the integral equation was carried out with Gauss-Jacobi integration formulation. Secondly, finite element model of the problem is constituted using ANSYS software and the two-dimensional analysis of the problem is carried out. The results show that contact areas and the contact stresses obtained from FEM provide boundary conditions of the problem as well as analytical results. Thirdly, the contact problem has been extended based on the MLP. The MLP with three-layer was used to calculate the contact distances. Material properties and loading states were created by giving examples of different values were used at the training and test stages of MLP. Program code was rewritten in C++. As a result, average deviation values such as 0.375 and 1.465 was obtained for FEM and MLP respectively. The contact areas and contact stresses obtained from FEM and MLP are very close to results obtained from analytical method. Finally, this study provides evidence that there is a good agreement between three methods and the stiffness parameters has an important effect on the contact stresses and contact areas.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.28 no.9
• /
• pp.1284-1296
• /
• 2004

In order to investigate effects of anisotropic fiber packing on stresses in composites, a Volume Integral Equation Method is applied to calculate the elastostatic field in an unbounded isotropic elastic medium containing multiple orthotropic inclusions subject to remote loading, and a Mixed Volume and Boundary Integral Equation Method is introduced for the solution of elastostatic problems in unbounded isotropic materials containing multiple anisotropic inclusions as well as one void under uniform remote loading. A detailed analysis of stress fields at the interface between the isotropic matrix and the central orthotropic inclusion is carried out for square, hexagonal and random packing of orthotropic cylindrical inclusions, respectively. Also, an analysis of stress fields at the interface between the isotropic matrix and the central orthotropic inclusion is carried out, when it is assumed that a void is replaced with one inclusion adjacent to the central inclusion of square, hexagonal and random packing of orthotropic cylindrical inclusions, respectively, due to manufacturing and/or service induced defects. The effects of random orthotropic fiber packing on stresses at the interface between the isotropic matrix and the central orthotropic inclusion are compared with the influences of square and hexagonal orthotropic fiber packing on stresses. Through the analysis of plane elastostatic problems in unbounded isotropic matrix with multiple orthotropic inclusions and one void, it will be established that these new methods are very accurate and effective for investigating effects of general anisotropic fiber packing on stresses in composites.

• Journal of Mechanical Science and Technology
• /
• v.14 no.2
• /
• pp.197-206
• /
• 2000

A boundary integral equation method in the shape design sensitivity analysis is developed for the elasticity problems with axisymmetric non-homogeneous bodies. Functionals involving displacements and tractions at the zonal interface are considered. Sensitivity formula in terms of the interface shape variation is then derived by taking derivative of the boundary integral identity. Adjoint problem is defined such that displacement and traction discontinuity is imposed at the interface. Analytic example for a compound cylinder is taken to show the validity of the derived sensitivity formula. In the numerical implementation, solutions at the interface for the primal and adjoint system are used for the sensitivity. While the BEM is a natural tool for the solution, more generalization should be made since it should handle the jump conditions at the interface. Accuracy of the sensitivity is evaluated numerically by the same compound cylinder problem. The endosseous implant-bone interface problem is considered next as a practical application, in which the stress value is of great importance for successful osseointegration at the interface. As a preliminary step, a simple model with tapered cylinder is considered in this paper. Numerical accuracy is shown to be excellent which promises that the method can be used as an efficient and reliable tool in the optimization procedure for the implant design. Though only the axisymmetric problem is considered here, the method can be applied to general elasticity problems having interface.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.3
• /
• pp.261-275
• /
• 2016

In this paper, an integro-differential equation which arises in oscillating magnetic fields is studied. The generalized fractional order Chebyshev orthogonal functions (GFCF) collocation method used for solving this integral equation. The GFCF collocation method can be used in applied physics, applied mathematics, and engineering applications. The results of applying this procedure to the integro-differential equation with time-periodic coefficients show the high accuracy, simplicity, and efficiency of this method. The present method is converging and the error decreases with increasing collocation points.

• Journal of the Korean Operations Research and Management Science Society
• /
• v.25 no.3
• /
• pp.1-15
• /
• 2000

In this paper, we propose a practical cut generation method based on the Chvatal-Gomory procedure for the (0, 1)-Knapsack problem with a variable capacity. For a given set N of n items each of which has a positive integral weight and a facility of positive integral capacity, a feasible solution of the problem is defined as a subset S of N along with the number of facilities that can satisfy the sum of weights of all the items in S. We first derive a class of valid inequalities for the problem using Chvatal-Gomory procedure, then analyze the associated separation problem. Based on the results, we develop an affective cut generation method. We then analyze the theoretical strength of the inequalities which can be generated by the proposed cut generation method. Preliminary computational results are also presented which show the effectiveness of the proposed cut generation method.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.24 no.3
• /
• pp.1-8
• /
• 1987

In this paper, the localized finite element method(LFEM) is applied to 3-dimensional ship motion problems in water of infinite depth. The LFEM used here is based on the functional constructed by Bai & Yeung(1974). To test the present numerical scheme, a few vertical axisymmetric bodies are treated by general 3-dimensional formulation. The computed results of hydrodynamic coefficients for a few vertical spheroids and vertical circular cylinders show good agreement with results obtained by others. The advantages of the present numerical method compared with the method of integral equation are as follows; (i) The cumbersome existence of irregular frequencies in the method of conventional integral equation is removed. (ii) The final matrix is banded and symmetric and the computation of the matrix elements is comparatively easier, whereas the size of the matrix in the present scheme is much larger. (iii) In the future research, it is possible to accommodate with the nonlinear exact free surface boundary condition in the localized finite element subdomain, whereas the linear solution is assumed in the truncated(far field) subdomain.

• The Transactions of the Korean Institute of Electrical Engineers C
• /
• v.53 no.12
• /
• pp.626-633
• /
• 2004

In this paper, we present a time domain combined field integral equation (TD-CFIE) formulation to analyze the transient electromagnetic response from three-dimensional dielectric objects. The solution method in this paper is based on the method of moments (MoM) that involves separate spatial and temporal testing procedures. A set of the RWG (Rao, Wilton, Glisson) functions Is used for spatial expansion of the equivalent electric and magnetic current densities and a combination of RWG and its orthogonal component is used as spatial testing. We also investigate spatial testing procedures for the TD-CFIE to select the proper testing functions that are derived from the Laguerre polynomials. These basis functions are also used for temporal testing. Use of this temporal expansion function characterizing the time variable enables one to handle the time derivative terms in the integral equation and decouples the space-time continuum in an analytic fashion. Numerical results computed by the proposed formulation are presented and compared with the solutions of the frequency domain combined field integral equation (FD-CFIE).