Microsoft Word - numero_49_art_23_2487 V. Matveenko et alii, Frattura ed Integrità Strutturale, 49 (2019) 225-242; DOI: 10.3221/IGF-ESIS.49.23 225 Focused on Russian mechanics contributions for Structural Integrity Construction of analytical eigensolutions for isotropic conical bodies and their application for estimation of stresses singularity V. Matveenko, I. Shardakov Institute of Continuous Media Mechanics of the Ural Branch of Russian Academy of Sciences, Russia mvp@icmm.ru shardakov@icmm.ru T. Korepanova Institute of Continuous Media Mechanics of the Ural Branch of Russian Academy of Sciences, Russia Perm Military Institute of the National Guard Troops of the Russian Federation, Russia ton@icmm.ru ABSTRACT. A complete set of eigensolutions is constructed for different variants of circular conical bodies: homogeneous cone with one lateral surface (solid cone), homogeneous cone with two lateral surfaces (hollow cone) and composite cone for different boundary conditions on the lateral surface. It has been shown that the constructed eigensolutions can be readily applied for estimation of the character of stress singularity at the vertices of conical bodies. The character of stress singularity at the vertex of the solid and hollow cones for different boundary conditions on the lateral surfaces has been defined by direct numerical simulations. Numerical results obtained for solid, hollow and compose cones under different boundary conditions on the lateral surfaces are discussed. KEYWORDS. Stress singularity; Conical bodies; Eigensolutions. Citation: Matveenko, V., Shardakov, I., Korepanova, T., Construction of analytical eigensolutions for isotropic conical bodies and their application for estimation of stresses singularity , Frattura ed Integrità Strutturale, 49 (2019) 225-242. Received: 17.04.2019 Accepted: 24.05.2019 Published: 01.07.2019 Copyright: © 2019 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. INTRODUCTION ne of the main findings of the classical elasticity theory is the possibility of existence of singular solutions due to the occurrence of infinite stresses at the points of surface non-smoothness, changes in the type of boundary conditions, contact of different materials and inside the body at the points of contact of dissimilar materials. An example of theoretical justification of this statement can be found in Kondratiev (1967), where it was shown that the solutions to the equations of the linear theory of elasticity in the vicinity of angular points can be expressed as O http://www.gruppofrattura.it/VA/49/2487.mp4 V. Matveenko et alii, Frattura ed Integrità Strutturale, 49 (2019) 225-242; DOI: 10.3221/IGF-ESIS.49.23 226 nα -1 1 2 nn n 1 τ r 0 < ...c < Reα < Reα < ... < Reαf r~ , ,    (1) or can be written in a more complicated form involving logarithmic components in the case of multiple spectrum points n . Here r is the distance to the angular point, nf is the function of an angular distribution of the stress field  in the vicinity of an angular point, which in 2D case depends on one polar angular variable  , 0c  , and in 3D case - on two spherical coordinates , , 1 / 2с    . Such a representation of the problem solution suggests that if there are n with Re 1n  , the stresses tend to infinity at r tends to zero. One of the trends in the construction of the solutions of form (1) concerns the analysis of regions with specific configuration. In two-dimensional problems these are plane wedges. For almost half-a century history of studying these problems the researchers have investigated nearly all types of wedge-shaped bodies: homogeneous, composite, isotropic, anisotropic, etc. In [2-3] the authors give an extended review of the literature concerning the construction and analysis of singular solutions in two-dimensional problems of the theory of elasticity. In three-dimensional problems one can differentiate between two types of regions: an edge of a 3D wedge (an edge is not necessarily rectilinear, the wedge angle can vary along the edge) and a vertex of a polyhedral wedge or a cone. The results of some publications have disengaged further interest in the problems of the first class. In these articles (for example, in [4]) it has been shown that the solution of the plane and anti-plane problems for wedges obtained in the planes perpendicular to the edge of a 3D wedge defines the type of stress singularity at the edge points, through which the plane passes. In recent years, considerable attention has been focused on studying the stress singularity at the vertices of polyhedral wedges or cones. These problems are solved by different algorithms of the finite element method, the boundary elements method or by applying the Mellin transformation to the two-dimensional boundary integral equations. Out of numerous papers, which have used the concepts of these numerical methods, one should set aside papers [5-12]. As in the other branches of the elasticity theory, the analytical methods are of considerable importance both for numerical simulation and testing of numerical methods. In three-dimensional problems the analytical methods are generally applied to circular cones. One of the pioneering works in this field Bazant and Keer [13] is concerned with the stress singularity problem for a solid cone under axisymmetric deformation and torsion, for which the boundary conditions are specified in terms of displacement and stresses. A detailed theoretical analysis of the construction of axisymmetric solutions for elastic cone is given in [14]. The method of construction of singular and regular solutions to the Laplace and Navier-Stokes equations written for the axisymmetric domain is described in [15]. The idea of the method is the expansion of the desired solution in the Fourier series in the cyclic coordinate for axisymmetric coordinate systems, which is followed by the numerical solution of a sequence of two-dimensional problems for series coefficients. In the papers that followed the above cited studies, the analytical solutions were constructed for some particular problems. For example, in [16-17] the authors presented the results for a composite cone under axisymmetric deformation. In this case, a composite cone is a body composed of two embedded cones, which have a common contact surface. The solutions were obtained for the perfect bonding and frictionless slip conditions. The axisymmetric problem for a circular cone made of transversally isotropic material was considered in [18]. A rather complete list of papers dealing with the investigation of circular cones by the analytical methods can be found in review article [3]. The most comprehensive research on this topic was carried out in [19], in which the authors constructed the analytical solution for a solid circular cone and presented numerical data on the character of stress singularity at the vertex of a circular solid cone with the boundary conditions at the lateral surface specified in terms of the displacements and stresses. Compared to the known results, the present work presents a full spectrum of analytical eigensolutions, which allow constructing the solutions not only for a homogeneous isotropic cone but also for other types of isotropic cone, as for example, for hollow and composite cones for different boundary conditions on the lateral surfaces. MATHEMATICAL STATEMENT OF THE PROBLEM OF CONSTRUCTING EIGENSOLUTIONS FOR SEMI- INFINITE CIRCULARCONICAL BODIES onsider a homogeneous circular cone (Fig. 1a), the vertex of which coincides with the center of spherical coordinates , ,r   , and the base is perpendicular to the axis 0  . The volume occupied by the cone is defined as 1 00 , , 0 2r            , and its boundary is defined by the coordinate surfaces 1 0,     . C V. Matveenko et alii, Frattura ed Integrità Strutturale, 49 (2019) 225-242; DOI: 10.3221/IGF-ESIS.49.23 227 The case when 1 0  corresponds to a solid cone. Figure 1: Hollow cone (a); composite cone (b). We state the problem of constructing eigensolutions, which will satisfy the homogeneous equations of equilibrium (1 ) 0S grad div rot rot  u u (2) (here 1/1 2S n= - , n is Poisson's ratio, u is the displacement vector) and homogeneous boundary conditions at the surfaces 1 0,q q q q= = for displacements 0, 0, 0ru u u    (3) and stresses 0, 0, 0r       (4) or mixed boundary conditions, which in the context of solid mechanics correspond to a perfect- slip boundary condition at the lateral surface 0, 0, 0ru      (5) For the examined body of revolution and boundary conditions (3)-(5), the eigensolutions can be represented as a Fourier series [20] in the circumferential coordinate                          0 0 0 1 1 1 , , sin , , v v sin , , cos r w k k k k k k u r u r u r k u r r r k u r w r r k                                                (6) Here the dependence on the radius is expressed according to (1), , ,ru u u  -are the components of the displacement vector along the , ,r   -axes, , ,r     are the components of the stress tensor,  is the characteristic index. V. Matveenko et alii, Frattura ed Integrità Strutturale, 49 (2019) 225-242; DOI: 10.3221/IGF-ESIS.49.23 228 If 1 0  , the region under consideration is bounded by only one coordinate surface 0  , and at 0  is assumed to meet the regularity conditions / 0, 0, 0ru u u     (7) In the framework of the proposed problem formulation we can consider a composite cone occupying a region (1) ( 2 )V V V  , where a subregion (1)V (subregion ( 2 )V ) is made of the material with the shear modulus ( 1) ( ( 2 ) ) and Poisson’s ratio (1) ( ( 2 ) ) and its geometry is defined by the relations 0 r   , 0 2   , 2 0    ( 1 2    ) (Fig. 1b). In particular cases 1 and 0 can be equal to 0 and  , respectively. For a composite cone eigensolutions (6) are constructed for each subregion. At the contact line 2  we can prescribe perfect bonding conditions ( 1) ( 2 ) ( 1) ( 2 ) ( 1) ( 2 ), , ;r ru u u u u u      ( 1) ( 2 ) ( 1) ( 2 ) ( 1) ( 2 ), ,r r             (8) or perfect slip conditions ( 1) ( 2 ) ( 1) ( 2 ) ( 1) ( 2 ) ( 1) ( 2 ); , 0r ru u                  (9) Upon substituting Eqns. (6) into the equilibrium Eqn. (2) and going to a new independent variable (1 cos )/ 2x   [20] we obtain the following equations for each harmonic of the Fourier series:                            22 1 2 2 2 4 1 1 1 2 4 1 1 v 1 v 0 2 21 1 k k k k k k xH x kd u x du x x x x u x dx x xdx x x H d x kw xH x x dxx x x x                         (10)                               22 2 1 1 12 1 1 3 4 1 1 1 2 4 1 1 1 2 1 1 0 2 4 1 k k k k k k xG x k Gd v x dv x G x x G x v x dx x xdx k G G k xdw xd G x x u x w x dx dx x x                       (11) ( )( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 22 1 2 1 13 4 1 1 1 1 2 4 1 1 1 2 1 0 2 4 12 1 k k k k k k xG x G kd w x dw x x x x w x dx x xdx G k dv x G k xkG u x v x dx x xx x 2 é ù- + +ê úë û- + - + + - é ù- + -ê ú + + + ⋅ =ê ú ê ú-- ë û (12) Here we introduce the following notation:        1 22 1 1 2 / 2 1 , 3 4 /( 2 1), H H                        1 2 32 1 / 1 2 , 1 , 2 4 4 / 1 2 .G G G               Using (6) we can reduce boundary conditions (3)-(5) and regularity conditions (7) to condition for V. Matveenko et alii, Frattura ed Integrità Strutturale, 49 (2019) 225-242; DOI: 10.3221/IGF-ESIS.49.23 229 displacements ( ) ( ) ( )0 , 0 , 0 , k k ku x v x w x= = = (13) stresses        1 1 0k k du x x x v x dx                                  2 1 1 1 11 0 2 1 2 1 k k k k dv x S S u x S x x dx S k S x v x w x x x x x                         (14)               1 2 1 0 2 1 2 1 k k k dw x x k x x w x v x dx x x x x                    1 0 ; 0 ; 0.k k k du x x x v x w x dx     (15) CONSTRUCTION OF PARTIAL SOLUTION FOR ZEROTH HARMONIC OF THE FOURIER SERIES solution for zeroth harmonic of the Fourier series is considered separately, because it does not follow explicitly from the algorithm for construction of partial solutions to the system of differential Eqns. (10)-(12) for any 0k  . At 0k  there are two problems: axisymmetric rotation and axisymmetric deformation. In the first problem, the component of the displacement vector 0w is defined by Eqn. (12). In the problem of axisymmetric deformation the components of the displacement vector 0 0,u v are defined by Eqns. (10), (11). The solutions for the function 0w are constructed in the form of the generalized power series    0 0 m m m w x R x         (16) where mR are the coefficients of the series,  is the characteristic index. The possibility of constructing a solution in the form of (16) was proved in [21]. For Eqn. (12) the point 0x  is a regular singular point. In this case, one of the partial solutions is written as a power series (16), for which the range of convergence lies in the interval 0 1x  , because the value 1x  is the zero of the high-order derivative. function closest to the point 0x  . The coefficients of the series mR and the characteristic index  are found by substituting (16) into (12). Setting to zero the expressions with like powers of x , yields a recurrence relation for mR :            1 2 2 2 1 2 2 1 4 1 2 2 1 1 4 2 1 0,( =0, if 0 ) m m m m m m R m m R m m R R m                                    (17) From the condition of existence of non-zero solution for 0R we obtain the characteristic equation: ( 2 1)( 2 1) 0    (18) for which 1 1/ 2  and 2 1/ 2   are the roots. A V. Matveenko et alii, Frattura ed Integrità Strutturale, 49 (2019) 225-242; DOI: 10.3221/IGF-ESIS.49.23 230 According to the theory of differential equations [21], there always exists a solution in the form of a generalized power series (16), that fits the largest root 1 . Substituting the value of the root 1 into (17), we obtain a recurrence relation for ( 1) mR            2 2 (1) ( 1) (1) (1) 1 2 0 2 2 1 2 2 3 2 , 0, 1 1 4 1 m m m m m m m R R R m R m m m m                    (19) Here in after, the superscript will denote the number of the partial solution. After completing these transformations we obtain the first partial solution for Eqn. (12) represented as a generalized power series   1 (1) (1) 2 0 0 m m m w x R x                (20) A quantity that plays an important part in the construction of the second, linearly independent partial solution in the form of the generalized power series is the difference in the roots of the characteristic equation [21], i.e. the number 1 2 1     . If the number g is not a positive integer there exists a second, linearly independent solution in the form of the generalized power series (16). However, if  is a positive integer, the existence of the second partial solution in the general case in the form of a generalized power series (16) is not assured. This is exactly the case we have, because 1 2 1.     To eliminate this ambiguity, we use an approach, which is based on a subsequent reduction of the initial differential equation and retention of a certain number of series terms. The performance of this technique is illustrated by way of constructing the second particular solution  ( 2 )0w x to Eqn. (12) at 0k  , which corresponds to axisymmetric deformation. Eqn. (12) at 0k  takes the following form:                22 10 0 02 4 1 1 1 1 2 0 4 1 xG x G kd w x dw x x x x w x dx x xdx           (21) The construction of the second particular solution proceeds as follows. At the first stage using the obtained first particular solution  (1)0w x (20) we reduce the second order differential Eqns. (21) to the first order equation replacing the desired solution in Eqn. (21) by the relation [21]    (1)0 0 ( )w x w x F x dx  (22) where ( )F x is a new unknown function. Upon this substitution we obtain the linear first-order differential equation for 1 0 ( ) ( ) ( ) ( ) 0 dF x P x P x F x dx   (23) where (1) 1 0 (1) (1) 2 0 0 0 ( ) (1 ) ( ) ( ) ( ) (1 ) ( ) 2( ) . P x x x w x dw x P x x w x x x dx       (24) V. Matveenko et alii, Frattura ed Integrità Strutturale, 49 (2019) 225-242; DOI: 10.3221/IGF-ESIS.49.23 231 For this equation we can also construct a solution [21] for F(x) in the form of a generalized power series ( ) 0 ( ) mm m F x N x         (25) In the traditional approach, the coefficients mN and index  of this series are determined by substituting (25) into Eqn. (23). By setting the expressions of the same power of x equal to zero we can obtain the recurrence relations for mN and from the condition of the existence of zero solution for 0N we get the characteristic equation, determining the index  . However, it should be noted that practical implementation of this procedure presents considerable difficulties, because the expressions in front of the function ( )F x and its derivatives in Eqn. (23) will appear in the form of in finite series. Moreover, the fact that that after construction of series (25) it will be necessary to proceed to the solution ( )0w x through the relation (22), implies that analytical transformations of infinite series will become increasingly complex. In this regard, we suggest using an algorithm, which allows us to circumvent the above difficulties in the construction of a particular solution. To illustrate our ides we will apply the proposed algorithm to Eqn. (23). 1. At the first stage we define the form of the characteristic equation for the index . For this purpose, we retain in the series (25) only the first term with the coefficient 0N , and in the expressions (24) in front of the function ( )F x and its derivatives in Eqn. (23) we retain only the lowest powers of the variable x , which, however are not less than the order of the highest derivative (in this case, for Eqn. (23) this is the first order). After these manipulations the coefficients (24) take the following form: (1) 1/ 2 1 0 (1) 1/ 2 (1) 3/ 2 2 (1) 1/ 2 (1) 1/ 2 0 0 1 0 1 ( ) (1 ) ( ) (1 2 )( ) 2( )(1/ 2 3/ 2 ). P x x x R x P x x R x R x x x R x R x         (26) 2. Upon substituting relations (25) and (26), in which we retain only the first term with the coefficient 0N , into Eqn. (23) and setting to zero the expression with the lowest power of x , we obtain 0 ( 2) 0N    (27) From (27) it follows that the characteristic index 2   . With a knowledge of this index we can formally represent the structure of the solution for ( )F x : 2 1 0 0 1 2 3( ) ...F x N x N x N x N x       (28) Using the obtained solution (28), we can find a formal representation of the second particular solution for  0w x based on the relation (22). In this case the number of the series terms retained in the series (28) and (20) should be no less than  (namely, no less than 2). Substitution of these segments of the series into relation (22) yields   3 (1) ( 2 ) ( 2 ) (1/ 2 ) 0 0 ( ) ( ) ln( )m m m w x F x dx R x T x x        (29) Thus, after generalization of this form of relation to the infinite number of terms we obtain the following form of the second particular solution:       1/ 2( 2 ) ( 2 ) ( 2 )0 0 lnm m m x R T x x         (30) In this expression the coefficients ( 2 )mR , ( 2 ) mT are still to be defined. V. Matveenko et alii, Frattura ed Integrità Strutturale, 49 (2019) 225-242; DOI: 10.3221/IGF-ESIS.49.23 232 3. Substituting relation (30) into the initial differential Eqn. (21) and equating to zero the expressions with the same powers of x , we obtain the recurrence relations for determination of the coefficients ( 2 )mR , ( 2 ) mT           2 ( 2 ) ( 2 ) ( 2 ) 1 2 1 2 3 2 3 2 2 5 2 1 4 1 m m m m m m m T T T m m m m                                          2 ( 2 ) 2 ( 2 ) ( 2 ) 1 1 ( 2 ) ( 2 ) 2 2 1 2 31 2 4 5 1 1 1 2 3 2 2 5 2 2 2 4 1 1 m m m m m m m mm m R T R T m m m m m m m m m R T m m m m                               (31)  ( 2 ) ( 2 ) ( 2 ) ( 2 )0 1 0 1( 1, 1, 1 , 1, 1)m T T R R        The proposed algorithm is valid for the linear ordinary differential equations of any order. Hence, the proposed method allows us to define step by step the types of the generalized power series for all partial solutions of the initial differential equation and in all partial solutions to set apart the regular from irregular solutions (in our case for 0x  ). Such a capability of the method is rather essential for constructing solutions to particular problems, for example, for hollow and composite cones. The form of the obtained solutions (1) ( 2 )0 0 , w w suggests that (1) 0w is a regular solution, and ( 2 ) 0w is irregular solution at 0x  . A general solution to differential Eqn. (12) takes the following form:      (1) ( 2 )0 1 0 2 0w x C w x C w x    (32) where 1C , 2C are the constants defined by the prescribed combination of the boundary conditions (3) - (5). To construct partial solutions to Eqns. (10), (11), corresponding to the axisymmetric deformation, we need to solve this system for 0v                      3 22 0 0 0 0 3 22 1 1 4 2 2 1 1 2 x x S d u x d u x du x v x x x x S S dxdx dx                    (33) and with respect to the function 0u , which yields a differential equation of the fourth order ( ) ( ) ( )( ) ( ) ( )( )( )[ ] ( ) ( )( ) ( ) ( )( )( ) ( ) 4 3 2 22 0 0 0 4 3 2 20 0 1 1 4 8 2 2 3 2 1 2 2 2 1 1 2 0 d u x d u x d u x x x x x x x x dx dx dx du x x u x dx a a a a a a a a - - - - + - + - - + - - + - - + - + = (34) This equation is the differential equation with a regular singular point. Therefore, the linearly independent partial solutions can be represented in the form of convergent generalized power series. Using the above procedure for constructing such series, we obtain four partial solutions (1) ( 2 ) (3) ( 4 )0 0 0 0, , , u u u u written as                   1(1) (1) ( 2) ( 2 ) 0 0 0 0 1( 3) ( 3) ( 3) ( 4 ) ( 4 ) ( 4 ) 0 0 0 0 , , ln , ln , m m m m m m m m m m m m m m u x A x u x A x u x A B x x u x A B x x                                 (35) V. Matveenko et alii, Frattura ed Integrità Strutturale, 49 (2019) 225-242; DOI: 10.3221/IGF-ESIS.49.23 233 where the coefficients (1) ( 2 ) ( 3 ) ( 4 ) ( 3 ) ( 4 ) , , , , , m m m m m mA A A A B B are determined from the recurrence relations available on E- resource (www. icmm.ru/compcoeff) and can be calculated using the suggested options. For this purposes, it is necessary to input the number of terms of a series m , Poisson's ratio  , the number of harmonics k and the desired  , which can be either complex or real. Substituting (35) into (33), we obtain partial solutions (1) ( 2 ) ( 3 ) ( 4 )0 0 0 0, , , v v v v for the function 0v ( ) ( ) ( )[ ]( ) ( ) ( ) ( )[ ]( ) ( ) ( ) ( )[ ]( ) ( ) ( )( ) ( ) ( ) ( )[ ]( ) ( )( ) ( ) (1) (1) 0 2 0 (2) (2) 0 2 0 (3) (3) (3) 0 2 0 (4) 1(4) (4) 0 2 0 1 , 1 2 1 , 1 2 1 1 ln , 1 2 1 ln 1 2 m m m m m m m m m m m m m m x x v x P x S x x v x P x S x x S v x P D x x xS x x v x P D x x S a a a a a a a a a a a a ¥ = ¥ = ¥ = ¥ - = - é ù= ê úë û- - + - é ù= ê úë û- - + ì üï ï- + é ùï ï= + +í ýê úï ïë û- - + ï ïî þ ì- é ù= + ⋅ê úë û- - + å å å å , üï ïï ïí ý ï ïï ïî þ (36) where the coefficients (1) (2) (3) (4) (3) (4) , , , , , m m m m m mP P P P D D are available on E-resource (www. icmm.ru/compcoeff) and can be calculated ibidem. A general solution for 0u and 0v can be written as                     (1) ( 2 ) ( 3) ( 4 ) 0 1 0 2 0 3 0 4 0 (1) ( 2 ) ( 3) ( 4 ) 0 1 0 2 0 3 0 4 0 ,u x C u x C u x C u x C u x v x C v x C v x C v x C v x                 (37) where 1 2 3 4, , ,C C C C are the constants defined by the prescribed combination of boundary conditions (3) - (5). CONSTRUCTION OF PARTIAL SOLUTION FOR NONZEROTH HARMONICS OF THE FOURIER SERIES he construction of partial solutions to Eqns. (10)-(12) involves some transformations [21], which yield a system of two differential equations with respect to , k kw v                                             4 3 2 4 3 2 1 04 3 2 2 3 2 2 0 3 2 1 02 3 2 0k k k k k k k k k k k d w x d w x d w x dw x f x f x f x f x f x w x dxdx dx dx d v x d w x d w x dw x x x v x x x x x w x dxdx dx dx                 (38) where 0 1 2 3 4 0 2 0 1 2 3, , , , , , , , , , f f f f f       are given by             2 220 1 1 1 1 2 3 2 1 1 1 1 2 16 f x x x x x k k k                            2 220 1 1 1 1 2 3 2 1 1 1 1 2 16 f x x x x x k k k                       2 21 1 11 2 1 4 3 1 2 2 f x x x x x x k            T V. Matveenko et alii, Frattura ed Integrità Strutturale, 49 (2019) 225-242; DOI: 10.3221/IGF-ESIS.49.23 234       22 2 22 1 1 4 18 1 13 2 f x x x x x k                    3 43 43 46 1 2 1 1f x x x x f x x x     (39)                                    22 2 0 2 2 0 2 2 1 2 3 2 3 2 3 1 1 1 1 ; 1 4 2 11 4 1 1 1 2 11 4 4 1 1 2 1 1 5 2 1 2 x x x k x x x x x x x x k k x x x x x k k x x x x x x x k k                                           Moreover, the results of these transformations is a set of equations, which relate the function ku to the functions , k kw v and their derivatives ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 2 2 2 1 4 1 4 2 1 2 1 2 2 k k k x x d w x dw x u x x x x x x x k S S dxdxa ì- ïï= - + - -í ï- + + ïî ( )( ) ( ) ( ) ( ) ( ) ( )( ) ( )24 1 1 1 1 2 1 2 1 2kk k dv x x x k S w x kSx x k x S v x dx a a üé ùïïé ù- - + + + + - - + - +ê úýê úë û ïê úë ûïþ (40) Eqn. (38) is independent of (39) and represents a linear differential equation of the fourth order with respect to the function kw . Eqn. (39) can be considered as a differential equation of the second order with respect to the function kv with right –hand side depending on kw . Such a peculiarity of the differential Eqns. (28) and the obtained relation (30) allows us to determine the sequence of constructing partial solutions for functions , , k k kw v u . The hang of this sequence can be outlines as follows. First, from the solution of Eqn. (38) we obtain four partiсular solutions ( 1) ( 2 ) ( 3 ) ( 4 ), , , k k k kw w w w , written as               1 1 (1) (1) ( 2 ) ( 2 )2 2 0 0 1 1 ( 3 ) ( 3 ) ( 3) ( 4 ) ( 4 ) ( 4 )2 2 0 0 , , ln , ln k k m m k m k m m m k k m m k m m k m m m m w x A x w x A x w x A B x x w x A B x x                                                                       (41) where the coefficients (1) ( 2 ) ( 3) ( 4 ) ( 3) ( 4 ), , , , , m m m m m mA A A A B B can be determined and calculated E-resource (www. icmm.ru/compcoeff). Then substituting the obtained particular solutions (1) ( 2 ) ( 3 ) ( 4 ), , , k k k kw w w w into the right-hand side of Eqn. (28. b) and solving it as an inhomogeneous equation we arrive at four partial solutions (1) ( 2 ) ( 3 ) ( 4 ), , ,k k k kv v v v , which can be written as ( ) ( ) 1 1 (1) (2)(1) (2)2 2 0 0 , , k k m m m mk k m m v x P x v x P x æ ö æ ö+ -÷ ÷¥ ¥ç ç+ +÷ ÷ç ç÷ ÷÷ ÷ç çè ø è ø = = é ù é ù ê ú ê ú= =ê ú ê ú ê ú ê ú ë û ë û å å ( ) ( ) 1 (3) (3) (3) 2 0 ln , k m m mk m v x P D x x æ ö- ÷¥ ç - ÷ç ÷÷çè ø = ì üï ïï ïï ïé ù= +í ýê úë ûï ïï ïï ïî þ å (42) V. Matveenko et alii, Frattura ed Integrità Strutturale, 49 (2019) 225-242; DOI: 10.3221/IGF-ESIS.49.23 235 ( ) ( ) 1 (4) (4) (4) 2 0 ln , k m m mk m v x P D x x æ ö+ ÷¥ ç - ÷ç ÷÷çè ø = ì üï ïï ïï ïé ù= +í ýê úë ûï ïï ïï ïî þ å where the coefficients (1) ( 2 ) ( 3) ( 4 ) ( 3) ( 4 ), , , , , m m m m m mP P P P D D can be determined and calculated E-resource (www. icmm.ru/compcoeff). Then, solving Eqns. (39) as a homogeneous one, we find two more partial solutions ( 5 ) ( 6 ), k kv v . The form of this differential equation indicates that the point 0x  is a regular singular point. The construction of partial solutions in the form of the generalized power series is accomplished in the framework of the proposed approach. The obtained partial solutions are expressed as   1 ( 5 ) ( 5 ) 2 0 , k m k m m v x P x                     1 ( 6 ) ( 6 ) ( 6 ) 2 0 ln k m k m m m v x P D x x                 (43) where the coefficients ( 5) ( 6 ) ( 6 ), ,m m mP P D can be determined and calculated E-resource (www. icmm.ru/compcoeff). At the next stage, we use the partial solutions (1) ( 2 ) ( 3 ) ( 4 ), , , k k k kw w w w , (1) ( 2 ) ( 3 ) ( 4 ) ( 5 ) ( 6 ), , , , , k k k k k kv v v v v v and the obtained relation (40) to derive six partial solutions ( 1) ( 2 ) ( 3 ) ( 4 ) ( 5 ) ( 6 ), , , , , k k k k k ku u u u u u , which are expressed as        1 (1) (1) 2 0 2 1 1 2 2 k m k m m x x u E x kx x S                           1 ( 2 ) ( 2 ) 2 0 2 1 1 2 2 k m k m m x x u E x kx x S                              1 ( 3) ( 3) ( 3) 2 0 2 1 ln 1 2 2 k m k m m m x x u E G x x kx x S                     (44)           1 ( 4 ) ( 4 ) ( 4 ) 2 0 2 1 ln 1 2 2 k m k m m m x x u E G x x kx x S                            1 ( 5 ) ( 5 ) 2 0 1 1 2 2 k m k m m x x u E x kx x S                              1 ( 6 ) ( 6 ) ( 6 ) 2 0 1 ln 1 2 2 k m k m m m x x u E G x x kx x S                     where the coefficients (1) ( 2 ) ( 3) ( 4 ) ( 5 ) ( 6 ) ( 3) ( 4 ) ( 6), , , , , , , , m m m m m m m m mE E E E E E G G G for any value of 0m  can be determined and calculated E-resource (www. icmm.ru/compcoeff). The general solutions for , , k k ku v w are given by                                             (1) ( 2 ) ( 3) ( 4 ) ( 5 ) ( 6 ) 1 2 3 4 5 6 (1) ( 2 ) ( 3) ( 4 ) ( 5 ) ( 6 ) 1 2 3 4 5 6 (1) ( 2 ) ( 3) ( 4 ) ( 5 ) ( 6 ) 1 2 3 4 5 6 (1 1 k k k k k k k k k k k k k k k k k k k k k k k u x C u x C u x C u x C u x C u x C u x v x C v x C v x C v x C v x C v x C v x v x C v x C v x C v x C v x C v x C v x w x C w                                              ) ( 2 ) ( 3) ( 4 )2 3 4k k kx C w x C w x C w x      (45) V. Matveenko et alii, Frattura ed Integrità Strutturale, 49 (2019) 225-242; DOI: 10.3221/IGF-ESIS.49.23 236 where 1 2 3 4 5 6, , , , ,C C C C C C are the constants defined by the prescribed combination of boundary conditions (3) –( 5). EXAMPLE OF CALCULATION OF CHARACTERISTIC INDICES FOR CONICAL BODIES he general solutions obtained for 0k  and 1k  under the prescribed combination of the boundary conditions were used to construct a homogeneous system of algebraic equations with respect to constants iC for the examined conical body. The coefficients of this system of equations depend on the vertex angles of conical bodies, elastic characteristics of the materials and characteristic index  . The indices  , defining the character of stress singularity at the vertices of conical bodies are calculated from the condition of the existence of non-zero solution to the system of linear algebraic equations. Solid cone Let us consider a solid cone ( 00 , 0 , 0r          ). For this case, in general solutions (32), (37), (45) it is necessary to retain terms with particular solutions (15) at 0x  (or 0  ). Subject to this condition the solutions take the following form: axisymmetric rotation ( 0k  )    (1) (1)0 0 0w x C w x  (46) axisymmetric deformation ( 0k  )             (1) (1) ( 2 ) ( 2 ) 0 0 0 0 0 (1) (1) ( 2 ) ( 2 ) 0 0 0 0 0 ,u x C u x C u x v x C v x C v x         (47) non-axisymmetric deformation ( 0k  )                       (1) (1) ( 2 ) ( 2 ) ( 5 ) ( 5 ) (1) (1) ( 2 ) ( 2 ) ( 5 ) ( 5 ) (1) (1) ( 2 ) ( 2 ) k k k k k k k k k k k k k k k k k k k u x C u x C u x C u x v x C v x C v x C v x w x C w x C w x                 (48) The substitution of these solutions in one of the boundary conditions (3)–(5) at 0q q= yields a system of linear homogeneous algebraic equations for the unknown quantities, (1) ( 2 ) (1) ( 2 ) ( 5 )0 0, , , ,k k kC C C C C , which have the following representation: axisymmetric rotation ( 0k  ) (1) (1,1) 0 0 0C   (49) axisymmetric deformation ( 0k  ) (1) (1,1) ( 2 ) (1,2 ) 0 0 0 0 (1) ( 2,1) ( 2 ) ( 2,2 ) 0 0 0 0 0, 0. C C C C         (50) non-axisymmetric deformation ( 0k  ) T V. Matveenko et alii, Frattura ed Integrità Strutturale, 49 (2019) 225-242; DOI: 10.3221/IGF-ESIS.49.23 237 (1) (1,1) ( 2 ) (1,2 ) ( 5 ) (1,5 ) (1) ( 2,1) ( 2 ) ( 2,2 ) ( 5 ) ( 2,5 ) (1) ( 5,1) ( 2 ) ( 5,2 ) ( 5 ) ( 5,5 ) 0, 0, 0. k k k k k k k k k k k k k k k k k k C C C C C C C C C                   (51) A particular form of the coefficients (1,1)0 , (1,1) 0 , (1,2 ) 0 , ( 2,1) 0 , ( 2,2 ) 0 , ( 1,1) k , … ( 5,5 ) k will be specified by the selected boundary conditions, but in any case these coefficients will depend on the finite series determining the form of particular solutions ( 1)0w , (1) (1) ( 2 ) ( 2 ) 0 0 0 0, , ,u u  , (1) (1) (1) ( 2 ) ( 2 ) ( 2 ) ( 5 ) ( 5 ), , , , , , ,k k k k k k k ku w u w u   the cone angle, the elastic characteristics of the material and the sought-for characteristic index  . The condition for the existence of nontrivial solutions to the systems of the algebraic Eqns. (49) - (51) are written as (1,1) 0 0  (52) (1,1) (1,2 ) 0 0 ( 2,1) ( 2,2 ) 0 0 det 0           (53) (1,1) (1,2 ) (1,5 ) 0 0 0 ( 2,1) ( 2,2 ) ( 2,5 ) 0 0 0 ( 5,1) ( 5,2 ) ( 5,5 ) 0 0 0 det 0                    (54) From these relations we obtain the values of the characteristic indices  for each of the examined cases (axisymmetric rotation, axisymmetric deformation and non-axisymmetric deformation). When implementing numerically the algorithm for evaluation of the characteristic indices one should bear in mind that the number of terms in the series of constructed particular solutions (20), (30), (35), (36), (41) – (44) is selected in such a way that a subsequent increase in their number do not change the values of the characteristic indices in the third decimal place. It should be noted that the most “unfavorable” cases required the retention of 300-400 series terms. Fig. 2 presents the values of Re 1n  , determining singular solutions for a solid cone under the boundary conditions expressed in terms of displacements and stresses. The values given on the graphs are identical to the data presented in [19]. Figure 2: Re 1n  as a function of the vertex angle of the solid cone with boundary conditions on the lateral surface for displacements (a) and stresses (b) (▲ - 0k  ,●- 1k  ,■ - 2k  ). It should be noted that the singular solutions for a solid cone under the boundary conditions in terms of stresses take place for the zeroth, first and second harmonics of the Fourier series and under the boundary conditions in terms of displacements- for the zeroth and first harmonics of the Fourier series. V. Matveenko et alii, Frattura ed Integrità Strutturale, 49 (2019) 225-242; DOI: 10.3221/IGF-ESIS.49.23 238 All results in this paper were obtained for a material with Poisson’s ratio 0.3  . Fig. 3 displays new data on the character of stress singularity at the vertex of the solid cone for the perfect-slip boundary conditions at the lateral surface. Here singular solutions occur at the zeroth first and second harmonics of the Fourier series and at the angle 0 less than 180°. Figure 3: Re 1n  as a function of the vertex angle of the solid cone under the perfect slip boundary conditions on the lateral surface (▲ - 0k  ,●- 1k  ,■ - 2k  ). Figure 4: Variation Re 1n  with the angle 0 at fixed values of the angle 1 of the hollow cone under zero-stress boundary conditions on the lateral surface (▲ - 0k  ,●- 1k  ,■ - 2k  ). The algorithm for evaluation of the characteristic indices described above was applied to a hollow cone with two conical boundary surfaces 0  and 1  subject to different homogeneous boundary conditions (13) - (15). Note that in search for a solution to this problem we had to use a full set of particular solutions ( ( 1) ( 2 )0 0, ,w w (1) (1) ( 2 ) ( 2 ) ( 3 ) ( 3 ) ( 4 ) ( 4 ) 0 0 0 0 0 0 0 0, , , , , , , ,u u u u    (1) (1) (1), , ,k k ku w (1) (1) (1) ( 2 ) ( 2 ) ( 2 ) ( 3 ) ( 3) ( 3), , , , , , , , ,k k k k k k k k ku w u w u w   ( 4 ) ( 4 ) ( 4 ), , ,k k ku w ( 5 ) ( 5 ) ( 6 ) ( 6 ), , ,k k k ku u  ) defined by relations (20), (30), (35), (36), (41) - (44). Since the number of particular solutions doubled, the order of determinants (52) - (54) doubled, too. V. Matveenko et alii, Frattura ed Integrità Strutturale, 49 (2019) 225-242; DOI: 10.3221/IGF-ESIS.49.23 239 Fig. 4 shows the dependence of the eigenvalues Re 1n  on the cone angle formed by the external surface 0 at different values of the internal cone angle 1 . The conical surfaces meet zero boundary conditions for stresses. Here the solid line corresponds to real eigenvalues and the dashed line - to complex eigenvalues. In Fig. 5, the dependence of eigenvalues Re 1n  on 0 is plotted for two values of the internal cone angle for the case when the conical surfaces meet zero boundary conditions for displacements. Figure 5: Variation Re 1n  with the angle 0 at fixed values of the angle 1 (▲ - 0k  ,●- 1k  ) of the hollow cone under zero- stress boundary conditions on the lateral surface. For a hollow cone different combinations of the boundary conditions can be realized at the internal and external conical surfaces. We have considered two variants. In the first case, the internal surface meets the zero-stress boundary conditions and the external surface – the condition of zero displacements. In the second case the condition of zero displacements is prescribed on the internal surface and on the internal surface – the condition of zero stresses. The variation of the stress singularity index Re 1n  with the cone angle 0 formed by the external surface at different values of the internal cone angle is shown in Fig. 6. The eigenvalues, at which singular stresses occur, are observed at the values of 0 higher than 80˚. Figure 6: Re 1n  as a function of 0 for different values of the cone angle 1 at zero stresses on the internal surface and zero displacements on the external lateral surfaces (▲ - 0k  ,●- 1k  ,■ - 2k  ). For the second variant of the boundary conditions the dependence of the values of Re 1n  on 0 is shown in Fig. 7. Composite cone In the case of a composite cone the algorithm for computation of characteristic indices presented in subparagraph “Solid cone” can be applied to different versions of a composite cone. First we consider a composite cone with a single boundary conical surface 0  and contact boundary 2  . A solution to this problem is found by using regular particular solutions for the internal sub-domain 1 2(0 , , 0 2 )r            and non-regular particular solutions defined by relations ((20), (30), (35), (36), (41) - (44)) for the external sub-domain 2 0(0 , ,r        V. Matveenko et alii, Frattura ed Integrità Strutturale, 49 (2019) 225-242; DOI: 10.3221/IGF-ESIS.49.23 240 0 2 )   . As second case we consider a composite cone with two boundary conical surfaces, and contact boundary. In this case, it is necessary to use a full set of particular ( ( 1) ( 2 )0 0, ,w w (1) (1) ( 2 ) ( 2 ) ( 3 ) ( 3 ) ( 4 ) ( 4 ) 0 0 0 0 0 0 0 0, , , , , , , ,u u u u    (1) (1) (1), , ,k k ku w (1) (1) (1) ( 2 ) ( 2 ) ( 2 ) ( 3 ) ( 3) ( 3), , , , , , , , ,k k k k k k k k ku w u w u w   ( 4 ) ( 4 ) ( 4 ), , ,k k ku w ( 5 ) ( 5 ) ( 6 ) ( 6 ), , ,k k k ku u  ) defined by relations ((20), (30), (35), (36), (41)-(44)) both for the internal and external sub-domains. Figure 7: Re 1n  as a function of 0 for different values of the cone angle 1 of the hollow cone at zero displacements on the internal surface and zero stresses on the external lateral surfaces (▲ - 0k  ,●- 1k  ,■ - 2k  ). Figure 8: Re 1n  at different values of 1 2log( / )G G for 2 0 1 260 , 120 , 0.3         (a). Re 1n  at different values of Poison`s coefficient for 1 2,  for 2 060 , 120     (b). (▲ - 0k  ,●- 1k  ,■ - 2k  ). As an example, let us consider a composite cone with zero stress conditions on the external lateral boundary surface and perfect contact conditions on the boundary between different cone parts made of dissimilar materials. Fig. 8a shows variation Re 1n  of with the ratio between shear modules of different materials 1 2/G G . The dependence of eigenvalues on Poisson's ratio for 1 2/ 1G G  is given in Fig. 8b. Solid lines denote eigenvalues for 1 0.3  as a function of 2 and dash and dot lines denote eigenvalues for 2 0.3  as a function of 1 . Of particular interest is a composite cone with the internal singular point. In this case, two types of conditions can be prescribed at the boundary between two different materials 2  (perfect bonding condition (8) or frictionless slip conditions (9). The stress singularity indices shown in Fig. 9 correspond to condition (8) for different shear modules ratios of different cone layers. It has been found that at the internal singular point the condition of stress singularity is realized at any ratio between shear modules of the inner and outer cones, except for the trivial case of 2 90   . V. Matveenko et alii, Frattura ed Integrità Strutturale, 49 (2019) 225-242; DOI: 10.3221/IGF-ESIS.49.23 241 Figure 9: Variation of Re 1n  with 1 2log( / )G G in condition of perfect contact between two different materials at 2 1 260 , 0.3      . (▲ - 0k  ,●- 1k  ,■ - 2k  ). CONCLUSION e have considered the analytical method of constructing eigensolutions for circular cones. It has been shown that the proposed analytical relations can be used to construct solutions and to evaluate the character of stress singularity for different conical bodies (solid, composite) and different types of boundary conditions on the lateral surfaces and contact surface of dissimilar materials. The numerical results presented in this paper provide information about the character of singular stresses at the vertex of a solid cone for displacement, stress and combined boundary conditions and at the vertex of a hollow cone for different boundary conditions on the internal and external lateral surfaces. The importance of the theoretical output of this research is that the obtained analytical solutions provide full information about singular solutions for different types of circular conical bodies made of isotropic elastic materials. All formulas and options for calculation, allowing readers to obtain independently numerical results for the problems considered in the paper, can be found on E-resource (www. icmm.ru/compcoeff). The results of this work may have considerable practical implications in various fields. Recent advances in the development of computational methods extend the possibilities of the in-depth analysis of the stress state, so that present computations of real objects provide a great deal of information on the stresses in the vicinity of singular points, which are often the stress concentrator zones. A search for effective ways of reducing the level of stress concentration, for example, by changing the geometry in the vicinity of singular points leads to the conclusion that the variants close to optimal are related with the character of singular solutions [22]. Such conclusions have found commercial applications, for example, for the development of techniques allowing engineers to increase the strength of adhesive or glue [23-25]. In this respect, the results presented in the paper can be used to find an appropriate variant of reducing the level of stress concentration in the vicinity of non-smooth points on the surface of three-dimensional elastic bodies. The obtained information has its own practical value since it allows us to estimate the character of stress behavior at the vertices of different conical natural objects and engineering structures. Moreover, the available data on the indices of stress singularity open up new possibilities for construction of singular finite elements [22]. REFERENCES [1] Kondratiev, V.A. (1967). Boundary value problems for elliptic equations in the regions with conical and angular points. Transactions of the Moscow Mathematical Society 16, pp. 209-292. [2] Sinclair, G.B. (2004). Stress singularities in classical elasticity - I: Removal, interpretation, and analysis. Applied Mechanics Reviews 57(4), pp. 251-297. [3] Sinclair, G.B. (2004). Stress singularities in classical elasticity - II: Asymptotic identification. Applied Mechanics Reviews 57(5), pp. 385-439. W V. Matveenko et alii, Frattura ed Integrità Strutturale, 49 (2019) 225-242; DOI: 10.3221/IGF-ESIS.49.23 242 [4] Mikhailov, S.E. (1979). Stress singularity in the vicinity of a rib of inhomogeneous composite anisotropic body and some applications to composites. Proceedings of the Academy of Sciences of the USSR. Solid Mechanics. 5, pp. 103- 110. [5] Apel, T., Mehrmann, V., Watkins, D. (2002). Structured eigenvalue methods for the computation of corner singularities in 3D anisotropic elastic structures. Computer Methods in Applied Mechanics and Engineering 191, pp. 4459-4473. DOI: 10.1016/S0045-7825(02)00390-0 [6] Babeshko, V.A., Glushkov, E.V., Glushkova, N.V., Lapina, O.N. (1991). Stress singularity in the vicinity of the vertex of elastic trihedral. Proceedings of the Russian Academy of Sciences 318(5), pp. 1113-1116. [7] Dimitrov, A., Andra, H., Schnack, E. (2002). Singularities near three-dimensional corners in composite laminates. International Journal of Fracture 115, pp. 361-375. DOI: 10.1002/nme.230 [8] Matveenko, V.P., Nakaryakova, T.O., Sevodina, N.V., Shardakov, I.N. (2006). Investigation of Singularity of Stresses at the Apex of a Cone with an Elliptic Base. Doklady Physics 51(11), pp. 630-633. DOI: 10.1134/S1028335806110140. [9] Matveenko, V.P., Nakaryakova, T.O., Sevodina, N.V., Shardakov, I.N. (2008). Stress singularity at the vertex of homogeneous and composite cones for different boundary conditions. Journal of Applied Mathematics and Mechanics 72, pp. 331-337. DOI: 10.1016/j.jappmathmech.2008.07.012. [10] Korepanova, T.O., Matveenko, V.P., Sevodina, N.V. (2010). Numerical analysis of stress singularity at the vertex of a cone with unsmooth lateral surface. Computer Mechanics of continuum Mechanics 3(3), pp. 68-76. [11] Korepanova, T.O., Matveenko, V.P., Sevodina, N.V. (2013). Numerical analysis of stress singularity at singular points of three-dimensional elastic bodies. Acta Mechanica 224, pp. 2045-2063. DOI 10.1007/s00707-013-0845-y. [12] Savruk, M.P. and Shkaraev, S.V. (2001). Stress singularities for three-dimensional corners using the boundary integral equation method. Theoretical Applied Mechanics 36, pp. 263-275. [13] Bazant, Z.P., Keer, L.M. (1974). Singularities of elastic stresses and of harmonic functions at conical notches or inclusions. International Journal of Solids and Structures 10(9), pp. 957-965. [14] Beagles, A.E. and Sänding, A.M. (1991). Singularities of rotationally symmetric solutions of boundary value problems for the Lame equations. ZAMM. Journal of Applied Mathematics and Mechanics. 71(11), pp. 423-431. DOI: 10.1002/zamm.19910711102. [15] Bernardi, C., Dauge, M., Maday, Y. (1999). Spectral Methods for Axisymmetric Domains. Series in Applied Mathematics 3. [16] Parihar, K.S., Keer, L.M. (1978). Elastic stress singularities at conical inclusions. International Journal of Solids and Structures 14, pp. 261-263. DOI: 10.1016/0020-7683(78)90036-7. [17] Picu, C.R. (1996). Stress singularities at vertices of conical inclusions with freely sliding interfaces, International Journal of Solids and Structures 33(17), pp. 2453-2457. DOI: 10.1016/0020-7683(95)00164-6. [18] Somaranta, N., Ting, T.C.T. (1986). Three dimensional stress singularities at conical notches and inclusions in transversely isotropic materials, Journal of Applied Mechanics 53, pp. 89-96. [19] Kozlov, V.A., Maz'ya, V.G., Rossmann, J. (2001). Spectral Problems Associated with Corner Singularities of Solutions to Elliptic Equations. Providence RI.N.J.: American Mathematical Society, ser. Mathematical Surveys and Monographs 85. DOI: http://dx.doi.org/10.1090/surv/085. [20] Nowacki, W. (1970). Theory of Elasticity. PWN. Warsaw. [21] Kamke, E. (1959). Differentialgleichungen: Lösungsmethoden und Lösungen. Gewöhnliche Differentialgleichungen, Leipzig. [22] Matveyenko, V.P. and Borzenkov, S.M. (1996). Semianalytical singular element and its application to stress calculation and optimization. International Journal for Numerical Methods in Engineering 39(10), pp. 1659-1680. DOI: 10.1002/(SICI)1097-0207(19960530)39 [23] He, D., Sawa, T., Karami, A. (2009). Stress analysis and strength evaluation of scarf adhesive joints with dissimilar adherends subjected to static tensile loadings, International Journal of Adhesion and Adhesives 3(8), pp. 1033-1044. [24] Lang, T.P., Mallick, P.K. (1998). Effect of spew geometry on stresses in single lap adhesive joints. International Journal of Adhesion and Adhesives 18, pp. 167-177. DOI: 10.1016/S0143-7496(97)00056-0. [25] Xu, L.R., Kuai, H., Sengupta, S. (2004). Dissimilar material joints with and without free-edge stress singularities: Part II. An integrated numerical analysis. Experimental mechanics 44(6), pp. 616-621. DOI: 10.1007/BF02428251. << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice