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

Authors: Roman Janjgava
Journal: Advances in Mathematical Physics
 (2016)

---

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


---

## Body

## 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 Section2. 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 Section3 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.Section4 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 Section5, 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

LetOxyz be the Cartesian coordinate system with unit vectors i,j,k. Consider the plane deformed state of the moment-elastic medium which is parallel to the Oxy-plane. In that case, the displacement vector u and the rotation vector w at the points of the medium have the form(1)u=ux,yi+vx,yj,w=ωx,yk.Homogeneous equations of static equilibrium in terms of displacement and rotation components are written as [10](2)μ+αΔu+λ+μ-α∂x∂xu+∂yv+2α∂yω=0,μ+αΔv+λ+μ-α∂y∂xu+∂yv-2α∂xω=0,ν+βΔω+2α∂xv-∂yu-4αω=0,where λ,μ are the Lamé constants; α,β,ν are the constants characterizing the microstructure of the considered medium; and ∂x≡∂/∂x,  ∂y≡∂/∂y, and Δ≡∂xx+∂yy.Relations of Hooke’s law for the moment-elastic medium have the form(3)σxx=λ∂xu+∂yv+2μ∂xu,σyy=λ∂xu+∂yv+2μ∂yv,σxy=μ+α∂xv+μ-α∂yu-2αω,σyx=μ+α∂yu+μ-α∂xv+2αω,σzz=λ∂xu+∂yv,μxz=ν+β∂xω,μyz=ν+β∂yω,μzx=ν-β∂xω,μzy=ν-β∂yω,where σxx,…,σyx are the stress tensor components and μxz,…,μzy 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:(4)θ=λ+2μ∂xu+∂yv,η=μ+α∂xv-∂yu-2αω.Using (4), the first two equations of system (2) can be written as follows:(5)∂xθ-∂yη=0,∂yθ+∂xη=0.From system (5) we see that θ and η are the self-conjugate harmonic functions(6)Δθ=0,Δη=0and represent them in the form(7)θ=a∂xφ+∂yψ,η=a∂xψ-∂yφ,where φ and ψ are arbitrary harmonic functions and a 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 (8)∂xu+∂yv=aλ+2μ∂xφ+∂yψ,∂xv-∂yu=aμ+α∂xψ-∂yφ+2αμ+αω.Let us substitute the second formula (8) in the third equation of system (2). After simple transformations we obtain the equation(9)Δω-4μαν+βμ+αω=-2αaν+βμ+α∂xψ-∂yφ.In view of the fact that φ is a harmonic function, the general solution of (9) can be written in the form (10)ω=bχ+a2μ∂xψ-∂yφ,where χ is an arbitrary solution of the Helmholtz equation(11)Δχ-4μαν+βμ+αχ=0and b is any nonzero real constant.Using representation (10) and keeping in mind that χ is the solution of (11), we write system (8) as follows:(12)∂xu+∂yv=aλ+2μ∂xφ+∂yψ,∂xv-∂yu=aμ∂xψ-∂yφ+ν+βb2μ∂xxχ+∂yyχ.The second equation (12) is satisfied identically if we take(13)u=∂xΦ+aμφ-ν+βb2μ∂yχ,v=∂yΦ+aμψ+ν+βb2μ∂xχ.By substituting formulas (13) in the first equation of system (12), we obtain the Poisson equation satisfied by function Φ: (14)ΔΦ=-λ+μaμλ+2μ∂xφ+∂yψ.Since φ and ψ are harmonic functions, we can easily write the general solution of the latter equation:(15)Φ=-λ+μa2μλ+2μxφ+∂xΨ+yψ+∂yΨ+aμΨ,where Ψ is an arbitrary harmonic function.The substitution of expression (15) in formula (13) gives(16)2μu=λ+3μaλ+2μφ∗-λ+μaλ+2μx∂xφ∗+y∂xψ∗-ν+βb∂yχ,2μv=λ+3μaλ+2μψ∗-λ+μaλ+2μx∂yφ∗+y∂yψ∗+ν+βb∂xχ,where(17)φ∗=φ+∂xΨ,ψ∗=ψ+∂yΨ.Obviously, φ∗ and ψ∗ are arbitrary harmonic functions. For simplicity, we denote them by φ and ψ.It is convenient to fix arbitrary constantsa and b as follows: (18)a=λ+2μλ+μ,b=1ν+β.Then formulas (16) and (10) take the form(19)2μu=λ+3μλ+μφ-x∂xφ-y∂xψ-∂yχ,2μv=λ+3μλ+μψ-y∂yψ-x∂yφ+∂xχ,2μω=2μν+βχ+λ+2μλ+μ∂xψ-∂yφ.System (19) is the representation of the general solution of system (2) by two arbitrary harmonic functions φ(x,y) and ψ(x,y) and an arbitrary solution χ(x,y) 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 (20)σxx=λ+2μλ+μ∂x-x∂xxφ+λλ+μ∂y-y∂xxψ-∂xyχ,σyy=λλ+μ∂x-x∂yyφ+λ+2μλ+μ∂y-y∂yyψ+∂xyχ,σxy=μλ+μ∂y-x∂xyφ+μλ+μ∂x-y∂xyψ-∂yyχ,σyx=μλ+μ∂y-x∂xyφ+μλ+μ∂x-y∂xyψ+∂xxχ,μxz=ν+βλ+2μ2μλ+μ∂xxψ-∂xyφ+∂xχ,μyz=ν+βλ+2μ2μλ+μ∂xyψ-∂yyφ+∂yχ,μzx=ν-βλ+2μ2μλ+μ∂xxψ-∂xyφ+ν-βν+β∂xχ,μzy=ν-βλ+2μ2μλ+μ∂xyψ-∂yyφ+ν-βν+β∂yχ.Let us assume that l and s are mutually normal unit vectors such that(21)l×s=k,and ϑ is the angle between vector l and the positive direction of the Ox-axis. Then, by virtue of the transformation formulas of first- and second-rank tensor components, the following formulas hold true: (22)ul=ucos⁡ϑ+vsin⁡ϑ,us=-usin⁡ϑ+vcos⁡ϑ,σll=σxxcos2⁡ϑ+σxy+σyxsin⁡ϑcos⁡ϑ+σyysin2⁡ϑ,σss=σxxsin2⁡ϑ-σxy+σyxsin⁡ϑcos⁡ϑ+σyycos2⁡ϑ,σls=σyy-σxxsin⁡ϑcos⁡ϑ+σxycos2⁡ϑ-σyxsin2⁡ϑ,σsl=σyy-σxxsin⁡ϑcos⁡ϑ-σxysin2⁡ϑ+σyxcos2⁡ϑ,μlz=μxzcos⁡ϑ+μyzsin⁡ϑ,μsz=-μxzsin⁡ϑ+μyzcos⁡ϑ,μzl=μzxcos⁡ϑ+μzysin⁡ϑ,μzs=-μzxsin⁡ϑ+μzycos⁡ϑ.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 φj and ψj are taken for each index j=1,2,…,N (N is some natural number) in the following manner: (23)φj,ψj=aj,bjln⁡x2+y2,and solution χj of the Helmholtz equation is represented as(24)χj=cjK0γx2+y2,where K0γx2+y2 is a modified Bessel function of second kind or a Macdonald function (of zeroth order) which tends to infinity at point (0,0) and vanishes at infinity [45] and γ=4μα/(ν+β)(μ+α) and  aj,bj,cj,  j=1,2,…,N, are the real coefficients to be defined.By formula (23) the first-order partial derivatives of the harmonic functions φj and ψj have the form(25)∂xφj,∂xψj=aj,bjxx2+y2,∂yφj,∂yψj=aj,bjyx2+y2,and the partial derivatives of the metaharmonic functions χj are written as (26)∂xχj=-cjγxK1γx2+y2x2+y2,∂yχj=-cjγyK1γx2+y2x2+y2.The second derivatives of these functions appear in formulas (20):(27)∂xxφj,∂xxψj=-∂yyφj,∂yyψj=aj,bjy2-x2x2+y22,∂xyφj,∂xyψj=-aj,bj2xyx2+y22,∂xxχj=cjγ22K0γx2+y2-y2-x2x2+y2K2γx2+y2,∂yyχj=cjγ22K0γx2+y2+y2-x2x2+y2K2γx2+y2,∂xyχj=∑j=1Ncjγ2xyx2+y2K2γx2+y2.If we substitute the corresponding expressions (23)–(26) in formulas (19), then for each index j=1,2,…,N we obtain solution (2μuj,2μvj,2μωj) of system (2) satisfying this system everywhere except the origin where it has a singularity:(28)2μujx,y=λ+3μλ+μlnx2+y2-x2x2+y2aj-xyx2+y2bj+yK1γx2+y2x2+y2cj,2μvjx,y=-xyx2+y2aj+λ+3μλ+μlnx2+y2-y2x2+y2bj-xK1γx2+y2x2+y2cj,2μωjx,y=λ+2μλ+μ1x2+y2-yaj+xbj+2μν+βK0γx2+y2cj.The values of stresses and moment stresses σxx,σyy,…,μyz are obtained if representations (25)–(27) are substituted in the respective formulas (20):(29)σxxjx,y=λ+2μλ+μxx2+y2-xy2-x2x2+y22aj-λλ+μyx2+y2-yy2-x2x2+y22bj-γ2xyx2+y2K2γx2+y2cj,σyyjx,y=λλ+μxx2+y2+xy2-x2x2+y22aj+λ+2μλ+μyx2+y2+yy2-x2x2+y22bj+γ2xyx2+y2K2γx2+y2cj,⋮μzyjx,y=ν-βλ+2μ2μλ+μy2-x2x2+y22aj-2xyx2+y22bj-ν-βν+βγyK1γx2+y2x2+y2cj.After that, eachjth function of (28) and (29) is shifted by value (ξj,ηj) (Figure 1). For this, in formulas (28) and (29) and also in the expressions for stresses and moment stresses variables x and y are replaced, respectively, by values x-ξj and y-ηj. Then functions uj(x-ξj,y-ηj),vj(x-ξj,y-ηj),ωj(x-ξj,y-ηj), as well as also functions σxxj(x-ξj,y-ηj),σyyj(x-ξj,y-ηj),…,μzyj(x-ξj,y-ηj), have singularities at points (ξj,ηj).Figure 1
Shift of a singular point.Let us consider the following sums:(30)u-=∑j=1Nujx-ξj,y-ηj,v-=∑j=1Nvjx-ξj,y-ηj,ω-=∑j=1Nωjx-ξj,y-ηj,σ-xx=∑j=1Nσxxjx-ξj,y-ηj,σ-yy=∑j=1Nσyyjx-ξj,y-ηj,⋮μ-zy=∑j=1Nμzyjx-ξj,y-ηj.For our algorithm approximate solutions of boundary value problems are sought in form (30). Functions (30) satisfy system of (2) everywhere except points (ξj,ηj), j=1,2,…,N, where they have a singularity. The unknown coefficients aj,bj,cj,  j=1,2,…,N, 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 L having almost everywhere the external normal l. In domain Ω, we seek such a solution of system (2) that satisfies the following boundary conditions on boundary L: (31)ulL=f1x,y,usL=f2x,y,ωL=gx,y,x,y∈L,where f1,f2, and g are the functions given on the boundary.In the first place, on the boundary we more or less uniformly distributeN points with coordinates (x1,y1),(x2,y2),…,(xN,yN). As to points (ξj,ηj), 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 (ξj,ηj) lying on some external contour and points (xi,yi) lying on the boundary contours.An approximate solution is sought in the form(32)u-l=∑j=1Nujx-ξj,y-ηjcos⁡ϑ+vjx-ξj,y-ηjsin⁡ϑ,u-s=∑j=1N-ujx-ξj,y-ηjsin⁡ϑ+vjx-ξj,y-ηjcos⁡ϑ,ω-=∑j=1Nωjx-ξj,y-ηj.We substitute the coordinates of points (xi,yi) marked on the boundary and the values of angles ϑi 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(33)∑j=1Nujxi-ξj,yi-ηjcos⁡ϑi+vjxi-ξj,yi-ηjsin⁡ϑi=f1xi,yi,∑j=1N-ujxi-ξj,yi-ηjsin⁡ϑi+vjxi-ξj,yi-ηjcos⁡ϑi=f2xi,yi,∑j=1Nωjxi-ξj,yi-ηj=gxi,yi,i=1,2,…,N.But according to formulas (28), the latter system is a system of 3N linear algebraic equations with 3N unknowns aj,bj,cj,  j=1,2,…,N:(34)∑j=1NA1ijaj+B1ijbj+C1ijcj=f1xi,yi,∑j=1NA2ijaj+B2ijbj+C2ijcj=f2xi,yi,∑j=1NA3ijaj+B3ijbj+C3ijcj=gxi,yi,where i=1,2,…,N and  Akij,Bkij,Ckij,  k=1,2,3, are the already known coefficients. For example, A1ij has the form(35)A1ij=12μλ+3μλ+μln⁡xi-ξj2+yi-ηj2-xi-ξj2xi-ξj2+yi-ηj2cos⁡ϑi-12μxi-ξjyi-ηjxi-ξj2+yi-ηj2sin⁡ϑi.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 L1,L2, of which L2 encircles L1, while L1 lies around the origin.The harmonic functionsφ and, ψ in the general solutions (19) and (20) are represented as follows:(36)φ=aln⁡r+a0+∑n=1Nr-na1ncos⁡nϑ+a2nsin⁡nϑ+rna3ncos⁡nϑ+a4nsin⁡nϑ,ψ=bln⁡r+b0+∑n=1Nr-nb1ncos⁡nϑ+b2nsin⁡nϑ+rnb3ncos⁡nϑ+b4nsin⁡nϑ,where r(x,y)=x2+y2 and(37)ϑx,y=arctan⁡yx,x>0,arctan⁡yx+π,x<0,y≥0,arctan⁡yx-π,x<0,y<0,π2,x=0,y>0,-π2,x=0,y<0.Function χ is written in the form(38)χ=cI0γr+c0K0γr+∑n=1NInγrc1ncosnϑ+c2nsinnϑ+Knγrc3ncos⁡nϑ+c4nsin⁡nϑ,where In(γr) and Kn(γr) are, respectively, a modified Bessel and a Macdonald function of order n [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 a,a0,b,b0,c,c0,aαn,bαn,cαn,  n=1,2,…,N,  α=1,2,3,4.Further, we more or less uniformly distribute2(2N+1) points with coordinates (x1,y1),(x2,y2),…,(x2(2N+1),y2(2N+1)) on contours L1 and L2 (Figure 3). In order that the given boundary conditions be satisfied at these points (Figure 3), we substitute the coordinates of the boundary points (xi,yi),  i=1,2,…,2(2N+1), 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 a,a0,b,b0,c,c0,aαn,bαn,cαn,  n=1,2,…,N,  α=1,2,3,4, we obtain a system of 6(2N+1) linear algebraic equations with 6(2N+1) 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 Section5, 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:λ=82.098765432 GPa, μ=35.185185185 GPa, α=7.037037036 GPa, and ν+β=0.7037037036 GN. These data correspond to the elastic characteristics of brass.Example 1.
Let domainΩ1 be a square with a circular hole Ω1=Ω10∖Ω-0, where Ω10=(x,y)∣-3<x<3,-3<y<3 and Ω0=(x,y)∣x2+y2<4 (Figure 4).
Consider the so-called Kirsch problem [19, 45] for domain Ω1 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:(39)σxx=1.0,σxy=0,μxz=0,x=±3,-3<y<3,σyy=0,σyx=0,μyz=0,y=±3,-3<x<3,σll=0,σls=0,μlz=0,x2+y2=4.This problem is solved approximately by the proposed method both for the Cosserat medium and for the classical elastic medium.
We arrange points(ξj,ηj),  j=1,2,…,72, on the square containing domain Ω1 (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 72×3=216 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(40)aj=-a25-j,bj=b25-j,cj=-c25-j,j=13,14,…,24,aj=-aj-24,bj=-bj-24,cj=cj-24,j=25,26,…,36,aj=a49-j,bj=-b49-j,cj=-c49-j,j=36,37,…,48,aj=-a110-j,bj=b110-j,cj=-c110-j,j=56,57,…,61,aj=-aj-12,bj=-bj-12,cj=cj-12,j=62,63,…,67,aj=a122-j,bj=-b122-j,cj=-c122-j,j=68,69,…,72.After solving the resulting system we compared the boundary values ​​of the obtained solutions with the given boundary conditions. If absolute errors of order 10-2 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 10-6.Figure 4
The considered domainΩ1.Figure 5
Unilateral tension of domainΩ1.Figure 6
DomainΩ1 encircled by singular points.In the case of the Kirsch problem it is of special interest to investigate the distribution of stressσss along the hole contour. A maximal value of stress σss both for the classical elastic medium and for the Cosserat medium is obtained, as expected, at ϑ=±π/2. For the classical medium, σssmax=13.7 and this value does not depend on the material (if the material is isotropic and homogenous). For the Cosserat medium, σssmax depends on the material and in our case its value is σssmax=12.5 (Figure 7).Figure 7
Graphs of the distribution of stressσss along the hole contour for the classical elastic medium and the Cosserat medium for domain Ω1.Example 2.
Now let us consider the Kirsch problem for domainΩ2=Ω20∖Ω-0,  Ω20=(x,y)∣-4<x<4,-4<y<4 (Figure 8).
As expected, in this case the values ofσssmax 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 σssmax 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 a/r, where a is the side of the square and r is the hole radius, increases to a certain value, values σssmax 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Ω2.Figure 9
Graphs of the distribution of stressσss along the hole contour for the classical medium and the Cosserat medium for domain Ω2.Example 3.
We consider domainV=-2.5<x<2.5,-2<y<2 (Figure 10). In domain V, it is required to find a solution of system (2) (where λ=2,  μ=1,  α=0, and ν+β=0) that satisfies the following conditions: (41)u=-0.38y2±1.5y+0.125,x=±2.5,-2≤y≤2,v=-0.14y2±0.5y-2.125,x=±2.5,-2≤y≤2,σyyy=2-σyyy=1=0.08,-2.5<x<2.5,σyxy=2=-0.08x-1.12,-2.5<x<2.5,u=0.02x2-1.2x-1.52,y=-2,-2.5<x<2.5,v=-0.34x2-0.4x-0.56,y=-2,-2.5<x<2.5.
The exact solution of this problem is as follows:(42)u=0.02x2-0.38y2+0.6xy,v=-0.34x2-0.14y2+0.2xy.
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(43)u=0.5∑n=136rn-156rcos⁡nϑ-n2ysin⁡n-1ϑ-n2xcos⁡n-1ϑan+56rsin⁡nϑ+n2ycos⁡n-1ϑ-n2xsin⁡n-1ϑbn-nsin⁡n-1ϑcn+ncos⁡n-1ϑdn,v=0.5∑n=136rn-156rsin⁡nϑ+n2xsin⁡n-1ϑ-n2xcos⁡n-1ϑan-56rcos⁡nϑ+n2ycos⁡n-1ϑ+n2xsin⁡n-1ϑbn-ncos⁡n-1ϑcn-nsin⁡n-1ϑdn.The stress tensor components σyy and σyx figuring in the initial condition are represented as finite series:(44)σyy=∑n=136nrn-2rcos⁡n-1ϑ+n-12ysin⁡n-2ϑ+n-12xcos⁡n-2ϑan+rsin⁡n-1ϑ-n-12ycos⁡n-2ϑ+n-12xsin⁡n-2ϑbn+n-1sin⁡n-2ϑcn-n-1cos⁡n-2ϑdn,σyx=∑n=136nn-12rn-2-ycos⁡n-2ϑ+xsin⁡n-2ϑan-ysin⁡n-2ϑ+xcos⁡n-2ϑbn-2cos⁡n-2ϑcn-2sin⁡n-2ϑdn.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 (a1,…,a36,  b1,…,b36,  c1,…,c36,  d1,…,d36). After solving this system by the above formulas we can easily define the displacement vector and stress tensor components.Figure 10
DomainV in which the nonlocal problem of plane elasticity is solved.In our study we used Maple 17 software. The numerical results are presented in Table1, where u- and v- denote the approximate values of the displacement vector components.Table 1
Numerical results for Example3.
(
x
,
y
)
u
-
(
x
,
y
) v-(x,y)
u
(
x
,
y
) v(x,y)
u
-
-
u v--v
(
0,0
)
-
7.205904960
⋅
1
0
-
9 -2.170915985⋅10-7
0 0
7.2059
⋅
1
0
-
9 2.171⋅10-7
(
1.0
,
1.0
)
0.2399999816 -0.2800003047
0.24 -0.28
1.84
⋅
1
0
-
8 3.047⋅10-7
(
-
1.5,1.5
)
-
2.159999943 -1.530000297
-
2.16 -1.53
5.7
⋅
1
0
-
8 2.97⋅10-7
(
1.2,1.3
)
0.3225999578 -0.4142003246
0.3226 -0.4142
4.22
⋅
1
0
-
8 3.346⋅10-7
(
-
1.8
,
-
2.3
)
0.5386005059 -1.014199600
0.5386 -1.0142
5.059
⋅
1
0
-
7 4⋅10-7
(
2.2
,
-
1.4
)
-
2.496000007 -2.535999981
-
2.4960 -2.5360
7
⋅
1
0
-
9 1.9⋅10-8
(
1.25
,
2.15
)
-
0.1128001215 -0.6409003984
-
0.1128 -0.64090
1.215
⋅
1
0
-
7 3.984⋅10-7From 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.



---

*Source: 3845362-2016-02-29.xml*