• Journal of Mechanical Science and Technology
• /
• v.19 no.8
• /
• pp.1611-1620
• /
• 2005

This paper shows the use of wavelet transformation combined with inverse acoustics to reconstruct the surface velocity of a noise source. This approach uses the boundary element analysis based on the measured sound pressure at a set of field points, the Helmholtz integral equations and wavelet transformation for reconstructing the normal surface velocity field. The reconstructed field can be diverged due to the small measurement errors in the case of nearfield acoustic holography (NAH) using an inverse boundary element method. In order to avoid this instability in the inverse problem, the reconstruction process should include some form of regularization for enhancing the resolution of source images. The usual method of regularization has been the truncation of wave vectors associated with small singular values, although the order of an optimal truncation is difficult to determine. In this paper, a wavelet transformation is applied to reduce the computation time for inverse acoustics and to enhance the reconstructed vibration field. The computational speed-up is achieved, with solution time being reduced to $14.5\%$.

• Journal of Mechanical Science and Technology
• /
• v.17 no.8
• /
• pp.1120-1132
• /
• 2003

A problem of a circular elastic inhomogeneity interacting with a crack under uniform loadings (mechanical tension and heat flux at infinity) is solved. The singular. integral equations for edge and temperature dislocation distribution functions are constructed and solved numeric-ally, to obtain the stress intensity factors. The effects of the material property ratio on the stress intensity factor (SIF) are investigated. The computed SIFs are used to predict the kink angle of the crack when the crack grows.

• Advances in materials Research
• /
• v.3 no.1
• /
• pp.237-257
• /
• 2014

In this paper, we adopted a two-dimensional analytical electro-elastic model to predict the stress distributions of the piezoelectric actuator in 3D case. The actuator was embedded in an elastic host structure under electrical loadings. The problem is reduced to the solution of singular integral equations of the first kind. The interfacial stresses and the axial normal stress in both plane stress state and plane strain state were obtained to study the actuation effects being transferred from the actuator to the host. The stress distributions of the PZT actuator in different length and different thickness were analyzed to guarantee the generality. The validity of the present model has been demonstrated by application of specific examples and comparisons with the corresponding results obtained from the Finite Element Method.

• Structural Engineering and Mechanics
• /
• v.73 no.5
• /
• pp.585-598
• /
• 2020

We develop design model within which nucleation and propagation of crack in a fibrous composite is described. It is assumed that under loading, crack initiation and fracture of material happens in the composite. The problem of equilibrium of a composite with embryonic crack is reduced to the solution of the system of nonlinear singular integral equations with the Cauchy type kernel. Normal and tangential forces in the crack nucleation zone are determined from the solution of this system of equations. The crack appearance conditions in the composite are formed with regard to criterion of ultimate stretching of the material's bonds. We study the case when near the fiber, the binder has several arbitrary arranged rectilinear prefracture zones and a crack with interfacial bonds. The proposed computational model allows one to obtain the size and location of the zones of damages (prefracture zones) depending on geometric and mechanical characteristics of the fibrous composite and applied external load. Based on the suggested design model that takes into account the existence of damages (the zones of weakened interparticle bonds of the material) and cracks with end zones in the composite, we worked out a method for calculating the parameters of the composite, at which crack nucleation and crack growth occurs.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.22 no.4
• /
• pp.779-789
• /
• 1998

In this paper, the thermal shock responses of collinear cracks in a layered medium are investigated based on the uncoupled, quasi-static plane thermoelasticity. The medium is modeled as a bonded structure composed of a surface layer and a semi-infinite substrate. Between these two dissimilar homogeneous constituents, a functionally graded interfacial zone exists with the nonhomogeneous features of continuously varying thermoelastic properties. Three cracks are assumed to be present in the layered medium, one in each one of the constituent materials, aligned collinearly normal to the nominal interfaces. A system of singular integral equations is solved, subjected to the forcing terms of equivalent transient thermal tractions acting on the locations of cracks via superposition. Main results presented are the transient thermal stress intensity factors to illustrate the parametric effects of various geometric and amterial combinations of the medium with the thermoelastically graded interfacial zone and the collinear cracks.

• Structural Engineering and Mechanics
• /
• v.86 no.4
• /
• pp.535-546
• /
• 2023

In the present paper, multiple interface cracks between a functionally graded orthotropic coating and an orthotropic half-plane substrate under concentrated loading are considered by means of the distribution dislocation technique (DDT). With the use of integration of Fourier transform the problem is reduced to a system of Cauchy-type singular integral equations which are solved numerically to compute the dislocation density on the surfaces of the cracks. The distribution dislocation is a powerful method to calculate accurate solutions to plane crack problems, especially this method is very good to find SIFs for multiple unequal cracks located at the interface. Hence this technique allows considering any number of interface cracks. The primary objective of this paper is to investigate the effects of the interaction of multiple interface cracks, load location, material orthotropy, nonhomogeneity parameters and geometry parameters on the modes I and II SIFs. Numerical results show that modes I/II SIFs decrease with increasing the nonhomogeneity parameter and the highest magnitude of SIF occurs where distances between the load location and crack tips are minimal.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.38-52
• /
• 1991

In this paper, a semi-Lagrangian method is used to solve the nonlinear hydrodynamics of a three-dimensional body beneath the free surface in the time domain. The boundary value problem is solved by using the boundary integral method. The geometries of the body and the free surface are represented by the curved panels. The surfaces are discretized into the small surface elements using a bi-cubic B-spline algorithm. The boundary values of $\phi$ and $\frac{\partial{\phi}}{\partial{n}}$ are assumed to be bilinear on the subdivided surface. The singular part proportional to $\frac{1}{R}$ are subtracted off and are integrated analytically in the calculation of the induced potential by singularities. The far field flow away from the body is represented by a dipole at the origin of the coordinate system. The Runge-Kutta 4-th order algorithm is employed in the time stepping scheme. The three-dimensional form of the integral equation and the boundary conditions for the time derivative of the potential Is derived. By using these formulas, the free surface shape and the equations of motion are calculated simultaneously. The free surface shape and fille forces acting on a body oscillating sinusoidally with large amplitude are calculated and compared with published results. Nonlinear effects on a body near the free surface are investigated.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.2
• /
• pp.207-216
• /
• 2003

In this paper, the finite element method for the elastodynamic axisymmetric fracture analysis is presented in matrix form through the application of the Galerkin method to the time integral equations of motion with no inertia forces. Isoparametric quadratic quadrilateral element and triangular crack tip singular elements with one-quarter node are used in the mesh division of the finite element model. To show the validity and accuracy of the proposed method, the infinite elastic medium with the penny shaped crack is solved first and compared with the analytical solution and the numerical results by the finite difference method and the boundary element method existing in the published literatures, and then the dynamic stress intensity factors of solid and hollow cylinders of finite dimensions haying penny-shaped cracks and internal and external circumferential tracks are computed in detail.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.10
• /
• pp.3048-3057
• /
• 1996

Analysis model containing two inclined surface cracks on semi-infinite elastic body is established and analyzed on the basis of linear fracture mechanics to examine mutual interference of two surface cracks. Muskhelishvili's complex stress functions are introduced and a set of singular integral equations is obtained for a dislocation density function. The stress intensity factors at crack tip are obtained by using the Gerasoulis'method. When two surface cracks are parallel and have the same length, the values of $K_1$and $\Delta K_11$(variation of $K_11$) for crack 1 and crack 2 decrease by the mutual interference of two surface cracks as the distance between the two surface cracks shortens. The effect of mutual interference is remarkable in high friction coefficient. In case that two surface cracks are parallel, the values of $K_1$and $\Delta K_11$for crack 2 decrease as the length ratio ot crack 2 to crack 1 becomes small. As the crack inclination angle rises, the value of $K_1$ and the mutual interference of $K_1$for crack 2 increase and the value of$\Delta K_11$ for crack 1 becomes smaller than that for crack 2.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.1
• /
• pp.117-124
• /
• 2009

The p-convergent boundary element method has been proposed to analyze two-dimensional potential problem on the basis of high order Legendre shape functions that have different property comparing with the shape functions in conventional boundary element method. The location of nodes corresponding to high order shape function are not defined along the boundary, called by nodeless node, similar to the p-convergent finite element method. As the order of shape function increases, the collocation point method is used to solve linear simultaneous equations. The collocation patterns of p-convergent boundary element method consist of non-symmetric hierarchial or symmetric non-hierarchical. As the order of shape function increases, the number of collocation point increases. The singular integral that appears in p-convergent boundary element has been calculated by special numeric quadrature technique and semi-analytical integration technique. The L-shape domain problem including singularity in the vicinity of reentrant comer is analyzed and the numerical results show that the relative error is smaller than $10^{-2}%$ range as compared with other results in literatures. In case of same condition, the symmetric p-collocation point pattern shows high accuracy of solution.