• Structural Engineering and Mechanics
• /
• v.57 no.6
• /
• pp.1143-1156
• /
• 2016

This paper presents a comparative study of finite element method (FEM) and analytical method for the plane problem of a layered composite containing an internal perpendicular crack in literature. The layered composite consists of two elastic layers having different elastic constants and heights. External load is applied to the upper elastic layer by means o a rigid punch and the lower elastic layer rests on two simple supports. Numerical simulations are realized by the world wide code ANYS software. Two dimensional analysis of the problem is carried out and the results are verified by comparison with solutions reported in literature. Main goal of the numerical simulation is to investigate the normal stress ${\sigma}_x$(0, y), stress intensity factors at the crack factor and the crack opening displacements.

• Computers and Concrete
• /
• v.12 no.3
• /
• pp.261-283
• /
• 2013

This paper aims to contribute to the three-dimensional generalization of numerical prediction of crack propagation through the formulation of finite elements with embedded discontinuities. The analysis of crack propagation in two-dimensional problems yields lines of discontinuity that can be tracked in a relatively simple way through the sequential construction of straight line segments oriented according to the direction of failure within each finite element in the solid. In three-dimensional analysis, the construction of the discontinuity path is more complex because it requires the creation of plane surfaces within each element, which must be continuous between the elements. In the method proposed by Chaves (2003) the crack is determined by solving a problem analogous to the heat conduction problem, established from local failure orientations, based on the stress state of the mechanical problem. To minimize the computational effort, in this paper a new strategy is proposed whereby the analysis for tracking the discontinuity path is restricted to the domain formed by some elements near the crack surface that develops along the loading process. The proposed methodology is validated by performing three-dimensional analyses of basic problems of experimental fractures and comparing their results with those reported in the literature.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.7
• /
• pp.1255-1265
• /
• 1992

Under uniform heat flow, the thermal stress intensity factors for the interfacial crack on a rigid cusp-type inclusion are determined by Hilbert problem expressed with complex variable. The thermal stress intensity factors are expressed in terms of the periodic function of heat flow angle. When the tip of the interfacial crack meets that of the cusp crack, the thermal stress intensity factors have singularities. The thermal stress intensity factors at the interfacial crack tip located in the distance from the cusp crack tip vary with the location of the interfacial crack tip. From the results of the analysis, the complex potential functions and the thermal stress intensity factors for the cusp-type inclusion without the interfacial crack are derived under the cusp surface boundary conditions insulated or fixed to zero relative temperature.

• Structural Engineering and Mechanics
• /
• v.25 no.5
• /
• pp.501-518
• /
• 2007

In this work, a new approach is developed for dynamic analysis of a composite beam with an interply crack, based on finite element solution of partial differential equations with the use of the COMSOL Multiphysics package, allowing for fast and simple change of geometric characteristics of the delaminated area. The use of COMSOL Multiphysics package facilitates automatic mesh generation, which is needed if the problem has to be solved many times with different crack lengths. In the model, a physically impossible interpenetration of the crack faces is prevented by imposing a special constraint, leading to taking account of a force of contact interaction of the crack faces and to nonlinearity of the formulated boundary value problem. The model is based on the first order shear deformation theory, i.e., the longitudinal displacement is assumed to vary linearly through the beam's thickness. The shear deformation and rotary inertia terms are included into the formulation, to achieve better accuracy. Nonlinear partial differential equations of motion with boundary conditions are developed and written in the format acceptable by the COMSOL Multiphysics package. An example problem of a clamped-free beam with a piezoelectric actuator is considered, and its finite element solution is obtained. A noticeable difference of forced vibrations of the delaminated and undelaminated beams due to the contact interaction of the crack's faces is predicted by the developed model.

• Structural Engineering and Mechanics
• /
• v.23 no.1
• /
• pp.63-74
• /
• 2006

In this paper, the scattering of harmonic elastic anti-plane shear waves by two collinear cracks in functionally graded materials is investigated by means of nonlocal theory. The traditional concepts of the non-local theory are extended to solve the fracture problem of functionally graded materials. To overcome the mathematical difficulties, a one-dimensional non-local kernel is used instead of a two-dimensional one for the anti-plane dynamic problem to obtain the stress field near the crack tips. To make the analysis tractable, it is assumed that the shear modulus and the material density vary exponentially with coordinate vertical to the crack. By use of the Fourier transform, the problem can be solved with the help of a pair of triple integral equations, in which the unknown variable is the displacement on the crack surfaces. To solve the triple integral equations, the displacement on the crack surfaces is expanded in a series of Jacobi polynomials. Unlike the classical elasticity solutions, it is found that no stress singularities are present at crack tips.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.11 s.254
• /
• pp.1384-1391
• /
• 2006

The problem of an infinite anisotropic material with a crack inclined with respect to the principal material axes is analyzed. The material is subjected to uniform biaxial load along its boundary. It is assumed that the material is homogeneous, but anisotropic. By considering the effect of the horizontal load, the distribution of stresses at the crack tip is analyzed. The problem of predicting critical stress in anisotropic solids which is a subject of considerable practical importance is examined and the effect of load biaxiality is made explicitly. The present results based on the normal stress ratio theory show significant effects of biaxial load, crack inclination angle and fiber orientation on the critical stress. The analysis is performed for a wide range of the crack angles and biaxial loads.

• Journal of Mechanical Science and Technology
• /
• v.17 no.2
• /
• pp.269-279
• /
• 2003

The anti-plane shear problem of bonded elastic materials containing a crack at an arbitrary angle to the graded interfacial zone is investigated in this paper The interfacial zone is modeled as a nonhomogeneous interlayer of finite thickness with the continuously varying shear modulus between the two dissimilar, homogeneous half-planes. Formulation of the crack problem is based upon the use of the Fourier integral transform method and the coordinate transformations of basic field variables. The resulting Cauchy-type singular integral equation is solved numerically to provide the values of mode 111 stress intensity factors. A comprehensive parametric study is then presented of the influence of crack obliquity on the stress intensity factors for different crack size and locations and for different material combinations, in conjunction with the material nonhomogeneity within the graded interfacial zone.

• Proceedings of the KWS Conference
• /
• 2004.05a
• /
• pp.264-266
• /
• 2004

for the kぉ mechanics analysis of a crack in a weld of dissimilar steels, residual stress analysis and fracture analysis must be performed simultaneously. The standard definition of the J-integral leads to a path dependent value in the presence of a residual stress field. And unlike cracks in homogeneous materials, a bimaterial interface crack always induces both opening and shearing modes of stress in the vicinity of the crack tip. Therefore, it is necessary to develope a path independent J-integral definition for a crack in a residual stress field generated by welding of dissimilar steels. This paper addresses the modification of the Rice-J-integral to produce a path independent J-integral when residual stresses due to welding of dissimilar steels and external forces are present. The residual stress problem is heated as an initial stain problem and the J-integral proposed for this class of problems is used And a program which can evaluate the 1-integral for a crack in a weld of dissimialr steels is developed using proposed J-integral definition.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.7
• /
• pp.1590-1596
• /
• 1995

The energy release rates for a dynamically growing crack in orthotropic materials are expressed explicitly in terms of dynamic stress intensity factors. The stress functions suitable for the problem are found and the evaluation of the J-integral for the theoretical singular crack tip fields yields energy release rates. The present results are simpler than the existing ones and can be reduced to the well known solutions in special cases. Examples of extracting stress intensity factors from the finite element solution using the present results are given for the dynamically growing crack problem of orthotropic materials.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2002.10a
• /
• pp.820-823
• /
• 2002

The asymptotic problem of a semi-infinite conducting crack parallel to the poling direction in ferroelectric ceramics subjected to electric fields is analyzed. The main mechanism for the conducting crack growth behavior is thought to be ferroelectric domain switching leading to the development of a process zone around the crack. The shape and size of the switching zone is shown to depend strongly on the relative magnitude on the ratio of the coercive electric field to the yield electric field. It is shown that the crack growth can be either enhanced or retarded depending on the ratio of the coercive electric field to yield electric field.