• Interaction and multiscale mechanics
• /
• v.2 no.2
• /
• pp.177-187
• /
• 2009

Machine or structural members subjected to fatigue loading will have a crack initiated during early part of their life. Therefore analysis of members with cracks and other discontinuities is very important. Finite element method has enjoyed widespread use in engineering, but it is not convenient for crack problems as the region very close to crack tip is to be discretized with very fine mesh. However, as the body force method (BFM), requires only the boundary of the discontinuity (crack or hole) to be discretized it is easy versatile technique to analyze such problems. In the present work fundamental solution for concentrated load x + iy acting in the semi-infinite plate at an arbitrary point $z_0=x_0+iy_0$ is considered. These fundamental solutions are in complex form ${\phi}(z)$ and ${\psi}(z)$ (England 1971). These potentials are known as Melan potentials (Ramakrishna 1994). A crack in the semi-infinite plate as shown in Fig. 1 is considered. This crack is divided into number of divisions. By applying pair of body forces on a division, the resultant forces on the remaining 'N'divisions are to be found for which ${\phi}_1(z)$ and ${\psi}_1(z)$ are derived. Body force method is applied to calculate stress intensity factor for crack in semi-infinite plate. Also for the case of crack emanating from circular hole in semi-infinite plate radial stress, hoop stress and shear stress are calculated around the hole and crack. Convergent results are obtained by body force method. These results are compared with FEM results.

• Structural Engineering and Mechanics
• /
• v.19 no.3
• /
• pp.281-295
• /
• 2005

An analytical method is presented to solve the elastodynamic problem of finitely long hollow cylinder subjected to torsional impact often occurs in engineering mechanics. The analytical solution is composed of a solution of quasi-static equation satisfied with the non-homogeneous boundary condition and a solution of dynamic equation satisfied with homogeneous boundary condition. The quasi-static solution is obtained directly by solving the quasi-static equation satisfied with the non-homogeneous boundary condition. The solution of the non-homogeneous dynamic equation is obtained by means of finite Hankel transform on the radial variable, r, Laplace transform on time variable, t, and finite Fourier transform on axial variable, z. Thus, the solution for finitely long, hollow cylinder subjected to torsion impact is obtained. In the calculating examples, the response histories and distributions of shear stress in the finitely long hollow cylinder subjected to an exponential decay torsion load are obtained, and the results have been analyzed and discussed. Finally, a dynamic finite element for the same problem is carried out by using ABAQUS finite element analysis. Comparing the analytical solution with the finite element solution, it can be found that two kinds of results obtained by means of two different methods agree well. Therefore, it is further concluded that the analytical method and computing process presented in the paper are effective and accurate.

• Journal of the Korean Geotechnical Society
• /
• v.24 no.6
• /
• pp.77-84
• /
• 2008

The existing equations for radial consolidation cannot account for the changes of well resistance with time and cannot predict the appropriate in-situ consolidation curve. In this study, small cylinder cell tests are performed to evaluate the discharge capacity of PVD. Also, a block sample of 1.2 m in diameter and 2.0 m in height was consolidated to observe the change in the drainage capacity with time for three types of PVD. From the test results on a block sample, the drainage curves normalized with initial drainage of each PVD are similar, regardless of the PVD type and the consolidation curve, which is predicted using solutions of radial consolidation based on the discharge capacity measured in a small cylinder cell tests, significantly overestimates the degree of consolidation. The term of well resistance in the radial consolidation solution was back-calculated to fit the consolidation curve of a large block sample and it is defined as the time dependent well resistance factor, L(t). The L(t) was found to be linearly proportional to the dimensionless time factor, Th. It was also shown that the consolidation curve evaluated by using L(t) provides more accurate prediction than the existing solution.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.4
• /
• pp.981-989
• /
• 1995

A new weight function approach to determine SIF(stress intensity factor) using single-layer potential has been presented. The crack surface displacement field was represented by one boundary integral term whose kernel was modified from Kelvin's fundamental solution. The proposed method enables the calculation of SIF using only one SIF solution without any modification for the crack geometries symmetric in two-dimensional plane such as a center crack in a plate with or without an internal hole, double edge cracks, circumferential crack or radial cracks in a pipe. The application procedure to those crack problems is very simple and straightforward with only one SIF solution. The necessary information in the analysis is two reference SIFs. The analysis results using present closed-form solution were in good agreement with those of the literature.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.56 no.3
• /
• pp.451-456
• /
• 2007

This paper presents a problem that reconfigure distribution power system using branch exchange method. Optimal reconfiguration problem calculates line loss, voltage condition about system states of all situations that become different according to line On/off status, and search for optimum composition of these. However, result is difficult to be calculated fast. Because radiated operation condition of system is satisfied using many connection and sectionalize switches in the distribution power system. Therefore, in this paper, optimization method for reducing system total loss and satisfying operating condition of radial and constraints condition of voltage is proposed using the fastest branch exchange. And optimal solution at branch exchange algorithm can be wrong estimated to local optimal solution according to initial operating state. Considering this particular, an initial operating point algorithm is added and this paper showed that optimal solution arrives at global optimal solution.

• Geomechanics and Engineering
• /
• v.29 no.1
• /
• pp.1-11
• /
• 2022

A semi-analytical solution to responses of overconsolidated (OC) unsaturated soils surrounding an expanding spherical cavity under constant suction condition is presented. To capture the elastoplastic hydro-mechanical property of OC unsaturated soils, the unified hardening (UH) model for OC unsaturated soil is adopted in corporation with a soil-water characteristic curve (SWCC) and two suction yield surfaces. Taking the specific volume, radial stress, tangential stress and degree of saturation as the four basic unknowns, the problem investigated is formulated by solving a set of first-order ordinary differential equations with the help of an auxiliary variable and an iterative algorithm. The present solution is validated by comparing with available solution based on the modified Cam Clay (MCC) model. Parametric studies reveal that the hydraulic and mechanical responses of spherical cavity expanding in unsaturated soils are not only coupled, but also affected by suction and overconsolidation ratio (OCR) significantly. More importantly, whether hydraulic yield will occur or not depends only on the initial relationship between suction yield stress and suction. The presented solution can be used for calibration of some insitu tests in OC unsaturated soil.

• Structural Engineering and Mechanics
• /
• v.67 no.6
• /
• pp.659-669
• /
• 2018

The Dynamic equations in the polar coordinates are drawn out using the MLPG method for the non-symmetric FG cylindrical shell. To simulate the mechanical properties of FGM, the nonlinear volume fractions for radial direction are used. The shape function applied in this paper is a form of the radial basis functions, by using this function all the requirements for an effective and suitable shape function are established. Hence in this study, the multiquadrics (MQ) radial basis functions are exploited as the shape function governing the problem. The MLPG method is combined with the the Newmark time approximation scheme to solve dynamic equations in the time domain. The obtained results by the MLPG method to be verified are compared with the analytical solution and the FEM. The obtained results through the MLPG method show a good agreement in comparison to other results and the MLPG method has high accuracy for dynamic analysis of the non-symmetric FG cylindrical shell. To demonstrate the capability of the present method to dynamic analysis of the non-symmetric FG cylindrical shell, it is analyzed dynamically with different volume fraction exponents under harmonic and rectangular shock loading. The present method shows high accuracy, efficiency and capability to dynamic analysis of the non-symmetric FG cylindrical shell with nonlinear grading patterns.

• Bulletin of the Korean Mathematical Society
• /
• v.50 no.3
• /
• pp.839-865
• /
• 2013

We deal with the existence of positive radial solutions concentrating on spheres for the following class of elliptic system $$\large(S) \hfill{400} \{\array{-{\varepsilon}^2{\Delta}u+V_1(x)u=K(x)Q_u(u,v)\;in\;\mathbb{R}^N,\\-{\varepsilon}^2{\Delta}v+V_2(x)v=K(x)Q_v(u,v)\;in\;\mathbb{R}^N,\\u,v{\in}W^{1,2}(\mathbb{R}^N),\;u,v&gt;0\;in\;\mathbb{R}^N,}$$ where ${\varepsilon}$ is a small positive parameter; $V_1$, $V_2{\in}C^0(\mathbb{R}^N,[0,{\infty}))$ and $K{\in}C^0(\mathbb{R}^N,[0,{\infty}))$ are radially symmetric potentials; Q is a $(p+1)$-homogeneous function and p is subcritical, that is, 1 < $p$ < $2^*-1$, where $2^*=2N/(N-2)$ is the critical Sobolev exponent for $N{\geq}3$.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.2
• /
• pp.243-253
• /
• 1999

The dual reciprocity boundary element method (DRBEM) is used to solve the Graetz problem of laminar flow inside circular duct. In this method the domain integral tenn of boundary integral equation resulting from source term of governing equation is transformed into equivalent boundary-only integrals by using the radial basis interpolation function, and therefore complicate domain discretization procedure Is completely removed. Velocity profile is obtained by solving the momentum equation first and then, using this velocities as Input data, energy equation Is solved to get the temperature profile by advancing from duct entrance through the axial direction marching scheme. DRBEM solution is tested for the uniform temperature and heat flux boundary condition cases. Local Nusselt number, mixed mean temperature and temperature profile inside duct at each dimensionless axial location are obtained and compared with exact solutions for the accuracy test Solutions arc in good agreement at the entry region as well as fully developed region of circular duct, and their accuracy are verified from error analysis.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.19 no.11
• /
• pp.1133-1143
• /
• 2009

Vibro-acoustic characteristics of a simplified disk-brake rotor containing a narrow radial slot are studied using a semi-analytical procedure. First, modal sound radiations for flexural and radial modes of a generic annular disk having identical key dimension and slot(with free boundaries) are defined using pre-developed analytical solutions based on the modal vibrations from finite element model. The analytical solutions are validated by fully computational methods. Second, sound radiation from a simplified brake rotor simulated using sound radiation solution of the generic disk based on the rotor eigensolutions computed using a finite element code. Predictions by the semi-analytical method matched well numerical calculations using finite element and boundary element method. Finally, sound radiation and vibration characteristics for the example rotor due to a harmonic excitation fixed to the rotor or rotating around the rotor are also obtained.