1. Introduction
The first problem of optimization of elastic rods was formulated by Lagrange, [
1]. The problem consists of finding the shape of an elastic rod of given volume that has largest value of the buckling force. The solution of the problem, with the simply supported boundary conditions, was obtained by Clausen in [
2]. Various versions of the optimal shape of a column problem were treated in a number of publications, see [
3,
4,
5,
6,
7,
8,
9,
10].
In this work we propose to solve the problem of the strongest column in a constant gravity field for the case when the upper end is fixed, as in [
10], and, additionally, loaded by a constant concentrated force at the top. Thus, we shall be able to reproduce both classical Clausen problem [
2] as well as the problem of heavy inverted column, treated in [
10]. Also we shall examine the bifurcation points of nonlinear equilibrium equations. Namely, we shall show that the lowest bifurcation point of linear and nonlinear problem coincide. This is important, since optimization will be performed at eigenvalues of the linearized problem. Also for differential equations describing the optimal shape of the column, a variational principle is formulated and two new first integrals are obtained. We also studied the invariance of the variational principle, by using the Noether’s theorem, see [
11,
12,
13], where integer and fractional order systems are treated. We showed that one of two first integrals follow from the Absolute invariance of the Hamilton’s action integral and the other from so called Gauge invariance. Pontryagin’s principle and variational methods, including Noether’s theory, represent powerful methods for the study of mechanical and physical systems. Our results demonstrate this on an example where new results are obtained by using these methods. For recent contribution to a nontrivial extension of the Noether’s theory see [
14]. In this work we will follow the notation of [
11,
15].
2. Formulation
In
Figure 1 we show the column with inextensible axis of length
At upper end
B the column is fixed with the possibility of sliding along the axis
A concentrated force of intensity
F is applied at the end
B. Equations describing behaviour of the column [
15] are
Here
H and
V are components of cross-sectional force along
x and
y axes,
M denotes the bending moment,
denotes the angle between the tangent to the column axis and the
x axis of a coordinate system
Also
S denotes the arc-length of the column axis. We use the specific distributed force given as
with
being the mass density,
g is the gravitational constant and
is area of the rod at an specific point of the rod axis. In order to apply the Bernoulli–Euler bending theory we assume that
Also
express inextensibility of the axis. Constitutive equation of the classical Bernoulli–Euler theory reads
In (
2) and (
3)
x and
y denote coordinates of an arbitrary point in system
. We use
E to denote modulus of elasticity and
I the second moment of the cross-section area. From
Figure 1 we conclude that
so that
The column has volume given by
If the column cross-sections are similar and similarly oriented, then
with
being a constant, equal
for circular cross-section. We introduce
so that (
1)–(
8) become
with
The volume of the column now becomes
For (
9) and (
10) we find a trivial solution
valid for all values of
and
and for any
We examine the values of
that lead to a nontrivial solution of (
9) and (
10). First we write (
9)
in operator form as
with
Then (
9)
subject to (
10)
is equivalent to
The solution
is trivial solution valid for all
The Fréchet derivative of
calculated at
is
We consider the linearized boundary value problem
, i.e.,
It is known that the necessary condition for the existence of nontrivial solution of (
9)
, (
10)
is that there is nontrivial solution of the linearized equation
subject to
Suppose that
with
fixed (gravitational force is not subject to changes) and
and
is a solution to (
12) and (
13). Then with
as a bifurcation parameter, we have:
Proposition 1. Giventhe boundary value problem (12) and (13) has only real eigenvalues, there are an infinite but countable number of them,,, and they can be ordered to satisfy The number of zeros of the eigenfunctionin the interval (0,1) isAlso, Proof of Proposition 1. Note that with
the conditions of the Theorem 4.3.1 of [
16] are satisfied. The result follows from the application of this Theorem. ☐
The condition that guarantees that (
9) and (
10) have nontrivial solution is formulated next. Our interest is to show that at the lowest eigenvalue of (
12) and (
13), that is for
, the system (
9) and (
10) with
arbitrary ,
, has a bifurcation point. We note that in the next Section
will be determined from the optimization procedure. We state this as:
Theorem 1. Letbe the lowest eigenvalue of the system (12) and (13). The nonlinear boundary value problem (11) has a bifurcation point at, i.e, there is continuously differentiable curve throughsuch thatand with Proof ofTheorem 1. We use the Crandall–Rabinowitz theorem, see [
17,
18], p.15. Thus, let
be eigenvector of (
12) and (
13) with given
and lowest eigenvalue
, i.e.,
We assume that
is normalized, so that
Note that
subject to
leads to unique eigenvector
. To prove that
is unique, observe that
If
from (
3) we conclude that
Then, the equilibrium equations for the rod lead to the conclusion that there is only trivial solution
Since
, the Theorem 5, p.73 of [
19], implies that
is unique. Next we determine
as
Therefore it follows that
Range
Thus, Theorem I.5.1 of [
18] applies and (
14) and (
15) follow. ☐
Remark 1. In principle the same results may be obtained if we treat bothandas bifurcation parameters. In that case a generalization of the Crandall–Rabinowitz theorem given in [18], p.161, must be used. However, for our purposes the result presented here suffices since it shows that for any,, the bifurcation point of (12) and (13) leads to the bifurcation of (11). The optimization problem is stated as: given
find
in (
12) and (
13) so that a nontrivial solution exists and
is
minimal. This will represent the strongest compressed inverted column. We proceed with determining optimal
3. Minimization of w for Given Load Parameters and
We rewrite the system (
12) and (
13) as
so that (
12) and (
13) transforms to
and
The optimization problem now becomes: determine the control
so that
with differential constraints (
17) and (
18).
We take U as a set of continuous nonnegative functions, defined on the interval , i.e.,
To solve the optimization problem we use the Pontryagin’s principle. For recent example of application of Pontryagin’s principle, see an example from biology [
20]. Here the Pontryagin’s function
becomes (see [
12])
where variables
, are determined from
and
From the condition
we obtain
or
We note that state variables
, and co-state variables
, have an important symmetry. Namely, by comparing
we conclude that for any
the solution
determines the co-state variables
by
This type of identification was used in [
9] and is applicable to both single and bi modal optimization of elastic rods. Since
, we conclude from (
23) that
We set
so that
From (
22) and (
24) follows
Note that (
26) is a necessary condition for
. We comment now on the sufficient conditions for the minimality of (
19). There are two main sufficiency theorems in optimal control theory. First one is the Mangasarian’s sufficiency theorem, see [
21], which requires that the objective function and constraints are convex jointly in state and control variables for the minimization problems. The second condition, known as the Arrow’s theorem, is applied as follows, see [
22]. The control variable (
23) is substituted in Pontryagin’s function
to obtain
Now, the condition
is guaranteed if in expression (
27) the function
is convex function with respect to
, when
, are fixed and positive. This is not the case here.
There are many other approaches to the problem of specifying the sufficient conditions for minimum of
For example, in [
23] the sufficient conditions involve Legendre–Clebsch condition and solution of an additional Riccati equation. Such a study is beyond the scope of our paper. Therefore, in the analysis that follow we will use only the necessary conditions for
given by (
22).
Solving (
20)
and (
21)
for
we obtain
so that (
25) leads to
Now, differentiating (
29) it follows
Thus, the optimal compressed inverted column is determined from
subject to
Our condition (
31) and (
32) reduces to condition presented in [
24]. Also, the optimality conditions presented in [
25,
26] are equivalent to our conditions. We note that in [
24,
25,
26] the results are obtained by methods different from ours. From (
29) and (
25) we obtain
5. Invariance Properties of the Integral (38)
Consider transformation of independent and dependent variables in (
38) given by
where
and
f are generators of the infinitesimal transformation group and
The case of multi-time transformation and corresponding version for Noether-type first integrals is presented in [
14]. Here, it is assumed that
and
f are continuously differentiable with respect to all variables. If action integral (
38) is invariant under the transformation (
47), then Noether’s theorem guarantees the existence of a first integral to the Euler–Lagrange’s system of equations (
34), see [
12], p. 137. Using Noether’s theorem for the present case, we state:
Theorem 3. If the generators of the infinitesimal transformation group,, f satisfywhereandis an arbitrary function continuously differentiable with respect to all variables (gauge function), then the system (34) has a first integral of the form We apply now the Theorem 3 to the Lagrangian (
41). Thus, we calculate
and substitute the result into (
48). The invariance condition then becomes
Consider two special cases generators of the infinitesimal transformation group and gauge function P:
Case 1: Suppose that
The condition (
50) is satisfied and the first integral (
49) is
given as (
40)
Note that this first integral (
51) may be writen as
Case 2: Suppose that
We shall determine
P so that (
50) is satisfied. This leads to the condition
or
By using (
46) to eliminate
we obtain
Finally, from (
29) we have
, so that substituting in the previous equation and integrating, we obtain
where
is a constant. The first integral (
49) now becomes
which is equal to (
40)
when we specify the value of constant.
6. Results of Numerical Solution to (36), (37)
Using variables (
16),
and
the system (
17), (
18), (
25) and (
28) transforms to
with boundary conditions
Thus, we have to choose
in order to satisfy
Conservation laws (
40) with
, become
In solving (
52), (
53) we used (
54) to monitor the accuracy of the integration. We solved (
52) and (
53) with
as a given parameter. The eigenvalue
is determined so that
, that is the volume of the column is given as
. In
Table 1 the results of computation are presented.
For values of
shown in
Table 1 the first integrals (
54) are constants up to the order of
. In numerical solution of (
52), (
53) and evaluation of the first integrals we used the computer package Mathcad 14.
In
Figure 2 the optimal cross-sectional area is shown that corresponds to the following values of parameter
The corresponding values of
are taken from
Table 1. We note that to
we have
This is a special case of optimally shaped light column, see [
15] p.216, and the exact value is
. Our numerical result given in the
Table 1 agrees with the exact value up to the order
. Also for the light optimal rod, the analytical result is
, while our numerical value is
Another special case presented in
Figure 2 corresponds to the case
and the compressive force is
. This shape is shown by curve
in
Figure 2.
We present also the shape of the optimally designed column in the post-critical state. Thus, we solve the following system
where
is given by (
23). The solution for
of (
55) is shown in
Figure 3 for
with
Also, in the same
Figure 3 we show the shape of the rod
according to linear theory.