J\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{J}$$\end{document}-, M- and L\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{L}$$\end{document}-integrals of line charges and line forces

In this work, the J\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{J}$$\end{document}-, M- and L\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{L}$$\end{document}-integrals of two line charges and two line forces in generalized plane strain are derived in the framework of three-dimensional (3D) electrostatics and three-dimensional compatible linear elasticity, respectively, in order to study the interaction between them. The key point in this derivation is achieved by expressing the J\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{J}$$\end{document}-, M- and L\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{L}$$\end{document}-integrals of point charges and point forces in three dimensions in terms of the corresponding three-dimensional Green functions, that is, the three-dimensional Green function of the Laplace operator and the three-dimensional Green tensor of the Navier operator, respectively. The major mathematical tool used in deriving the J\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{J}$$\end{document}-, M- and L3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L_3$$\end{document}-integrals of line charges and line forces from the corresponding J\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{J}$$\end{document}-, M- and L\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{L}$$\end{document}-integrals of point charges and point forces is the method of embedding or method of descent of Green functions to two dimensions (2D) from the corresponding Green functions in 3D. The analytical expressions of J\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{J}$$\end{document}-, M- and L3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L_3$$\end{document}-integrals of line charges and line forces in antiplane and plane strain are derived and discussed. The J\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{J}$$\end{document}-integral is the electrostatic part of the Lorentz force (electrostatic interaction force) between two line charges in electrostatics and the Cherepanov force between two line forces in elasticity. The M-integral of two line sources (charges or forces) equals half the electrostatic interaction energy in electrostatics and half the elastic interaction energy in elasticity of these two line sources, respectively. The L3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L_3$$\end{document}-integral of two line sources (charges or forces) is the z-component of the configurational vector moment or the rotational moment representing the total torque about the z-axis caused by the interaction of the two line sources. An important outcome is that the J\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{J}$$\end{document}-integral is twice the negative gradient of the M-integral, leading to the result that the J\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{J}}$$\end{document}-integral for line charges and line forces is a conservative force (with the corresponding interaction energy playing the role of the potential energy) and consequently an irrotational vector field. Finally, the obtained J\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{J}$$\end{document}-, M- and L3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L_3$$\end{document}-integrals being functions of the distance and of the angle, are able to describe the physical behaviour of the interaction of two line sources (charges and forces), since they represent the fundamental and necessary quantities, that is, the interaction force, interaction energy and total torque produced by the interaction of these two line sources.

The corresponding conservation integrals are the J-, Land M-integrals introduced by Budiansky and Rice [11]. Budiansky and Rice [11] were the first to give a physical interpretation to the J-, Land M-integrals as the energy-release rates per unit cavity translation, rotation and expansion, respectively. The J-integral for cracks is widely used and known as the J-integral of Rice [39], which is the driving force acting on the crack tip and it is a material or configurational force [33]. The J-integral of dislocations is equivalent to the well-known Peach-Koehler force that is the interaction force between two dislocations [23,40,45] and it is also a configurational force. In the framework of configurational or Eshelbian mechanics, Lazar and Kirchner [29] showed that, in the general case of anisotropic linear first strain gradient elasticity for non-homogeneous and incompatible media with dislocations and body forces present, the J-integral gives the sum of configurational forces of different origin, namely the Peach-Koehler force for dislocations, the Cherepanov force due to body forces and the Eshelby force due to inhomogeneities, the M-integral is equivalent to the total configurational or material work and the L-integral is equivalent to the sum of configurational or material vector moments also of different origin (see also [3]). An extensive introduction about the J-, M-and L-integrals can be found in Agiasofitou and Lazar [2], Lazar and Agiasofitou [28] and the references therein. A recent monograph devoted in the field of configurational mechanics, presenting among others, various vistas of the material or configurational forces and highlighting its advantages to tackle problems in continuum defect mechanics, providing a dissipation-consistent approach, has been given by Steinmann [42].
Special attention should be paid to the M-integral and its physical interpretation. Rice [40] was the first to raise the question "What is M for a dislocation?" and studying the M-integral when centered on a dislocation line in the framework of two-dimensional, compatible linear elasticity, he found that the M-integral equals to the "dislocation energy factor". In fracture mechanics, first Chen [12] wrote explicitly that the M-integral equals twice the change of the total potential energy owed to single cracking of a central crack in a plane elastic body. In the framework of three-dimensional incompatible elasticity of dislocations, the physical interpretation of the M-and L-integrals for dislocations was given for the first time in Agiasofitou and Lazar [2]. Namely, the M-integral between two straight dislocations (per unit dislocation length) is equal to the half of the interaction energy of the two dislocations (per unit dislocation length) depending on the distance and on the angle, plus twice the corresponding pre-logarithmic energy factor and the L 3 -integral between two straight dislocations is the z-component of the configurational vector moment or the rotational moment about the z-axis caused by the interaction of the two dislocations [2]. Moreover, the J-, M-and L-integrals of body charges and point charges in electrostatics, and the J-, M-and L-integrals of body forces and point forces in elasticity have been given by Lazar and Agiasofitou [28]. In particular, it is proven that the J-integral of body charges in electrostatics represents the electrostatic part of the Lorentz force, and the J-integral of body forces in elasticity represents the Cherepanov force, both of them are configurational forces. The M-integral of two point sources (charges or forces) equals half the electrostatic interaction energy in electrostatics and half the elastic interaction energy in elasticity of these two point sources. The L-integral represents the configurational vector moment or total torque between two body or point sources (charges or forces). The physical interpretation of the M-integral given above in [2,28] is in full agreement with the physical interpretation of the M-integral in fracture mechanics as we have discussed above given in [12].
It becomes clear that, in general, the J-integral is equivalent to the interaction force (electrostatic part of the Lorentz force in electrostatics, Peach-Koehler force of dislocations, Cherepanov force of point forces) and it is a configurational force. The M-integral is interaction energy (electrostatic interaction energy in electrostatics, interaction energy between two dislocations, elastic interaction energy between two point forces in elasticity) and it is not representing a "force" as misleadingly stated in the literature as "self-similar expansion force" [41] or "generalized force" [17]. Finally, the L-integral represents the total torque and it is a configurational vector moment [2,27,28]. This physical interpretation of the examined integrals has been also shown to be true in more complicated cases such as for electro-elastic dislocations in piezoelectric materials [3] and even for the complex structure of quasicrystals [26].
A crucial point that is noteworthy to be highlighted and which can shed light to the misunderstandings around the physical interpretation of the M-integral is that the analytical expression of the M-integral for a two-dimensional body is different than that for a three-dimensional body. This happens because the M-integral depends on the space dimension (see related discussion in [2]), since it is inherently related to the group of scaling transformations and the scaling or canonical dimension of the field variable of elasticity, that is of the displacement field, is different in 2D than that in 3D (see, e.g., [2,28]). For physical phenomena such as straight dislocations and line forces that have to be treated in three dimensions (so n = 3) due to their physics, if one reduces the study of these problems to two dimensions (so n = 2), then the M-integral is just a constant (see, e.g., [4,5,10,20,40,41]) loosing the important contribution which offers the functional dependence on the distance and on the angle between the examined defects or sources like dislocations, charges, point forces, etc. (see [2,28]).
A line force is a line singularity in a three-dimensional medium. Lazar and Maugin [31] found closed analytical solutions for the displacement fields of line forces, free of singularities, in the framework of simplified first strain gradient elasticity. In addition, generalized plane strain is a two-dimensional deformation in a threedimensional body (see, e.g., [34]). Therefore, the problem of the interaction of two line forces has to be studied in the framework of three-dimensional compatible linear elasticity, considering for the calculation of the Mintegral the space dimension equals three (n = 3). Regarding previous results in the literature as far as the J-, M-and L 3 -integrals of line forces are concerned, the J-, M-and L 3 -integrals of two-dimensional elastic line singularities (among them also line forces) in the framework of plane strain linear elasticity have been derived by Seo et al. [41] but only in a marginal case without providing the general analytical expressions of the examined integrals. Moreover, the M-integral in [41] being a constant (since it has been considered n = 2) cannot capture the interaction between the two line forces in contrast to the M-integral derived in this work (for n = 3), which is a function of the distance and of the angle and in this way the physics of the physical phenomenon can be appropriately described (see related discussion in the Appendix A). It is important to pay attention to this matter, since the results and statements in [41] about the M-integral of line forces have been taken over in a recent review article [20] about the M-integral.
This article is organized as follows. In Sect. 2, the J (AB) -, M (AB) -and L (AB) -integrals of point charges and point forces in 3D are given for the first time in terms of the three-dimensional Green function G (3) (R) of the Laplace operator and the three-dimensional Green tensor G (3) (R) of the Navier operator, respectively, giving the opportunity to use the method of embedding or method of descent of Green functions to 2D from 3D in Sect. 3  -integrals are given in terms of the corresponding pre-logarithmic energy factor of line forces which is introduced here for the first time in analogy to the pre-logarithmic energy factor of dislocations. The result that the J (AB) -integral is twice the negative gradient of the M (AB) -integral gives rise to prove that the J (AB) -integral for line charges and line forces is a conservative force in Sect. 6. The results are summarized in the conclusions in Sect. 7. The Appendix A serves to make the necessary comparisons with the already given results in the literature and to highlight important aspects concerning the M (AB) -integral.

J-, M-and L-integrals of point charges and point forces in 3D
In this Section, we present in brief the J (AB) -, M (AB) -and L (AB) -integrals of point charges in electrostatics and point forces in elasticity in 3D that have been derived in Lazar and Agiasofitou [28]. We recall only the results that are necessary for the purpose of this work. The reader is addressed to [28] for an extensive presentation and detailed derivations. In particular, the J (AB) -, M (AB) -and L (AB) -integrals of point charges and point forces in 3D with the "source field" applied at the origin derived in [28], are used and further expressed in terms of the corresponding three-dimensional Green functions providing the appropriate expressions for the forthcoming calculations.

Electrostatics: point charges in 3D
Let us consider an unbounded, linear, homogeneous and isotropic medium. In electrostatics, the field variable is the electrostatic potential ϕ = ϕ(x). The electric field strength vector E k = E k (x) is defined as the negative gradient of the electrostatic potential ϕ: if the electrostatic Bianchi identity, which has the meaning of electrostatic compatibility condition, is fulfilled where i jk is the Levi-Civita tensor and ∂ j denotes partial differentiation with respect to the spatial coordinates x j , j = 1, 2, 3.
Considering the interaction between an electric field E i (x) or the corresponding electrostatic potential ϕ (B) = ϕ (B) (x) as source field and a body charge density q (A) = q (A) (x) as receiver field, then the J (AB) k -, M (AB)and L (AB) k -integrals for continuous body charge density distributions are given by [28] where i, j, k = 1, 2, 3. We recall here that the J-integral of electrostatics represents the electrostatic force, which is the electrostatic part of the Lorentz force since is the electrostatic part of the Lorentz force density. The latter is the electrostatic force density on a body charge density q in presence of an electric field E k . Therefore, it can be seen in Eq.
It should be noted that for the derivation of the M (AB) -integral given by Eq. (4), the space dimension has been taken equal to three, n = 3 (see [28]). We consider now that the body contains two point charge densities, q (A) located at point A with position vector x (A) and q (B) located at point B with position vector x (B) . The point B is the "source-point" (point of sender) where the charge density q (B) , playing the role of a "source charge", is located; and A is the "field-point" (point of receiver) where the charge density q (A) , playing the role of a test charge measuring the electrostatic force, energy and torque, is located. The two body charge densities can be modelled in terms of the point chargesq (A) andq (B) and the three-dimensional Dirac delta function δ(x) as follows In the case that the interaction is between a point charge (source charge)q (B) located at the origin of the coordinate system O, that is x ( where R = |x (A) | = x 2 +ȳ 2 +z 2 and ε denotes the permittivity of material or dielectric constant of matter. It can be seen in Eq. (11) that the electrostatic J (AB) -integral is the electrostatic force F L(AB) acting on the point chargeq (A) that is placed at the position x (A) due to the point chargeq (B) placed at the position x (B) ≡ O and it is provided by the famous Coulomb law, which is the experimentally established law for the electrostatic part of the Lorentz force between two point charges in isotropic media. The electrostatic force along the line from B to A is repulsive ifq (B) andq (A) have the same sign, and attractive if their signs are opposite. Equation (12) shows that the M (AB) -integral is half the electrostatic interaction energy U (AB) of the two point charges (see, e.g., [21,28]). Equation (13) gives a zero value for the L (AB) -integral, which is reasonable since there is no torque between two scalar fields.
The important observation that can be made at this point is that the non-vanishing J (AB)and M (AB)integrals of point charges in 3D, Eqs. (11) and (12), can be written in terms of the three-dimensional Green function G (3) (R) of the Laplace operator as follows where the three-dimensional Green function of the Laplace operator is given by (see, e.g., [9,28])

Elasticity: point forces in 3D
Let us consider an unbounded, linear, homogeneous and isotropic elastic body. In elasticity, the field variable is the displacement vector u = u(x). The elastic distortion tensor β ik = β ik (x) can be written as gradient of a displacement vector field u i : if the compatibility condition of elasticity jlk ∂ l β ik = 0 is fulfilled. Considering the interaction between an elastic distortion β (B) ik (x) or the corresponding displacement field u (B) i (x) as source field and a body force density f (A) i (x) as receiver field, then the J (AB) -, M (AB)and L (AB) -integrals in 3D for continuous body force density distributions are given by [28] In Eq. (18), it can be seen that the J To recall, the Cherepanov force is given by where is the so-called Cherepanov force density [13], which is the configurational force density acting on a body force density f i in presence of an elastic distortion β ik .
The body force densities f can be modelled with the point forcesf where μ is the shear modulus and ν is the Poisson ratio. Here, U (AB) denotes the elastic interaction energy and T (AB) k denotes the total torque vector of the two point forces. It has to be noticed that for the derivation of the M (AB) -integral in Eqs. (19) and (27), the space dimension has been taken equal to three, n = 3 (see [28]).
The crucial step here is to observe that in the above expressions of the J i j (R) of the Navier operator by the following simple formulas: where the three-dimensional Green tensor of the Navier operator is given by (see, e.g., [28,32,35])

Method of embedding or method of descent to 2D from 3D
Line charges and line forces can be viewed as infinitely long line sources along the entire z-axis in 3D. By symmetry, the problem of line sources is independent ofz and is a 2D problem which is embedded in the 3D space (see, e.g., [9]). Therefore, the solution of line sources depending on (x,ȳ) can be obtained from the solution of point sources depending on (x,ȳ,z) using the projection from 3D to 2D. In this manner, the J-, M-and L-integrals of line charges and line forces can be constructed from the J-, M-and L-integrals of point charges and point forces using the method of embedding [9] or the so-called method of descent following Hadamard [19], because one is descending to 2D from 3D via the following relations [9,43] Eqs. (33)- (35) give the projection of the three-dimensional J k (x,ȳ,z)-, M(x,ȳ,z)-and L k (x,ȳ,z)-integrals onto thexȳ-plane and these are the J k (x,ȳ)-, M(x,ȳ)-and L k (x,ȳ)-integrals per unit length l z embedded in the 3D space. For reasons of brevity, per unit length l z is omitted to be written explicitly in the formulas of the J k (x,ȳ)-, M(x,ȳ)-and L k (x,ȳ)-integrals. However, it is implied that it is per unit length l z . It has to be noticed that the second term in the M-integral in Eqs. (4) and (19) arises by considering the space dimension 1 n = 3 and survives without being influenced by the projection, which is very reasonable since line sources "live" in 3D.  29), (30) and (31), in terms of the three-dimensional Green tensor function G (35) is reduced to the projection or descending of the corresponding Green functions from 3D to 2D, which is a well-known technique in the theory of partial differential equations (see, e.g., [9,43]) and it will be applied in the forthcoming Sects. 4 and 5.
The method of projection from 3D to 2D has also been used in dislocation theory for the calculation of solutions of straight dislocations in isotropic materials (see, e.g, [6,35]), whereas the method of the projectionslice theorem has been used for dislocations in anisotropic materials in [30].

J-, M-and L 3 -integrals of line charges
In this Section, we consider two infinitely long straight line charges parallel to z-direction in an unbounded, homogeneous and isotropic body.
Substituting the J (AB) where is the two-dimensional Green function of the Laplace operator. Using Eq. (16), Eq. (39) is explicitly written (see, e.g., [9]) wherer = x 2 +ȳ 2 and L denotes the length of the line charge in z-direction. Note that the argument of a logarithm must be dimensionless (see, e.g., [9]) and this is guaranteed by dividingr with the length L in the two-dimensional Green function (40) due to the projection from 3D to 2D. In fact, the limits of the integral from −∞ to ∞ in Eq. (39) are replaced by finite numbers from −L/2 to L/2 and the final expression (40) is valid for large L/r (see [9]). From a mathematical point of view, L is a cut-off parameter.  14), (15) and (13), with the corresponding ones of line charges, Eqs. (36), (37) and (38), it is important to make the following observation. In electrostatics, the projection of the J (AB) k -, M (AB) -and L (AB) k -integrals from 3D to 2D is reduced to the projection of the Green function of the Laplace operator from 3D to 2D: leading in this way to a known and standard quantity in the theory of partial differential equations, G (2) (x,ȳ), (see, e.g., [9,43]). Indeed, substituting the two-dimensional Green function (40) Here,r = x 2 +ȳ 2 is the distance between the two line charges (see Fig. 1). Equations (41) and (42) give the non-vanishing components of the electrostatic part of the Lorentz force F L(AB) , which is the electrostatic interaction force between the two line charges and is plotted in Fig. 2. It can be seen that the electrostatic force F L(AB) is repulsive ifq (A)q(B) > 0 ( Fig. 2(a)) and is attractive ifq (A)q(B) < 0 ( Fig. 2(b)). The M (AB)integral, Eq. (43), equals half the electrostatic interaction energy (per unit length) of the two line charges where the electrostatic interaction energy is given by (see [28]) with the electrostatic potential of a line charge In a cylindrical coordinate system with basis (e r , e ϕ , e z ), the two components J It can be seen that the electrostatic force between two parallel line charges is a purely radial force, since F which means that the electrostatic force is a central force (in 2D). If

J-, M-and L 3 -integrals of line forces
In this Section, we consider two infinitely long straight line forces (parallel to z-direction) of uniform strength per unit length in an unbounded, homogeneous and isotropic elastic body. Substituting the J where is the two-dimensional Green tensor of generalized plane strain. Using Eq. (32), we can obtain from Eq. (57) the two-dimensional Green tensor function of plane strain (see, e.g., [35]) and the two-dimensional Green function of antiplane strain (i = j = 3) where ( k -integrals from 3D to 2D is reduced to the projection of the three-dimensional Green tensor function to the two-dimensional Green tensor function of generalized plane strain: G ,ȳ), i, j = 1, 2, 3 leading to known and standard quantities in the mathematical theory of elasticity: G (2) i j (x,ȳ), i, j = 1, 2 and G (2) 33 (x,ȳ) (see, e.g., [35]). Note that in the projection: G (3) 33 (x,ȳ,z) −→ G (2) 33 (x,ȳ), an unimportant constant term has been neglected which may be disregarded because it satisfies the homogeneous Laplace equation (see [43]).
In  (56). It should be clarified here that the line direction is the z-direction but the strength can be in x-, y-or z-direction.

Line forces in antiplane strain
We examine first the case of two parallel line forces with strengthsf Eqs. (60) and (61) give the non-vanishing components of the Cherepanov force F C(AB) (per unit length l z ) of the two line forces, which is plotted in Figs. 3(a) and 4(a). It can be seen that the Cherepanov force F C(AB) is attractive iff > 0 ( Fig. 3(a)) and is repulsive iff < 0 ( Fig. 4(a)). The M (AB) -integral, Eq. (62), equals half the elastic interaction energy U (AB) per unit length l z of the two line forces where the elastic interaction energy (per unit length) for line forces in antiplane strain is given by (see [28]) with the corresponding displacement field In Eq. (63), we find that the L (AB) 3 -integral is zero, which is an expected result since for two parallel line forces with strengths in z-direction no torque T (AB) 3 can be produced. In analogy to dislocation theory where the pre-logarithmic energy factor is a well-known quantity (see, e.g., [7,8,40]), we introduce here the prefactor of the logarithmic term in the M (AB) -integral, Eq. (62), as the pre-logarithmic energy factor of line forces with strengths in z-direction: The matrix in an orthogonal coordinate system basis (m, n, τ ) so that the z-axis (and therefore the line of the force) is parallel to τ = m × n, is well-known in the integral formalism of anisotropic elasticity (see, e.g., [7,8]). In isotropic elasticity, Q i j is a negative definite diagonal matrix and reads as Therefore, Eq. (62) can be rewritten as In cylindrical coordinates, the two components J  It can be seen in Eqs. (71) and (72) that the Cherepanov force between two parallel line forces in z-direction in antiplane strain is a purely radial force: This means that in the examined case, the Cherepanov force is a central force (in 2D).
> 0, that is, if the forces are parallel. This behaviour is similar to magnetostatics, where the magnetostatic part of the Lorentz force between two parallel wires carrying the currents I 1 and I 2 is attractive if the currents flow in the same direction and repulsive if the currents flow in opposite directions (see, e.g., [21]).
To sum up, the J (AB) k -, M (AB) -and L (AB)

3
-integrals for parallel line forces in antiplane strain are given by or in terms of the corresponding pre-logarithmic energy factor K (AB) 33 read as

Line forces in plane strain
We continue with the case of two parallel line forces whose strengthsf where i, j, k = 1, 2. It is interesting to observe that both M (AB) -and L is given by (see, e.g., [28]) where the displacement vector of the line forcef Therefore, the M (AB) -integral in Eq. (83) equals half of the elastic interaction energy (per unit length) of the two line forces It is interesting to note that the elastic interaction energy of the two line forces, Eq. (85), has a similar form to the elastic interaction energy (per unit length) of two parallel edge dislocations with general Burgers vectors given in [25,36].
Regarding the physical interpretation of the L-integral, Agiasofitou and Lazar [3] studying electro-elastic dislocations in piezoelectric materials considering additionally inhomogeneities in presence of body forces and body charges showed that the L-integral is the total configurational or material vector moment representing the total torque, since it is able to capture every produced configurational vector moment density: the configurational vector moment density produced by the total configurational force density plus an intrinsic or field vector moment density due to body force density vector and the eigendistortion tensor plus the configurational vector moment density due to the inherent material anisotropy of the piezoelectric materials. Hence, the L (AB) 3 -integral, Eq. (84), represents the total torque T (AB) 3 of two parallel line forces in plane strain and it can be decomposed into its orbital and intrinsic parts as follows 2 The orbital part, which we call the orbital torque, is given as torque produced by the Cherepanov force and the intrinsic (or spin) part, which we call the intrinsic torque, is due to the displacement vector field u (B) i (which is a spin 1 field) and the line forcef If we substitute Eq. (82) into Eq. (89), the orbital torque reads as and if we substitute Eq. (58) into Eq. (90), the intrinsic torque reads as Remark 1 In order to gain more insight into the L-integral and the above representation, it is interesting to recall its derivation based on Noether's theorem as it has been done in Agiasofitou and Lazar [1] using the prolongation method of Olver [37], where it can be seen that the L-integral arises from the symmetry group of rotations of both, independent and dependent variables in space, giving rise to two kinds of torque, the orbital torque and the intrinsic one, respectively. The intrinsic torque is related to the vectorial nature of the displacement field.
In the following, we examine all possible combinations of the directions of strengths, applying Eqs. -integrals for plane strain.

Line forces with strengths in x-direction
We consider first the case of two parallel line forces with strengthsf The J  (94), give the non-vanishing components of the Cherepanov force F C(AB) (per unit length l z ) of the two line forces with strengths in x-direction, which is plotted in Figs. 3(b) and 4(b). It can be seen that the Cherepanov force F C(AB) is attractive iff > 0 (Fig. 3(b)) and repulsive iff < 0 (Fig. 4(b)). The M (AB) -integral (95) gives half the interaction energy (per unit length l z ) of the two line forces with strengths in x-direction and is plotted in Figs. 5(a) and 6(a). The is applied in one of the coordinate axis then no torque is produced. Using the decomposition of the total torque (88) into the orbital part (91) and the intrinsic part (92), the orbital torque for the case under study reduces to and the intrinsic torque becomes In cylindrical coordinates, the J In Eq. (99), it can be seen that the radial force component F are antiparallel and is attractive iff are parallel. Moreover, it is easy to see that due to the non-vanishing tangential force component F C(AB) ϕ , Eq. (100), the Cherepanov force between two parallel line forces with strengths in x-direction is not a purely radial force in contrast to the case of antiplane strain. In addition, F -integral, using its decomposition to orbital and spin parts, it can be seen that Therefore, it turns out that the L (AB) 3 -integral of two parallel line forces with strengths in x-direction provides the total torque, which is produced finally only by the tangential Cherepanov force F C(AB) ϕ , and the rotational direction of the line forcef = 0 if ϕ = {0, π/2, π, 3π/2}, which means that there is no torque at these positions as one can also see in Figs. 5(b) and 6(b).
The prefactor of the logarithmic term in the M (AB) -integral, Eq. (101), is the pre-logarithmic energy factor of line forces with strengths in x-direction: Eqs. (99)-(102) can be rewritten in terms of the pre-logarithmic energy factor K (AB) 11 , Eq. (104), as follows

Line forces with strengths in y-direction
The next case that has to be examined is the case of two parallel line forces with strengthsf The J It is easy to see that the Cherepanov force between two line forces with strengths in y-direction is not a purely radial force due to the non-vanishing component F 2 . It can be easily checked that the relation (103) also holds for parallel line forces with strengths in y-direction. Indeed, and finally Therefore, the L (AB) 3 -integral for two parallel line forces with strengths in y-direction provides the total torque and it is produced only by the tangential Cherepanov force

Line forces with perpendicular strengths
The last case is that of two parallel line forces with strengthsf The J -integral, Eq. (129), represents the total torque about the z-axis caused by the interaction of the two line forces. Using the decomposition of the total torque (88) into its orbital part (91) and its intrinsic part (92), the orbital torque reduces to and the intrinsic torque becomes In cylindrical coordinates, the J  < 0 (see also Fig. 4(d)). 3 On the other hand, the M (AB) -integral is just a function of sin 2ϕ and it no longer contains a logarithmic part as one can see in Eq. (134). Moreover, it is remarkable that the case of line forces with perpendicular strengths is the only case where the M (AB) -integral representing the elastic interaction energy does not include a logarithmic part and for that reason there is no pre-logarithmic energy factor, K (AB) 12 = 0. The latter is in accordance with the formula (137) that we will see in the next subsection, due to the fact that Q 12 = 0. Finally, in the case of line forces with perpendicular strengths the biggest contribution is given by the L (AB) 3 -integral, 3 Note that in these positions the orbital torque M (AB) 3 =r F C(AB) ϕ produced by the tangential Cherepanov force F C(AB) ϕ is also zero (see also Eq. (136)).
Eq. (135), and this is reasonable since due to the geometry of the line forces, the produced total torque is significant.
Furthermore, it is interesting to decompose the L (AB) 3 -integral in its orbital and intrinsic parts as follows In Eq. (136), it can be seen that the L (AB) 3 -integral of two parallel line forces with perpendicular strengths provides the total torque, which is not only produced by the tangential Cherepanov force F C(AB) ϕ as in the previous cases, since additional terms even of logarithmic type appear. Note that in this case, M (r , ϕ) depends not only on ϕ but also onr . -integrals of line forces in antiplane and in all cases of plane strain are functions of the radiusr and angle ϕ capturing the interaction between the two line forces. It is significant to emphasize this result concerning the M-integral, since it is in contrast to the results and statements of Seo et al. [41], where the M-integral is a constant and not a function ofr and ϕ and consequently "do not recognize any elastic field produced by other interacting defects", which is not true as it is shown in this paper (see also Appendix A).
which shows that the J (AB) -integral is nothing but twice the negative gradient of the M (AB) -integral On the other hand, it has been shown throughout the paper that for line charges and for all cases of line forces it holds also the relation which in combination with Eq. (148) yields Eq. (150) provides the important result that the J (AB) -integral representing the electrostatic part of the Lorentz force F L(AB) of two line charges in electrostatics and the Cherepanov force F C(AB) of two line forces in elasticity is a conservative force and the corresponding (electrostatic or elastic) interaction energy U (AB) plays the role of the potential energy for the corresponding force in electrostatics and elasticity, respectively. As a consequence, the curl of J (AB) -integral is zero, that is, which means that the J (AB) -integral is an irrotational (curl-less or curl-free) vector field. Equations (143)-(151) are valid for line charges and for all cases of line forces and they are fundamental relations. It is of great significance to notice that the fundamental relations between the J (AB) -and M (AB)integrals given in this work for line charges and line forces in isotropic case are also valid for straight screw and edge dislocations with parallel or antiparallel Burgers vectors in isotropic materials as found in [2]; for parallel electro-elastic screw dislocations in hexagonal piezoelectric materials as showed in [3]; and for parallel screw dislocations in one-dimensional hexagonal quasicrystals, even if they have a complicated structure, as found in [26].
Remark 5 It is noted that for line forces in plane strain with parallel or antiparallel strengths, the following relations also hold between the L It should be emphasized that for line forces with perpendicular strengths, the relations (152) and (153) are not valid.

Conclusions
In this work, we have determined the J ( -integrals of line forces in plane strain are given for (x,ȳ) = (a, 0) with a > 0 in order to make the necessary comparisons with already given results in the literature concerning line forces (see Seo et al. [41]). Obviously, this example is marginal and restricted, since in this case the angle ϕ is zero and one looses the influence of the angle dependence, which as it is shown here is important. However, only this case has been considered in [41]. The obtained results are listed in Table 2.
Comparing the obtained results in Table 2 and the corresponding results in Table 2 in Seo et al. [41], the following noteworthy observations can be made: • Concerning the M-integral, there is a substantial discrepancy between the results. The M-integral in [41] is calculated for "two-dimensional elastic line singularities in the framework of plane strain linear elasticity" and for line forces it is constant in all cases. That means on the one hand, it cannot capture the interaction between the two line forces and on the other hand that it is an unphysical result, since as it has been shown in this paper the M-integral represents the elastic interaction energy between the two line forces. Moreover, the M-integral, in the last two cases in Table 2 in [41] for line forces with perpendicular strengths and with parallel strengths in y-direction, is given by the exact same formula, which is a rather strange result. A crucial point that has to be emphasized as we have already mentioned in the Introduction is that the Mintegral in [41] is calculated taking the space dimension equal to two (n = 2) and not three as it should be, since a line force is a line singularity in a three-dimensional medium and plane strain is a two-dimensional deformation in a three-dimensional body, therefore the space dimension should be taken equal to three. As it can be seen here in Table 2, the M-integral for parallel line forces with parallel/antiparallel strengths in x-or y-direction which has been calculated taking the space dimension equal to three (n = 3), is a logarithmic function depending on the distance a, describing in this way the interaction of the two line forces. • Concerning the L 3 -integral, the result in Table 2 is in agreement with the result given in [41], only if we consider that there is a typing mistake (wrong placement of the parenthesis) in [41].
Finally, it should be noted that in Seo et al. [41] a comparison between the J (AB) -, M (AB) -and L (AB) 3 -integrals of line forces in plane strain and the J-, M-and L 3 -integrals of edge dislocations has been made. However, such a comparison could lead to misleading conclusions, since these two problems have different physics. A dislocation is a topological defect and in particular is a line defect in a three-dimensional crystal and has to be treated in the framework of three-dimensional incompatible elasticity taking into account plastic fields, whereas a line force is nothing but the source term in the force equilibrium condition and has to be treated in the framework of three-dimensional compatible elasticity. Even if there are similarities between the two problems, there are still important differences from a physical point of view. The J-, M-and L 3 -integrals of straight screw and edge dislocations have been derived in Agiasofitou and Lazar [2] in the framework of three-dimensional incompatible elasticity of dislocations and there are discrepancies with the results reported in [41] in the framework of plane strain elasticity. Comparisons of the related results are out of the aim of this paper.