Advances in Mathematical Physics
Volume 2016 (2016), Article ID 3845362, 12 pages
http://dx.doi.org/10.1155/2016/3845362

Research Article

The Approximate Solution of Some Plane Boundary Value Problems of the Moment Theory of Elasticity

Roman Janjgava

Ilia Vekua Institute of Applied Mathematics, Ivane Javakhishvili Tbilisi State University, University Street 2, 0186 Tbilisi, Georgia

Received 14 August 2015; Revised 19 January 2016; Accepted 31 January 2016

Academic Editor: Mikhail Panfilov

Copyright © 2016 Roman Janjgava. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We consider a two-dimensional system of differential equations of the moment theory of elasticity. The general solution of this system is represented by two arbitrary harmonic functions and solution of the Helmholtz equation. Based on the general solution, an algorithm of constructing approximate solutions of boundary value problems is developed. Using the proposed method, the approximate solutions of some problems on stress concentration on the contours of holes are constructed. The values of stress concentration coefficients obtained in the case of moment elasticity and for the classical elastic medium are compared. In the final part of the paper, we construct the approximate solution of a nonlocal problem whose exact solution is already known and compare our approximate solution with the exact one. Supposedly, the proposed method makes it possible to construct approximate solutions of quite a wide class of boundary value problems.

1. Introduction

An elastic medium which at every point is characterized not only by the displacement vector but also by the rotation vector is called the Cosserat medium after the Cosserat brothers who considered the deformation of such a medium for the first time as far back as 1909 [1]. The corresponding theory is often called the moment or asymmetrical theory of elasticity because its deformation and stress tensors are not symmetrical. The development of this theory along various lines attained its peak in the sixties–seventies of the last century [2–12]. Though, in subsequent years, more general models of a micropolar elastic medium were proposed and studied [13–15], the investigation of problems related to the Cosserat theory is still going on [16–36]. The popularity of the theory can possibly be explained by its physical clarity and a relative mathematical simplicity.

In the present paper, we consider the static case of plane deformation of the elastic Cosserat medium. The corresponding basic equilibrium equations are given in Section 2. For constructing approximate solutions of boundary value problems, a semianalytical method is proposed, which is based on using the representation of the general solution by solutions of simpler equations of mathematical physics.

In Section 3 we construct the general solution of the corresponding system of equations using two harmonic functions and a solution of the Helmholtz equation. As is known, the method of representation of the general solution by the well-studied functions of mathematical physics is widely applied to the construction of analytical (exact) solutions of boundary value problems, but one succeeds in obtaining analytical solutions of boundary value problems only for some particular types of boundary conditions when the considered domain has a regular configuration.

Section 4 deals with the application of the obtained general solution to the construction of approximate solutions of boundary value problems, which allows us to consider a wide class of boundary value problems without imposing special restrictions on the domain geometry. When the functions appearing in the general solution are expanded into series in Green’s functions and the boundary conditions are satisfied at the points of the domain boundary, we obtain a system of linear algebraic equations for the unknown coefficients of the series. After solving this system and again using the general representation of the solution, we can write an analytical approximate solution which satisfies the basic system of equations up to the domain boundary. The described algorithm is similar to the method of fundamental solutions [37–39] and the boundary elements method [40–44]. However it is simpler than the method of fundamental solutions for which the integral equation is considered on some contour that encircles the domain [39]. As different from the boundary elements method, the boundary does not have to be approximated by boundary elements and there is no need to use local coordinate systems [40], which also simplifies the construction of an approximate solution. But the main difference of our approach from the above-mentioned two methods is that in our case it becomes possible not to use Green’s functions at all. At the end of Section 4, the harmonic functions and the Helmholtz function in the general solution are represented as finite series obtained as a result of the separation of variables in a polar system of coordinates. Further the Cartesian coordinates of the points on the boundary are substituted in the representations of the functions whose values are given on the boundary and the obtained expressions are equated to the respective values of the boundary conditions at the same points. As a result we obtain a system of linear algebraic equations and after solving it we can easily write an approximate solution of the considered boundary value problem.

In Section 5, we construct by the proposed method the approximate solutions of the boundary value problems of stress concentration in the rectangular domains with circular holes. The coefficients of stress concentration along the hole contour are calculated and we compare the stress concentration values obtained in the case of moment elasticity and for the classical elastic medium. For the classical elastic medium we also construct the approximate solution of the nonlocal problem with the known exact solution and compare our approximate solution with the exact one.

2. Basic Relations of the Plane Moment Theory of Elasticity

Let be the Cartesian coordinate system with unit vectors . Consider the plane deformed state of the moment-elastic medium which is parallel to the -plane. In that case, the displacement vector and the rotation vector at the points of the medium have the form

Homogeneous equations of static equilibrium in terms of displacement and rotation components are written as [10]

where are the Lamé constants; are the constants characterizing the microstructure of the considered medium; and ,  , and 

Relations of Hooke’s law for the moment-elastic medium have the form

where are the stress tensor components and are the moment stress tensor components. All other stress tensor and moment stress tensor components are assumed in the case of plane deformation to be equal to zero.

3. The General Solution of System (2)

The general solution of system (2) can be represented by two arbitrary harmonic functions and an arbitrary solution of the Helmholtz equation.

We introduce the following notation:

Using (4), the first two equations of system (2) can be written as follows:

From system (5) we see that and are the self-conjugate harmonic functions

and represent them in the form

where and are arbitrary harmonic functions and is any nonzero real constant.

Remark. In representations (7) only one harmonic function could have been used, but we have introduced two harmonic functions to obtain the general solution which is more convenient for our purpose.

With (7) taken into consideration, from relations (4) we have

Let us substitute the second formula (8) in the third equation of system (2). After simple transformations we obtain the equation

In view of the fact that is a harmonic function, the general solution of (9) can be written in the form

where is an arbitrary solution of the Helmholtz equation

and is any nonzero real constant.

Using representation (10) and keeping in mind that is the solution of (11), we write system (8) as follows:

The second equation (12) is satisfied identically if we take

By substituting formulas (13) in the first equation of system (12), we obtain the Poisson equation satisfied by function :

Since and are harmonic functions, we can easily write the general solution of the latter equation:

where is an arbitrary harmonic function.

The substitution of expression (15) in formula (13) gives

where

Obviously, and are arbitrary harmonic functions. For simplicity, we denote them by and

It is convenient to fix arbitrary constants and as follows:

Then formulas (16) and (10) take the form

System (19) is the representation of the general solution of system (2) by two arbitrary harmonic functions and and an arbitrary solution of the Helmholtz equation.

Let us now express stresses and moment stresses through functions , and . To this end, formulas (19) are substituted in relations (3) to obtain

Let us assume that and are mutually normal unit vectors such that

and is the angle between vector and the positive direction of the -axis. Then, by virtue of the transformation formulas of first- and second-rank tensor components, the following formulas hold true:

Substituting formulas (19) and (20) in relations (22) we can express normal and tangent displacements as well as stresses and moment stresses acting on the area element of arbitrary orientation through functions ,  , and

4. Use of the General Solution for the Construction of Approximate Solutions of Boundary Value Problems

In this section it is shown how the general solution derived in the preceding section can be used for constructing approximate solutions of boundary value problems. The technique to be described here may be called a semianalytical method in view of the fact that approximate solutions are written in analytical form. In this sense, it is analogous to the method of fundamental solutions [37–39] and the boundary elements method [40–44], but there are differences which have been mentioned in Introduction.

So, we start the description of the algorithm of constructing the approximate solution of a boundary value problem when Green’s functions are used.

In the general solutions (19) and (20), the harmonic functions and are taken for each index ( is some natural number) in the following manner:

and solution of the Helmholtz equation is represented as

where is a modified Bessel function of second kind or a Macdonald function (of zeroth order) which tends to infinity at point and vanishes at infinity [45] and and  ,  , are the real coefficients to be defined.

By formula (23) the first-order partial derivatives of the harmonic functions and have the form

and the partial derivatives of the metaharmonic functions are written as

The second derivatives of these functions appear in formulas (20):

If we substitute the corresponding expressions (23)–(26) in formulas (19), then for each index we obtain solution of system (2) satisfying this system everywhere except the origin where it has a singularity:

The values of stresses and moment stresses are obtained if representations (25)–(27) are substituted in the respective formulas (20):

After that, each th function of (28) and (29) is shifted by value (Figure 1). For this, in formulas (28) and (29) and also in the expressions for stresses and moment stresses variables and are replaced, respectively, by values and . Then functions , as well as also functions , have singularities at points

Figure 1: Shift of a singular point.

Let us consider the following sums:

For our algorithm approximate solutions of boundary value problems are sought in form (30). Functions (30) satisfy system of (2) everywhere except points , , where they have a singularity. The unknown coefficients ,  , figuring in formulas (28) and (29), and therefore in formulas (30) too, have to be chosen so that the boundary conditions be satisfied.

For simplicity, let us consider the algorithm of construction of an approximate solution using the example of a boundary value problem for a finite domain with displacements and rotation given on the boundary. Let be a simply or multiply connected domain with a sufficiently smooth boundary having almost everywhere the external normal . In domain , we seek such a solution of system (2) that satisfies the following boundary conditions on boundary :

where , and are the functions given on the boundary.

In the first place, on the boundary we more or less uniformly distribute points with coordinates . As to points , we arrange them on some contour that lies around the external boundary of and also on the contours within the internal boundaries if the domain is multiply connected as shown in Figure 2.

Figure 2: Domain with singular points lying on some external contour and points lying on the boundary contours.

An approximate solution is sought in the form

We substitute the coordinates of points marked on the boundary and the values of angles in formulas (32) and equate the resulting expressions to the respective values of the boundary functions at these points. As a result we obtain the following system of the equations

But according to formulas (28), the latter system is a system of linear algebraic equations with unknowns ,  

where and  ,  , are the already known coefficients. For example, has the form

After solving system (34), we substitute the obtained coefficient values in the respective expressions (30) and find all the desired displacement and rotation components, as well as all stress and moment stress components.

In an analogous manner we can construct approximate solutions of boundary value problems when stresses and moment stresses are given on the domain boundary. We can also consider all kinds of mixed boundary value problems.

In the algorithm described above we have used Green’s functions for the construction of the approximate solution. But, as has already been mentioned in Introduction, the use of Green’s functions is not obligatory at all. The harmonic functions and a solution of the Helmholtz equation can be represented as finite series obtained as a result of the separation of variables in a polar system of coordinates but, at the same time, boundary value problems are solved in a Cartesian coordinate system.

Let us give a brief description of the latter approach to the construction of an approximate solution of a boundary value problem. It is assumed for the sake of simplicity that we consider a doubly connected domain bounded by the simple closed contours , of which encircles , while lies around the origin.

The harmonic functions and, in the general solutions (19) and (20) are represented as follows:

where and

Function is written in the form

where and are, respectively, a modified Bessel and a Macdonald function of order [45].

Substituting formulas (36) and (38) in the general solutions (19) and (20), we obtain the analytical type of the desired approximate solution. However this solution contains the unknown coefficients ,  ,  

Further, we more or less uniformly distribute points with coordinates on contours and (Figure 3). In order that the given boundary conditions be satisfied at these points (Figure 3), we substitute the coordinates of the boundary points ,  , in the respective formulas of the approximate solution and equate the resulting expressions to the known values of the boundary conditions at these points. Thus, for coefficients ,  ,  , we obtain a system of linear algebraic equations with unknowns. By solving this system we define the values of all the unknown coefficients in the formulas for an approximate solution and thereby accomplish the construction of the approximate solution.

Figure 3: The considered doubly connected domain .

In Section 5, by the above technique not requiring the use of Green’s functions we solve a nonlocal problem of the Bitsadze-Samarskii type [46–49]. We want to show by the above example that the proposed method can be used for constructing approximate solutions not only of classical boundary value problems but also of some nonclassical ones.

As to the accuracy of approximate solutions constructed by the proposed method, we can say the following. Since in our case an approximate solution of a boundary value problem has the analytical form, we can easily calculate the values of this solution on the domain boundary. If the result obtained is sufficiently close to the given boundary conditions, then provided that the problem is well posed, we can conclude that the constructed solution is sufficiently close to an exact solution since the equilibrium equations are satisfied exactly within the domain.

5. Examples

We take the following elastic constants of the Cosserat medium:  GPa,  GPa,  GPa, and  GN. These data correspond to the elastic characteristics of brass.

Example 1. Let domain be a square with a circular hole , where and (Figure 4).
Consider the so-called Kirsch problem [19, 45] for domain when the constant normal stresses are applied to the two opposite sides of the square, and the remaining part of the boundary is free from stresses (Figure 5). Thus we solve the following boundary value problem:

This problem is solved approximately by the proposed method both for the Cosserat medium and for the classical elastic medium.
We arrange points ,  , on the square containing domain (48 points) and on the circumference within the internal contour (24 points) as is shown in Figure 6. On the domain boundary we also mark 72 points at which the given boundary conditions are satisfied. In the general case, for the unknown coefficients we obtain linear algebraic equations with the same number of unknowns. If the problem is symmetric with respect to the coordinate axes, then the number of unknown coefficients and the number of equations can be reduced. In that case, in the equations we preserve only the coefficients corresponding to the first quadrant and express the remaining coefficients through them according to the formulas

After solving the resulting system we compared the boundary values ​​of the obtained solutions with the given boundary conditions. If absolute errors of order are regarded as small from the practical point of view (errors are especially noticeable near the angular points), then the resulting solution can be regarded as satisfactory. If a higher exactness is needed, it can be achieved by increasing the number of points on the boundary where the boundary conditions are satisfied or by a more appropriate choice of the contour on which the singular points are arranged. For some problems admitting an exact solution, we have managed to construct the approximate solutions with absolute errors of order not higher than

Figure 4: The considered domain .

Figure 5: Unilateral tension of domain .

Figure 6: Domain encircled by singular points.

In the case of the Kirsch problem it is of special interest to investigate the distribution of stress along the hole contour. A maximal value of stress both for the classical elastic medium and for the Cosserat medium is obtained, as expected, at . For the classical medium, and this value does not depend on the material (if the material is isotropic and homogenous). For the Cosserat medium, depends on the material and in our case its value is (Figure 7).

Figure 7: Graphs of the distribution of stress along the hole contour for the classical elastic medium and the Cosserat medium for domain .

Example 2. Now let us consider the Kirsch problem for domain ,   (Figure 8).
As expected, in this case the values of diminish both for the classical medium and for the Cosserat medium and the difference between them becomes even more insignificant (see Figure 9). We know from the literature that in the case of an infinite domain with a circular hole value for the Cosserat medium is smaller than the stress concentration coefficient for the classical elastic medium [20] which is equal to three. Therefore, as the ratio of , where is the side of the square and is the hole radius, increases to a certain value, values for the classical medium and for the Cosserat medium get closer to each other.
Finally, let us we consider the nonlocal problem for a rectangular domain for the classical elastic medium.

Figure 8: The considered domain .

Figure 9: Graphs of the distribution of stress along the hole contour for the classical medium and the Cosserat medium for domain .

Example 3. We consider domain (Figure 10). In domain , it is required to find a solution of system (2) (where ,  ,  , and  that satisfies the following conditions:

The exact solution of this problem is as follows:

The boundary contour of the considered domain is divided by points into 72 equal segments. 19 points are also distributed uniformly on a segment within the domain with given nonlocal conditions. An approximate solution is sought in the form

The stress tensor components and figuring in the initial condition are represented as finite series:

Next we substitute the coordinates of the points marked on the boundary and within the domain into the last four formulas and thus satisfy at these points the corresponding boundary and nonlocal conditions. As a result we obtain the system consisting of 144 linear algebraic equations and containing 144 unknowns ,  ,  ,  . After solving this system by the above formulas we can easily define the displacement vector and stress tensor components.

Figure 10: Domain in which the nonlocal problem of plane elasticity is solved.

In our study we used Maple 17 software. The numerical results are presented in Table 1, where and denote the approximate values of the displacement vector components.

Table 1: Numerical results for Example 3.

From the numerical results we see that the proposed technique gives the good approximate solution for the considered nonlocal mixed boundary value problem of plane elasticity.

6. Conclusion

In this paper, we consider the static case of plane deformation for the elastic Cosserat medium. The constructed general solution of the corresponding system of differential equations is represented by two harmonic functions and a solution of the Helmholtz equation. Based on this solution, we propose the algorithm for the approximate solution of a wide class of boundary value problems.

The same boundary value problems can be solved by the method of Green’s functions or boundary elements. Our technique, based on representing the solution through two harmonic functions and a solution of the Helmholtz equation, is an alternative method, different from them.

To show the validity of the proposed algorithm for moment elasticity, we construct the approximate solutions of boundary value problems of stress concentration in a rectangular domain with circular holes. The stress concentration coefficients are calculated along the hole contour and their values are compared with those for the classical elastic case.

For the classical elastic medium we also obtain the approximate solution of a nonlocal problem and compare it with the known exact solution.

We believe that the proposed algorithm can be applied for constructing an approximate solution of a wide class of boundary value problems for systems of partial differential equations whose general solutions can be represented by means of the solutions of relatively simple equations of mathematical physics.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

References

E. Cosserat and F. Cosserat, Theorie des Corps Deformables, Hermann, Paris, France, 1909.

C. Truesdell and R. A. Toupin, “The classical field theories,” in Handbuch der Physik, Bd. III/1, S. Flügge, Ed., Springer, Berlin, Germany, 1960.

G. Grioli, “Elasticità assimetrica,” Annali di Matematica Pura ed Applicata, Serie IV, vol. 50, pp. 389–417, 1960.

E. Kuvshinskii and E. L. Aero, “Continuum theory of asymmetric elasticity. Equilibrium of isotropic body,” Fizika Tverdogo Tela, vol. 6, no. 9, pp. 2689–2699, 1964 (Russian).

R. D. Mindlin, “Influence of couple-stress on stress concentrations,” Experimental Mechanics, vol. 3, no. 1, pp. 1–7, 1963.

V. A. Pal'mov, “Fundamental equations of the theory of asymmetric elasticity,” Journal of Applied Mathematics and Mechanics, vol. 28, pp. 496–505, 1964 (Russian).

A. C. Eringen and E. S. Suhubi, “Nonlinear theory of simple micro-elastic solids—I,” International Journal of Engineering Science, vol. 2, no. 2, pp. 189–203, 1964.

W. Novacki, “Couple stresses in the theory of thermoelasticity,” Bulletin of the Polish Academy of Sciences: Technical Sciences, vol. 14, no. 8, 1966.

A. E. Green and P. M. Naghdi, “The linear theory of an elastic Cosserat plate,” Mathematical Proceedings of the Cambridge Philosophical Society, vol. 63, no. 2, pp. 537–550, 1967.

W. Novacki, “On the completeness of stress functions in asymmetric elasticity,” Bulletin of the Polish Academy of Sciences: Technical Sciences, vol. 14, no. 7, 1968.

V. D. Kupradze, T. G. Gegelia, M. O. Basheleishvili, and T. V. Burchuladze, Three-Dimensional Problems of the Mathematical Theory of Elasticity and Thermoelasticity, North-Holland, Amsterdam, Netherlands; Nauka, Moscow, Russia, 1979.

W. Novacki, Theory of Asymmetric Elasticity, Polish Scientific, Warsaw, Poland, 1986.

D. Natroshvili, R. Gachechiladze, A. Gachechiladze, and I. G. Stratis, “Transmission problems in the theory of elastic hemitropic materials,” Applicable Analysis, vol. 86, no. 12, pp. 1463–1508, 2007.

D. Natroshvili, L. Giorgashvili, and I. G. Stratis, “Mathematical problems of the theory of elasticity of chiral materials,” Applied Mathematics, Informatics and Mechanics, vol. 8, no. 1, pp. 47–103, 127, 2003.

R. Janjgava and B. Gulua, “About plane theory for hemitropic elastic materials,” Reports of Enlarged Session of the Seminar of VIAM, vol. 19, no. 1, pp. 38–42, 2004.

T. V. Meunargia, Development of a Method of I. N. Vekua for Problems of the Three-Dimensional Moment Theory of Elasticity, TSU Publishing, Tbilisi, Georgia, 1987 (Russian).

P. Schiavone, “Fundamental sequences of functions in the approximation of solutions to exterior problems in the bending of thin micropolar plates,” Applicable Analysis, vol. 35, no. 1–4, pp. 263–274, 1990.

J. D. Kaplunov and H. Lippmann, “Elastic-plastic torsion of a Cosserat-type rod,” Acta Mechanica, vol. 113, no. 1–4, pp. 53–62, 1995.

É. A. Bulanov, “On the moment theory of elasticity. Plane deformation. Part 3,” Strength of Materials, vol. 30, no. 5, pp. 516–521, 1998.

М. А. Kulesh, V. P. Мatveenko, and I. N. Shardakov, “Construction and analysis of exact analytical solution of kirsch problem within the cosserat continuum and pseudo-continuum,” Journal of Applied Mechanics and Technical Physics, vol. 42, no. 4, pp. 145–154, 2001 (Russian).

E. Providas and M. A. Kattis, “Finite element method in plane Cosserat elasticity,” Computers & Structures, vol. 80, no. 27–30, pp. 2059–2069, 2002.

M. A. Kulesh, V. P. Matveenko, and I. N. Shardakov, “Parametric analysis of analytical solutions to one- and two-dimensional problems in couple-stress theory of elasticity,” ZAMM: Zeitschrift für Angewandte Mathematik und Mechanik, vol. 83, no. 4, pp. 238–248, 2003.

A. Zeghadi, S. Forest, A. Gourgues, and O. Bouaziz, “Cosserat continuum modelling of grain size effects in metal polycrystals,” PAMM, vol. 5, no. 1, pp. 79–82, 2005.

H. W. Zhang, H. Wang, and G. Z. Liu, “Quadrilateral isoparametric finite elements for plane elastic Cosserat bodies,” Acta Mechanica Sinica, vol. 21, no. 4, pp. 388–394, 2005.

S. Potapenko, “A generalized fourier approximation in anti-plane Cosserat elasticity,” Journal of Elasticity, vol. 81, no. 2, pp. 159–177, 2005.

M. Bîrsan, “On Saint-Venant's principle in the theory of Cosserat elastic shells,” International Journal of Engineering Science, vol. 45, no. 2–8, pp. 187–198, 2007.

M. Bîrsan, “Inequalities of Korn's type and existence results in the theory of Cosserat elastic shells,” Journal of Elasticity, vol. 90, no. 3, pp. 227–239, 2008.

S. O. Sarkisyan, “Boundary value problems of the asymmetric theory of elasticity for thin plates,” Journal of Applied Mathematics and Mechanics, vol. 72, no. 1, pp. 129–147, 2008.

J. Jeong, H. Ramezani, I. Munch, and P. Neff, “A numerical study for linear isotropic Cosserat elasticity with conformally invariant curvature,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 89, no. 7, pp. 552–569, 2009.

H. Altenbach and V. A. Eremeyev, “On the linear theory of micropolar plates,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 89, no. 4, pp. 242–256, 2009.

A. Riahi and J. H. Curran, “Full 3D finite element Cosserat formulation with application in layered structures,” Applied Mathematical Modelling, vol. 33, no. 8, pp. 3450–3464, 2009.

J. Altenbach, H. Altenbach, and V. A. Eremeyev, “On generalized Cosserat-type theories of plates and shells: a short review and bibliography,” Archive of Applied Mechanics, vol. 80, no. 1, pp. 73–92, 2010.

E. F. Grekova, “Linear reduced cosserat medium with spherical tensor of inertia, where rotations are not observed in experiment,” Mechanics of Solids, vol. 47, no. 5, pp. 538–543, 2012.

V. V. Korepanov, V. P. Matveenko, and I. N. Shardakov, “Finite element analysis of two- and three-dimensional static problems in the asymmetric theory of elasticity as a basis for the design of experiments,” Acta Mechanica, vol. 223, no. 8, pp. 1739–1750, 2012.

S. De Cicco and D. Ieşan, “A theory of chiral Cosserat elastic plates,” Journal of Elasticity, vol. 111, no. 2, pp. 245–263, 2013.

P. Neff, M. Bîrsan, and F. Osterbrink, “Existence theorem for geometrically nonlinear Cosserat micropolar model under uniform convexity requirements,” Journal of Elasticity, vol. 121, no. 1, pp. 119–141, 2015.

V. D. Kupradze and M. A. Aleksidze, “The method of functional equations for the approximate solution of certain boundary value problems,” USSR Computational Mathematics and Mathematical Physics, vol. 4, no. 4, pp. 82–126, 1964.

M. A. Aleksidze, Solution of Boundary Value Problems by Expansion in Non-Orthogonal Functions, Nauka, Moscow, Russia, 1978 (Russian).

M. A. Aleksidze, Fundamental Functions in Approximate Solutions of Boundary Value Problems, Nauka, Moscow, Russia, 1991 (Russian).

S. L. Crouch and A. M. Starfield, Boundary Element Methods in Solid Mechanics, George Allen & Unwin, London, UK, 1983.

C. A. Brebbia, J. C. F. Telles, and L. C. Wrobel, Boundary Element Teqhniques. Theory and Applications in Engineering, Springer, New York, NY, USA, 1984.

P. K. Banerjee, The Boundary Element Methods in Engineering, McGraw-Hill College, London, UK, 1994.

G. Beer, I. Smith, and C. Duenser, The Boundary Element Method with Programming: For Engineers and Scientists, Springer, Vienna, Austria, 2008.

J. T. Katsikadelis, Boundary Elements: Theory and Applications, Elsevier science, Oxford, UK, 2002.

N. N. Lebedev, Special Functions and Their Applications, Prentice-Hall, 1965.

A. V. Bitsadze and A. A. Samarskii, “On some simple generalizations of linear elliptic problems,” Doklady Akademii Nauk SSSR, vol. 185, no. 4, pp. 739–740, 1969 (Russian).

A. V. Bitsadze, “On the theory of nonlocal boundary value problems,” Doklady Akademii Nauk SSSR, vol. 277, no. 1, pp. 17–19, 1984 (Russian).

D. G. Gordeziani, “On one method for solving the boundary value problem of Bitsadze-Samarskii,” Seminar of Institute of Applied Mathematics of Tbilisi State University, Annotation of Reports, vol. 2, pp. 39–40, 1970 (Russian).

D. G. Gordeziani, On Methods for Solving a Class of Nonlocal Boundary Value Problems, TSU Press, Tbilisi, Georgia, 1981 (Russian).