• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.6 s.177
• /
• pp.1557-1564
• /
• 2000

When the half infinite crack in the orthotropic material strip with a large anisotropic ratio(E11>>E22) propagates with constant velocity, dynamic stress component $\sigma$y occurre d along the $\chi$ axis is derived by using the Fourier transformation and Wiener-Hopf technique, and the dynamic stress intensity factor is derived. The dynamic stress intensity factor depends on a crack velocity, mechanical properties and specimen hight. The normalized dynamic stress intensity factors approach the maximum values when normalized time(=Cs/a) is about 2. They have the constant values when the normalized time is greater than or equal to about 2, and decrease with increasing a/h(h: specimen hight, a: crack length) and the normalized crack propagation velocity( = c/Cs, Cs: shear wave velocity, c: crack propagation velocity).

• Tribology and Lubricants
• /
• v.38 no.3
• /
• pp.73-83
• /
• 2022

This paper is concerned with an analysis of a surface edge crack emanated from a sharp contact edge. For a geometrical model, a square wedge is in contact with a half plane whose materials are identical, and a surface perpendicular crack initiated from the contact edge exists in the half plane. To analyze this crack problem, it is necessary to evaluate the stress field on the crack line which are induced by the contact tractions and pseudo-dislocations that simulate the crack, using the Bueckner principle. In this Part I, the stress filed in the half plane due to the contact is re-summarized using an asymptotic analysis method, which has been published before by the author. Further focus is given to the stress field in the half plane due to a pseudo-edge dislocation, which will provide a stress solution due to a crack (i.e. a continuous distribution of edge dislocations) later, using the Burgers vector. Essential result of the present work is the corrective functions which modify the stress field of an infinite domain to apply for the present one which has free surfaces, and thus the infiniteness is no longer preserved. Numerical methods and coordinate normalization are used, which was developed for an edge crack problem, using the Gauss-Jacobi integration formula. The convergence of the corrective functions are investigated here. Features of the corrective functions and their application to a crack problem will be given in Part II.

• Journal of applied mathematics & informatics
• /
• v.8 no.1
• /
• pp.59-69
• /
• 2001

The paper deals with the interaction between three Griffith cracks propagating under antiplane shear stress at the interface of two dissimilar infinite elastic half-spaces. The Fourier transform technique is used to reduce the elastodynamic problem to the solution of a set of integral equations which has been solved by using the finite Hilbert transform technique and Cooke’s result. The analytical expressions for the stress intensity factors at the crack tips are obtained. Numerical values of the interaction efect have been computed for and results show that interaction effects are either shielding or amplification depending on the location of each crack with respect to other and crack tip spacing. AMS Mathematics Subject Classification : 73M25.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.7 s.178
• /
• pp.1863-1870
• /
• 2000

When an impact stress is applied on the external boundary of double cantilever beam of orthotropic material which crack length is greater than specimen hight and anistropic ratio is very high, dyna mic energy release rate is derived, and the relationship between dynamic energy release rate and crack propagating velocity is studied. Dynamic energy release rate to static energy release rate is decreased with increasment of crack propagating velocity. The relationships between dynamic energy release rate and vertical strain have a similar pattern with those between static energy release rate and vertical strain. When normalized time(Cstla) is greater than or equal to 2, dynamic energy release rate approaches to a constant value.

• Steel and Composite Structures
• /
• v.50 no.1
• /
• pp.77-87
• /
• 2024

In this study, the two-dimensional crack problem was investigated by using the finite element method (FEM)-based ANSYS package program and the artificial neural network (ANN)-based multilayer perceptron (MLP) method. For this purpose, a half-infinite functionally graded (FG) layer with a crack pressed through two rigid blocks was analyzed using FEM and ANN. Mass forces and friction were neglected in the solution. To control the validity of the crack problem model exercised, the acquired results were compared with a study in the literature. In addition, FEM and ANN results were checked using Root Mean Square Error (RMSE) and coefficient of determination (R²), and a well agreement was found. Numerical solutions were made considering different geometric parameters and material properties. The stress intensity factor (SIF) was examined for these values, and the results were presented. Consequently, it is concluded that the considered non-dimensional quantities have a noteworthy influence on the SIF. Also FEM and ANN can be logical alternative methods to time-consuming analytical solutions if used correctly.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.7 s.178
• /
• pp.1740-1747
• /
• 2000

Under repeated rolling, initial plastic deformation introduces residual stresses which render the steady cyclic state purely elastic. This is called the process of shakedown. Many studies have been done about the shakedown in semi-infinite half space using calculated Hertizian pressure. In this paper shakedown processes in a shaft are studied by finite element analyses of a two dimensional(plane strain) model with elastic-linear-kinematic-hardening-plastic material subjected to repeated, frictionless rolling contact. Symmetric and non-symmetric pressure distributions are obtained using a simplified model of the bearing-shaft assembly. The rolling contact is simulated by repeatedly translating both pressure distributions along the surface of the shaft. By the influence of the non-symmetric pressure, larger residual radial tensile stress is generated in the immediate subsurface layer, which may make a crack propagate and, the subsurface undergoes a zigzag plastic deformation during the shakedown process, which may lead to a crack initiation.

• Steel and Composite Structures
• /
• v.45 no.4
• /
• pp.501-511
• /
• 2022

In this study, a two-dimensional model of the contact problem has been examined using the finite element method (FEM) based software ANSYS and based on the multilayer perceptron (MLP), an artificial neural network (ANN). For this purpose, a functionally graded (FG) half-infinite layer (HIL) with a crack pressed by means of two rigid blocks has been solved using FEM. Mass forces and friction are neglected in the solution. Since the problem is analyzed for the plane state, the thickness along the z-axis direction is taken as a unit. To check the accuracy of the contact problem model the results are compared with a study in the literature. In addition, ANSYS and MLP results are compared using Root Mean Square Error (RMSE) and coefficient of determination (R2), and good agreement is found. Numerical solutions are made by considering different values of external load, the width of blocks, crack depth, and material properties. The stresses on the contact surfaces between the blocks and the FG HIL are examined for these values, and the results are presented. Consequently, it is concluded that the considered non-dimensional quantities have a noteworthy influence on the contact stress distributions, and also, FEM and ANN can be efficient alternative methods to time-consuming analytical solutions if used correctly.

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

• Tribology and Lubricants
• /
• v.38 no.3
• /
• pp.84-92
• /
• 2022

In Part I, developed was a method to obtain the stress field due to an edge dislocation that locates in an elastic half plane beneath the contact edge of an elastically similar square wedge. Essential result was the corrective functions which incorporate a traction free condition of the free surfaces. In the sequel to Part I, features of the corrective functions, F_kij,(k = x, y;i,j = x,y) are investigated in this Part II at first. It is found that F_xxx(ŷ) = F_xyx(ŷ) where ŷ = y/η and η being the location of an edge dislocation on the y axis. When compared with the corrective functions derived for the case of an edge dislocation at x = ξ, analogy is found when the indices of y and x are exchanged with each other as can be readily expected. The corrective functions are curve fitted by using the scatter data generated using a numerical technique. The algebraic form for the curve fitting is designed as F_kij(ŷ) = $\frac{1}{\hat{y}^{1-{\lambda}}I+yp}$$\sum_{q=0}^{m}{\left}$$\left[A_q\left(\frac{\hat{y}}{1+\hat{y}} \right)^q \right]$ where λ_I=0.5445, the eigenvalue of the adhesive complete contact problem introduced in Part I. To investigate the exponent of F_kij, i.e.(1 - λ_I) and p, Log|F_kij|(ŷ)-Log|(ŷ)| is plotted and investigated. All the coefficients and powers in the algebraic form of the corrective functions are obtained using Mathematica. Method of analyzing a surface perpendicular crack emanated from the complete contact edge is explained as an application of the curve-fitted corrective functions.