• Structural Engineering and Mechanics
• /
• v.66 no.3
• /
• pp.329-341
• /
• 2018

This paper considers the smooth receding contact problem between a homogeneous half-plane and a composite laminate composed of an inhomogeneously coated elastic layer. The inhomogeneity of the elastic modulus of the coating is approximated by an exponential function along the thickness dimension. The three-component structure is pressed together by either a concentrated force or uniform pressures applied at the top surface of the composite laminate. Both semianalytical and finite element analysis are performed to solve for the extent of contact and the contact pressure. In the semianalytical formulation, Fourier integral transformation of governing equations and boundary conditions leads to a singular integral equation of Cauchy-type, which can be numerically integrated by Gauss-Chebyshev quadrature to a desired degree of accuracy. In the finite element modeling, the functionally graded coating is divided into homogeneous sublayers and the shear modulus of each sublayer is assigned at its lower boundary following the predefined exponential variation. In postprocessing, the stresses of any node belonging to sublayer interfaces are averaged over its surrounding elements. The results obtained from the semianalytical analysis are successfully validated against literature results and those of the finite element modeling. Extensive parametric studies suggest the practicability of optimizing the receding contact peak stress and the extent of contact in multilayered structures by the introduction of functionally graded coatings.

• Structural Engineering and Mechanics
• /
• v.67 no.3
• /
• pp.277-281
• /
• 2018

The analytical solution of two functionally graded layers with Volterra type screw dislocation is investigated under anti-plane shear impact loading. The energy dissipation of FGM layers is modeled by viscous damping and the properties of the materials are assumed to change exponentially along the thickness of the layers. In this study, the rate of gradual change ofshear moduli, mass density and damping constant are assumed to be same. At first, the stress fields in the interface of the FGM layers are derived by using a single dislocation. Then, by determining a distributed dislocation density on the crack surface and by using the Fourier and Laplace integral transforms, the problem are reduce to a system ofsingular integral equations with simple Cauchy kernel. The dynamic stress intensity factors are determined by numerical Laplace inversion and the distributed dislocation technique. Finally, various examples are provided to investigate the effects of the geometrical parameters, material properties, viscous damping and cracks configuration on the dynamic fracture behavior of the interacting cracks.

• Structural Engineering and Mechanics
• /
• v.64 no.1
• /
• pp.71-79
• /
• 2017

In this study, interaction of several interface cracks located between a functionally graded material (FGM) layer and an elastic layer under anti-plane deformation based on the distributed dislocation technique (DDT) is analyzed. The variation of the shear modulus of the functionally graded coating is modeled by an exponential and linear function along the thickness of the layer. The complex Fourier transform is applied to governing equation to derive a system of singular integral equations with Cauchy type kernel. These equations are solved by a numerical method to obtain the stress intensity factors (SIFs) at the crack tips. The effects of non-homogeneity parameters for exponentially and linearly form of shear modulus, the thickness of the layers and the length of crack on the SIFs for several interface cracks are investigated. The results reveal that the magnitude of SIFs decrease with increasing of FG parameter and thickness of FGM layer. The values of SIFs for FGM layer with exponential form is less than the linear form.

• Journal of the Korean Society for Precision Engineering
• /
• v.18 no.12
• /
• pp.95-103
• /
• 2001

The buckling of an orthotropic layer bonded to an orthotropic half-space with an interface crack subjected to compressive load under plane strain is analyzed. General solution to the stability equations describing the buckling behavior of both the layer and the half-space is expressed in terms of displacement functions. The displacement functions are represented by the solution of Cauchy-type singular integral equations, which are numerically solved. Numerical results of the critical buckling loads are presented fur various geometric parameters and material properties of both the layer and half-space.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.4
• /
• pp.1275-1289
• /
• 1996

The plane elasticity problem of collinear cracks in a layered medium is investigated. The medium is modeled as bonded structure constituted from a surface layer and a semi-infinite substrate. Along the bond line between the two dissimilar homegeneous constituents, it is assumed that as interfacial zone having the functionally graded, nonhomogeneous elastic modulus exists. The layered medium contains three collinear cracks, one in each constituent material oriented perpendicular to the nominal interfaces. The stiffness matrix formulation is utilized and a set of homogeneous conditions relevant to the given problem is readily satisfied. The proposed mixed boundary value problem is then represented in the form of a system of integral equations with Cauchy-type singular kernels. The stress intensity factors are defined from the crack-tip stress fields possessing the standard square-root singular behavior. The resulting values of stress intensity factors mainly address the interactions among the cracks for various crack sizes and material combinations.

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