• 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.

• 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.

• Journal of the Korean Institute of Landscape Architecture
• /
• v.33 no.4 s.111
• /
• pp.97-107
• /
• 2005

This study was carried out to construct accurate and scientific system of assessment indicators in selection of National Trust conservation areas, which was new concept of domestic environment movement and offer the raw data of new analytic method by introducing the fuzzy theory and weight for overcoming the uncertainty of ranking decision. To transform the Likert's scale granted to assessment indicators into the type of triangular fuzzy number(a, b, c), there was conversion to each minimum(a), median(b), and maximum(c) in applying membership function, and in using the center of gravity and eigenvalue, there was to decide the ranking. The rankings of converted values applied a mean importance and weight were confirmed that they were generally changed. Therefore, the ranking decision was better to accomplish objective and rational ranking decision by applying weight that was calculated in grouping of indicator than to judge the singular concept and to be useful in assessment of diverse National Trust site. In the future, because AHP, which was general method of calculating weight, was lacked, there was to understand the critical point to fix a pertinent weight, and to carry out the study applying engineering concept like fuzzy integral using $\lambda-measure$.

• 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.