1. Nonlocal Functions
Modern differential calculus was developed by I. Newton and G. Leibnitz in the 17th century. It is demonstrated in the following cubic polynomial:
To introduce the derivative, we assume that
and
have some small increments,
and
. Then:
Subtracting
in Equation (1) and dividing by
, we get:
To obtain the derivative, we need to eliminate
from the right-hand side of this equation. Naturally, we cannot put
because the resulting operation makes no sense. We also cannot neglect
because the operation becomes approximate in this case and is not acceptable in mathematics. Thus, we need to introduce the infinitesimal quantity
, which is infinitely small but not zero. Then, we arrive at the following equation for the derivative:
The other definition of the derivative that does not require the introduction of infinitesimal quantities was proposed by J. Landen in 1764 [
1]. It introduces an interval
and assumes that
. We determine the derivative of the function in Equation (1) as:
Taking
, we get:
This result coincides with Equation (2) but does not require the introduction of infinitesimal quantities.
We further use the Landen definition of the derivative to construct nonlocal derivatives.
To introduce the nonlocal function, we consider a conventional function
as shown in
Figure 1.
We assume that this function describes a real physical process and, according to the definition, is smooth and has traditional derivatives of any order. We introduce in the vicinity of point
the local coordinate
, such that
(
Figure 1). Since function
is smooth, we decompose it in the vicinity of point
into the Taylor series:
where
, and the series is restricted to the terms presented in Equation (3). We introduce the nonlocal function
as the average value of function
on the interval
, i.e.,:
Substituting Equation (3), we arrive at [
2]:
Thus, the nonlocal function depends not only on the value of the original function at a point but on the value of the second derivative as well. According to Landen, we should introduce the nonlocal derivative as:
In contrast to the original Landen definition, the parameter
is small but finite. Substituting Equation (3) and taking into account Equation (4), we get:
Thus, the nonlocal derivative of the original function is the classical derivative of the nonlocal function.
In the following sections, the proposed approach is demonstrated on the classical problems of mathematical physics—the problems of a string and a membrane loaded with concentrated forces.
2. A String Loaded with a Concentrated Force
We consider a string loaded with the axial tensile force
and the transverse force
P as in
Figure 2.
The equilibrium equation has the following form:
For a relatively small deflection,
, we can take
, and Equation (6) reduces to:
To obtain the classical solution, we put
, and Equation (7) becomes:
The solution of this equation is:
Satisfying the boundary condition,
(
Figure 2), we finally get the classical solution for the string deflection:
This solution has two points of discontinuity at
and
, at which the deflection derivative does not exist. For a real string, we have the obvious conditions of
, which cannot be satisfied with Equation (9). The experimental deflection of a steel string with diameter
, length
, and elastic modulus
loaded with forces
and
is shown in
Figure 3 with dots.
The experimental facility for investigating the transverse displacements of a stretched string under the action of a transverse local force is shown in
Figure 4. The equipment provides a constant controlled string tension. During the test, the transverse displacements of the string are measured with an accuracy of 0.1 mm using an optical system.
As can be seen, the classical solution in Equation (10) presented with the dashed line does not coincide with the experiment.
We obtain the nonlocal solution of the problem. Using
Figure 1 and Equation (4), we introduce the nonlocal angle:
where angle
is averaged over some intervals
. Changing
to
in Equation (7), we arrive at the following equilibrium equation:
Since
, we get:
This equation allows us to clear out the procedure used to construct the nonlocal solution. The left-hand part contains the ordinary Helmhotz-type operator acting on the string deflection, whereas its right-hand part is the classical solution of the problem specified by Equation (9).
This conclusion is valid in the general case [
2] in which the original function
can be found from the following equation:
where
is the Laplace operator, and
is the classical solution. If
is a scalar function, the Laplace operator is invariant, and Equation (13) is valid in any coordinate frame.
Equation (11) includes parameter
, which is not known yet. To determine this parameter, we should take into account that Equation (10) describes the mathematical model of the string, which is the one-dimensional manifold. The real string is characterized with an elastic modulus and the shape and dimensions of the cross-section, which do not enter into the classical solution, Equation (10). We consider a more adequate physical model—a beam shown in
Figure 5.
In this model, the axial force
acting in the beam cross-section is supplemented with transverse force
and bending moment
(
Figure 5). A beam is described by the following equations:
which can be reduced to:
where
is the bending stiffness. Matching this equation to Equation (11), we can conclude that
. Since
and
cannot be negative, parameter
is imaginary. This result looks natural because the obtained results are based on a model of a homogeneous continuum that ignores the material’s actual microstructure. Thus, we cannot expect that this model allows us to determine the actual value of the microstructural parameter. To avoid imaginary values, we present Equation (12) in the following form:
in which:
Thus, parameter depends on the bending stiffness and the axial force and does not depend on the transverse force and the string length. Putting , we get , and Equation (14) reduces to the classical Equation (9).
The solution of Equation (14) is:
This solution specifies two boundary-layer effects in the vicinity of points
and
. For a long string, we can neglect the interaction of these effects and present Equation (16) as:
We consider the boundary condition at
, i.e.,
. Taking
and using this condition to determine
, we arrive at the solution that is valid for
:
We consider the part of the string corresponding to
. To apply the boundary conditions
, we put
and get:
Using the boundary conditions, we can find
and
and present Equations (17) and (18) as:
For
, we can neglect the exponential terms in comparison with unity, and both solutions yield one and the same result. For the experimental string, the obtained solution is shown in
Figure 3 with the solid line. As can be seen, it is in good agreement with the experiment.
3. A Circular Membrane Loaded with a Concentrated Force
A circular membrane loaded with in-plane forces
and a concentrated force
applied at the membrane center (
Figure 6) is a classical singular problem of mathematical physics [
3].
The membrane deflection satisfies the following equation [
3]:
For a concentrated force,
, in which
is the delta function. The solution of Equation (19) that satisfies the boundary condition
is [
3]:
As can be seen, if . Moreover, the deflection is infinitely high for any force irrespective of how small it can be. Naturally, this result does not have a physical meaning.
To obtain the nonlocal solution of the problem, we apply Equation (13). Taking into account that
is a scalar function, we have [
2]:
where
is specified by Equation (20). To determine parameter
, we assume, as earlier, that the membrane is a three-dimensional object—a circular plate with a finite thickness
, rather than a two-dimensional mathematical manifold. The plate deflection satisfies the following equation [
4]:
where
,
is the elastic modulus, and
is Poisson’s ratio. Matching Equations (21) and (22), we can, as earlier, conclude that
. Introducing parameter
in accordance with Equation (15), we arrive at:
Using the following notations:
we finally have:
The general solution of this equation can be written in terms of modified Bessel functions
,
and the corresponding particular solution as:
Using the boundary condition
, we can determine the constant
and get:
As can be seen, this particular solution is singular at
. However, the Macdonald function
has the same type of singularity and can be used to eliminate the singularity of the solution in Equation (24). We decompose the Bessel functions into the power series [
5]:
where
is the Euler constant. Substituting series (25) in Equation (24) and putting
, we can conclude that the solution becomes regular if we take
. Finally, we arrive at:
At the membrane center, the deflection is finite, i.e.,:
The derivative of the deflection,
is zero at the membrane center because
and
.
We undertake the analysis of the obtained results. First, the singular classical solution in Equation (20) looks consistent since the right-hand side of Equation (19) includes the delta function, which is also singular. Thus, a singular action results in a singular solution. Being natural in mathematics, it is not acceptable in physics, in which the concentrated force and the membrane deflection are, in principle, different functions. A concentrated force does not exist in reality and is singular according to its definition, whereas a deflection is a physical variable that can be directly measured and cannot be infinitely high. Second, the consistency of the classical solution is proved in mathematics since it can be presented as a limit of a system of regular functions [
6]. Thus, the mathematically consistent solution is not physically consistent. Finally, as follows from the foregoing derivation, the logarithmic-type singularity in the right-hand side of Equation (23) is compensated for by the same type of singularity in a fundamental solution of the corresponding homogeneous equation. If we formally change the type of singularity in the right-hand side of Equation (23) (e.g., introduce a power-type singularity), it will not be eliminated by the logarithmic fundamental solution. Thus, the proposed method works only for the equations that describe physical problems.
The obtained solution was verified experimentally. The experimental membrane was made of a polymeric film and had the following parameters:
,
,
,
,
H. For these parameters, Equation (15) yields
mm. The calculated value of the maximum deflection—
is in good agreement with the experimental result—
. The description of the experiment is presented in [
7].
Dependence of the deflection on the radial coordinate is shown in
Figure 7 with the solid line along with experimental results (dots). The line corresponds to both solutions (classical and nonlocal)—the difference between them can be seen only in the vicinity of the membrane center (
Figure 8).
As follows from
Figure 8, the difference between the curves is observed for
, whereas the membrane radius is
.