• Structural Engineering and Mechanics
• /
• v.78 no.5
• /
• pp.585-597
• /
• 2021

This paper presents a comparative study of analytical method, finite element method (FEM) and Multilayer Perceptron (MLP) for analysis of a contact problem. The problem consists of a functionally graded (FG) layer resting on a half plane and pressed with distributed load from the top. Firstly, analytical solution of the problem is obtained by using theory of elasticity and integral transform techniques. The problem is reduced a system of integral equation in which the contact pressure are unknown functions. The numerical solution of the integral equation was carried out with Gauss-Jacobi integration formulation. Secondly, finite element model of the problem is constituted using ANSYS software and the two-dimensional analysis of the problem is carried out. The results show that contact areas and the contact stresses obtained from FEM provide boundary conditions of the problem as well as analytical results. Thirdly, the contact problem has been extended based on the MLP. The MLP with three-layer was used to calculate the contact distances. Material properties and loading states were created by giving examples of different values were used at the training and test stages of MLP. Program code was rewritten in C++. As a result, average deviation values such as 0.375 and 1.465 was obtained for FEM and MLP respectively. The contact areas and contact stresses obtained from FEM and MLP are very close to results obtained from analytical method. Finally, this study provides evidence that there is a good agreement between three methods and the stiffness parameters has an important effect on the contact stresses and contact areas.

• Structural Engineering and Mechanics
• /
• v.45 no.4
• /
• pp.495-519
• /
• 2013

This paper investigates properties of integral operators in complex variable boundary integral equation in plane elasticity, which is derived from the Somigliana identity in the complex variable form. The generalized Sokhotski-Plemelj's formulae are used to obtain the BIE in complex variable. The properties of some integral operators in the interior problem are studied in detail. The Neumann and Dirichlet problems are analyzed. The prior condition for solution is studied. The solvability of the formulated problems is addressed. Similar analysis is carried out for the exterior problem. It is found that the properties of some integral operators in the exterior boundary value problem (BVP) are quite different from their counterparts in the interior BVP.

• Structural Engineering and Mechanics
• /
• v.25 no.4
• /
• pp.383-403
• /
• 2007

The contact problem for an elastic layer resting on an elastic half plane is considered according to the theory of elasticity with integral transformation technique. External loads P and Q are transmitted to the layer by means of two dissimilar rigid flat punches. Widths of punches are different and the thickness of the layer is h. All surfaces are frictionless and it is assumed that the layer is subjected to uniform vertical body force due to effect of gravity. The contact along the interface between elastic layer and half plane will be continuous, if the value of load factor, ${\lambda}$, is less than a critical value, ${\lambda}_{cr}$. However, if tensile tractions are not allowed on the interface, for ${\lambda}$ > ${\lambda}_{cr}$ the layer separates from the interface along a certain finite region. First the continuous contact problem is reduced to singular integral equations and solved numerically using appropriate Gauss-Chebyshev integration formulas. Initial separation loads, ${\lambda}_{cr}$, initial separation points, $x_{cr}$, are determined. Also the required distance between the punches to avoid any separation between the punches and the layer is studied and the limit distance between punches that ends interaction of punches, is investigated. Then discontinuous contact problem is formulated in terms of singular integral equations. The numerical results for initial and end points of the separation region, displacements of the region and the contact stress distribution along the interface between elastic layer and half plane is determined for various dimensionless quantities.

• Structural Engineering and Mechanics
• /
• v.26 no.5
• /
• pp.591-615
• /
• 2007

An improvement is introduced to solve the plane problems of linear elasticity by reciprocal theorem for orthotropic materials. This method gives an integral equation with complex kernels which will be solved numerically. An artificial boundary is defined to eliminate the singularities and also an algorithm is introduced to calculate multi-valued complex functions which belonged to the kernels of the integral equation. The chosen sample problem is a plate, having a circular or elliptical hole, stretched by the forces parallel to one of the principal directions of the material. Results are compatible with the solutions given by Lekhnitskii for an infinite plane. Five different orthotropic materials are considered. Stress distributions have been calculated inside and on the boundary. There is no boundary layer effect. For comparison, some sample problems are also solved by finite element method and to check the accuracy of the presented method, two sample problems are also solved for infinite plate.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.24 no.1
• /
• pp.93-102
• /
• 2020

The concept of temperature distribution in inhomogeneous semi-infinite solids is examined by making use of direct integration method. The analysis is done on the solution of the in-plane steady state heat conduction problem under certain boundary conditions. The method of direct integration has been employed, which is then reduced to Volterra integral equation of second kind, produces the explicit form analytical solution. Using resolvent- kernel algorithm, the governing equation is solved to get present solution. The temperature distribution obtained and calculated numerically and the relation with distribution of heat flux generated by internal heat source is shown graphically.

• Structural Engineering and Mechanics
• /
• v.51 no.6
• /
• pp.923-937
• /
• 2014

The rational finite element method is different from the standard finite element method, which is constructed using basic solutions of the governing differential equations as interpolation functions in the elements. Therefore, it is superior to the isoparametric approach because of its obvious physical meaning and accuracy; it has successfully been applied to the isotropic elasticity problem. In this paper, the formulation of rational finite elements for plane orthotropic elasticity problems is deduced. This method is formulated directly in the physical domain with full consideration of the requirements of the patch test. Based on the number of element nodes and the interpolation functions, different approaches are applied with complete polynomial interpolation functions. Then, two special stiffness matrixes of elements with four and five nodes are deduced as a representative application. In addition, some typical numerical examples are considered to evaluate the performance of the elements. The numerical results demonstrate that the present method has a high level of accuracy and is an effective technique for solving plane orthotropic elasticity problems.

• Advances in Computational Design
• /
• v.9 no.1
• /
• pp.73-80
• /
• 2024

The aim of this paper is to remove the rigid body motion in the interior boundary value problem (BVP) of plane elasticity by solving the interior and exterior BVPs simultaneously. First, we formulate the interior and exterior BVPs simultaneously. The tractions applied on the contour in two problems are the same. After adding and subtracting the two boundary integral equations (BIEs), we will obtain a couple of BIEs. In the coupled BIEs, the properties of relevant integral operators are modified, and those integral operators are generally invertible. Finally, a unique solution for boundary displacement of interior region can be obtained.

• Structural Engineering and Mechanics
• /
• v.9 no.2
• /
• pp.209-214
• /
• 2000

The main purpose of this paper is to develop the results introduced in Artan (1996) and to find a general nonlocal linear elastic solution for Boussinesq problem. The general nonlocal solution given Artan (1996) is valid only when the distance to the boundary is greater than one atomic measure. The nonlocal stress field presented in this paper is valid for the whole half plane.

• Structural Engineering and Mechanics
• /
• v.64 no.4
• /
• pp.505-513
• /
• 2017

In this paper, a receding contact problem for an elastic layer resting on a half plane is considered. The layer is pressed by two rectangular stamps placed symmetrically. It is assumed that the contact surfaces are frictionless and only compressive traction can be transmitted through the contact surfaces. In addition the effect of body forces is neglected. Firstly, the problem is solved analytically based on theory of elasticity. In this solution, the problem is reduced into a system of singular integral equations in which half contact length and contact pressures are unknowns using boundary conditions and integral transform techniques. This system is solved numerically using Gauss-Jacobi integral formulation. Secondly, two dimensional finite element analysis of the problem is carried out using ANSYS. The dimensionless quantities for the contact length and the contact pressures are calculated under various stamp size, stamp position and material properties using both solutions. The analytic results are verified by comparison with finite element results.

• Structural Engineering and Mechanics
• /
• v.84 no.1
• /
• pp.35-48
• /
• 2022

This paper investigates the artificial neural network (ANN) to predict the dimensionless parameters for contact pressures and contact lengths under the rigid punch, the initial separation loads, and the initial separation distances of a contact problem. The problem consisted of two elastic infinitely layers (EL) loaded by means of a rigid cylindrical punch and resting on a half-infinite plane (HP). Firstly, the problem was formulated and solved theoretically using the Theory of Elasticity (ET). Secondly, the contact problem was extended based on the ANN. External load, the radius of punch, layer heights, and material properties were created by giving examples of different values used at the training and test stages of ANN. Finally, the accuracy of the trained neural networks for the case was tested using 134 new data, generated via ET solutions to determine the best network model. ANN results were compared with ET results, and well agreements were achieved.