Next Article in Journal
On the Stability with Respect to H-Manifolds for Cohen–Grossberg-Type Bidirectional Associative Memory Neural Networks with Variable Impulsive Perturbations and Time-Varying Delays
Previous Article in Journal
Improved Compact Cuckoo Search Algorithm Applied to Location of Drone Logistics Hub
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Shape and First Integrals for Inverted Compressed Column

by
Enes Kacapor
1,
Teodor M. Atanackovic
2,* and
Cemal Dolicanin
1
1
Department of Mathematical Sciences, State University of Novi Pazar, Vuka Karadzica 9, 36300 Novi Pazar, Serbia
2
Faculty of Technical Sciences, University of Novi Sad, Trg D. Obradovica 6, 21000 Novi Sad, Serbia
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(3), 334; https://doi.org/10.3390/math8030334
Submission received: 5 February 2020 / Revised: 12 February 2020 / Accepted: 12 February 2020 / Published: 3 March 2020

Abstract

:
We study optimal shape of an inverted elastic column with concentrated force at the end and in the gravitational field. We generalize earlier results on this problem in two directions. First we prove a theorem on the bifurcation of nonlinear equilibrium equations for arbitrary cross-section column. Secondly we determine the cross-sectional area for the compressed column in the optimal way. Variational principle is constructed for the equations determining the optimal shape and two new first integrals are constructed that are used to check numerical integration. Next, we apply the Noether’s theorem and determine transformation groups that leave variational principle Gauge invariant. The classical Lagrange problem follows as a special case. Several numerical examples are presented.

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 L . At upper end B the column is fixed with the possibility of sliding along the axis x . A concentrated force of intensity F is applied at the end B. Equations describing behaviour of the column [15] are
d H d S = q 0 A ; d V d S = 0 ; d M d S = V cos ϑ + H sin ϑ .
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 x O y . Also S denotes the arc-length of the column axis. We use the specific distributed force given as q 0 = ρ g A with ρ being the mass density, g is the gravitational constant and A > 0 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 A C [ 0 , L ] . Also
d x d S = cos ϑ ; d y d S = sin ϑ ,
express inextensibility of the axis. Constitutive equation of the classical Bernoulli–Euler theory reads
M = E I d ϑ d S .
In (2) and (3) x and y denote coordinates of an arbitrary point in system x O y . We use E to denote modulus of elasticity and I the second moment of the cross-section area. From Figure 1 we conclude that
y ( 0 ) = x ( 0 ) = 0 ; M ( 0 ) = 0 ; V ( 0 ) = 0 ; ϑ ( L ) = 0 ; H ( L ) = F ,
so that
H ( S ) = [ S L q 0 A ( ξ ) d ξ + F ] ; V ( S ) = 0 .
The column has volume given by
W = 0 L A ( S ) d S .
If the column cross-sections are similar and similarly oriented, then
I ( S ) = α A 2 ( S ) ,
with α being a constant, equal α = 1 4 π for circular cross-section. We introduce
t = S L ; a = A L 2 ; λ 1 = q 0 L α E ; λ 2 = F α E L 2 ; u = x L ; v = y L ; w = W L 3 ,
so that (1)–(8) become
( a 2 ϑ ˙ ) · + [ λ 1 t 1 a ( ξ ) d ξ + λ 2 ] sin ϑ = 0 ; u ˙ = cos ϑ ; v ˙ = sin ϑ , d ( · ) / d t = ( · ) · ,
with
u ( 0 ) = 0 ; v ( 0 ) = 0 ; lim t 0 a 2 ( t ) ϑ ˙ ( t ) = 0 ; ϑ ( 1 ) = 0 .
The volume of the column now becomes
w = 0 1 a ( t ) d t .
For (9) and (10) we find a trivial solution
ϑ 0 = 0 ; u 0 = t ; v 0 = 0 ,
valid for all values of λ 1 and λ 2 and for any a > 0 . We examine the values of ( λ 1 , λ 2 ) = K R 2 that lead to a nontrivial solution of (9) and (10). First we write (9) 1 in operator form as
M ( K , ϑ ) = ( a 2 ϑ ˙ ) · + [ λ 1 t 1 a ( ξ ) d ξ + λ 2 ] sin ϑ , ϑ U ,
with
U = { ϑ : ϑ C 2 ( 0 , 1 ) ; lim t 0 a 2 ( t ) ϑ ˙ ( t ) = 0 ; ϑ ( 1 ) = 0 } .
Then (9) 1 subject to (10) 3 , 4 is equivalent to
M ( K , ϑ ) = 0 .
The solution ϑ 0 = 0 is trivial solution valid for all K R 2 . The Fréchet derivative of M ( K , ϑ ) calculated at ϑ 0 = 0 is
D M ϑ ( K , 0 ) ϑ = B ( K ) ϑ = ( a 2 ϑ ˙ ) · + [ λ 1 t 1 a ( ξ ) d ξ + λ 2 ] ϑ .
We consider the linearized boundary value problem B ( K ) ϑ = 0 , i.e.,
( a 2 ϑ ˙ ) · + [ λ 1 t 1 a ( ξ ) d ξ + λ 2 ] ϑ = 0 .
It is known that the necessary condition for the existence of nontrivial solution of (9) 1 , (10) 3 , 4 is that there is nontrivial solution of the linearized equation
( a 2 ϑ ˙ ) · + [ λ 1 t 1 a ( ξ ) d ξ + λ 2 ] ϑ = 0 ,
subject to
lim t 0 a 2 ( t ) ϑ ˙ ( t ) = 0 ; ϑ ( 1 ) = 0 .
Suppose that K = ( λ , λ 2 ) with λ 1 = λ > 0 fixed (gravitational force is not subject to changes) and λ 2 R and ϑ 0 is a solution to (12) and (13). Then with λ 2 as a bifurcation parameter, we have:
Proposition 1.
Given λ 1 = λ the boundary value problem (12) and (13) has only real eigenvalues, there are an infinite but countable number of them, λ 2 n , n = 1 , 2 , , and they can be ordered to satisfy
0 < λ 21 < λ 22 < λ 23 <
The number of zeros of the eigenfunction ϑ 2 n in the interval (0,1) is n 1 . Also,
λ 2 n n 2 π 2 ( 0 1 d t a ( t ) ) 2 as n .
Proof of Proposition 1.
Note that with a 2 = p , λ 1 t 1 a ( ξ ) d ξ = q , w = 1 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 K = ( λ , λ 21 ) , the system (9) and (10) with arbitrary a ( t ) > 0 , t ( 0 , 1 ) , has a bifurcation point. We note that in the next Section a ( t ) will be determined from the optimization procedure. We state this as:
Theorem 1.
Let λ 21 = Λ be the lowest eigenvalue of the system (12) and (13). The nonlinear boundary value problem (11) has a bifurcation point at K = ( λ , Λ ) , i.e, there is continuously differentiable curve through ( ( λ , Λ ) , 0 ) such that
{ ( ( λ , Λ ( s ) ) , ϑ ( s ) ) , s ( ε , ε ) , ( ( λ , Λ ( 0 ) ) , ϑ ( 0 ) ) = ( ( λ , Λ ) , 0 ) } ,
and with
M ( K = ( λ , Λ ( s ) ) , ϑ ( s ) ) = 0 , s ( ε , ε ) .
Proof ofTheorem 1.
We use the Crandall–Rabinowitz theorem, see [17,18], p.15. Thus, let ϑ 0 be eigenvector of (12) and (13) with given λ and lowest eigenvalue Λ , i.e., K ( λ , Λ ) . We assume that ϑ 0 is normalized, so that 0 1 | ϑ 0 ( t ) | 2 d t = 1 . Note that
( a 2 ϑ ˙ ) · + [ λ t 1 a ( ξ ) d ξ + Λ ] ϑ = 0 ,
subject to
lim t 0 a 2 ( t ) ϑ ˙ ( t ) = 0 ; ϑ ( 1 ) = 0 ,
leads to unique eigenvector ϑ = ϑ 0 . To prove that ϑ 0 is unique, observe that a ( 1 ) 0 . If a ( 1 ) = 0 , from (3) we conclude that M ( 1 ) = 0 . Then, the equilibrium equations for the rod lead to the conclusion that there is only trivial solution ϑ 0 = 0 . Since a ( 1 ) 0 , the Theorem 5, p.73 of [19], implies that ϑ 0 is unique. Next we determine B ( K ) ϑ 0 Λ as
D 2 M ϑ , Λ ( K , 0 ) ϑ 0 = ϑ 0 .
Therefore it follows that D 2 M ϑ , Λ ( K , 0 ) ϑ 0 Range B ( K = ( λ , Λ ) ) . 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 both λ 1 and λ 2 as 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 a ( t ) > 0 , t ( 0 , 1 ) , the bifurcation point of (12) and (13) leads to the bifurcation of (11).
The optimization problem is stated as: given K 0 = ( λ , Λ ) find a * ( t ) in (12) and (13) so that a nontrivial solution exists and
w = 0 1 a * ( t ) d t
is minimal. This will represent the strongest compressed inverted column. We proceed with determining optimal a * ( t ) .

3. Minimization of w for Given Load Parameters λ and Λ

We rewrite the system (12) and (13) as
ϑ = x 1 ; a 2 ϑ ˙ = x 2 ; x 3 = t 1 a ( ξ ) d ξ ,
so that (12) and (13) transforms to
x ˙ 1 = x 2 a 2 ; x ˙ 2 = λ x 1 x 3 Λ x 1 ; x ˙ 3 = a ,
and
x 1 ( 1 ) = 0 ; x 2 ( 0 ) = 0 ; x 3 ( 1 ) = 0 .
The optimization problem now becomes: determine the control a * ( t ) U so that
min a U I = min a U 0 1 a ( t ) d t = 0 1 a * ( t ) d t ,
with differential constraints (17) and (18).
We take U as a set of continuous nonnegative functions, defined on the interval [ 0 , 1 ] , i.e., U = { u : u C 1 ( [ 0 , 1 ] , R + ) ; u ( t ) 0 } .
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 H becomes (see [12])
H = a + p 1 x 2 a 2 + p 2 ( λ x 1 x 3 Λ x 1 ) + p 3 ( a ) ,
where variables p i , i = 1 , 2 , 3 , are determined from
p ˙ 1 = H x 1 = p 2 ( λ x 3 + Λ ) ; p ˙ 2 = H x 2 = p 1 a 2 ; p ˙ 3 = H x 3 = λ x 1 p 2 ,
and
p 1 ( 0 ) = 0 ; p 2 ( 1 ) = 0 ; p 3 ( 0 ) = 0 .
From the condition min a U H we obtain
H a = 1 2 p 1 x 2 a 3 p 3 = 0 ,
or
a = ( 2 p 1 x 2 1 p 3 ) 1 / 3 .
We note that state variables x i , i = 1 , 2 , and co-state variables p i , i = 1 , 2 , have an important symmetry. Namely, by comparing
x ˙ 1 = x 2 a 2 ; x ˙ 2 = x 1 ( λ x 3 + Λ ) ; x 1 ( 1 ) = 0 ; x 2 ( 0 ) = 0 ; p ˙ 1 = p 2 ( λ x 3 + Λ ) ; p ˙ 2 = p 1 a 2 ; p 1 ( 0 ) = 0 ; p 2 ( 1 ) = 0 ,
we conclude that for any c 0 the solution x 1 , x 2 determines the co-state variables p 1 , p 2 by
p 1 = c x 2 ; p 2 = c x 1 .
This type of identification was used in [9] and is applicable to both single and bi modal optimization of elastic rods. Since a 0 , we conclude from (23) that c > 0 . We set c = 1 so that
p 1 = x 2 ; p 2 = x 1 .
Thus, (23) becomes
a = ( 2 x 2 2 1 p 3 ) 1 / 3 .
From (22) and (24) follows
2 H a 2 = 6 x 2 2 a 4 0 .
Note that (26) is a necessary condition for min a U H . 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 H to obtain
H = 2 p 1 x 2 ( 2 x 2 2 1 p 3 ) 2 / 3 + p 2 ( λ x 1 x 3 Λ x 1 ) .
Now, the condition min a U H is guaranteed if in expression (27) the function H is convex function with respect to x i , i = 1 , 2 , 3 , when p i , i = 1 , 2 , 3 , 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 H . 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 min a U H given by (22).
Solving (20) 3 and (21) 3 for p 3 we obtain
p 3 = λ 0 t x 1 2 ( ξ ) d ξ ,
so that (25) leads to
1 + λ 0 t ϑ 2 ( ξ ) d ξ = 2 a ϑ ˙ 2 .
Now, differentiating (29) it follows
( a ϑ ˙ 2 ) · λ 2 ϑ 2 = 0 .
Thus, the optimal compressed inverted column is determined from
( a 2 ϑ ˙ ) · + λ ϑ t 1 a ( ξ ) d ξ + Λ ϑ = 0 ; ( a ϑ ˙ 2 ) · λ 2 ϑ 2 = 0 ,
subject to
lim t 0 a 2 ( t ) ϑ ˙ ( t ) = 0 ; ϑ ( 1 ) = 0 .
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
a ( 0 ) ϑ ˙ 2 ( 0 ) = 1 2 ; a ( 0 ) = 0 .

4. Variational Principle and First Integrals for (31), (32)

Let
z ˙ = a ( t ) ,
with z ( 1 ) = 0 . Then z ( t ) is positive decreasing function. With (33) the system (31) becomes
( z ˙ 2 ϑ ˙ ) · + λ ϑ z + Λ ϑ = 0 ; ( z ˙ ϑ ˙ 2 ) · + λ 2 ϑ 2 = 0 ,
subject to
lim t 0 z ˙ 2 ( t ) ϑ ˙ ( t ) = 0 ; ϑ ( 1 ) = 0 ; w * = z ( 0 ) ; z ( 1 ) = 0 .
Without loss of generality, we analyze (34) and (35) for w * = 1 . With w * = 1 we have
( z ˙ 2 ϑ ˙ ) · + λ ϑ z + Λ ϑ = 0 ; ( z ˙ ϑ ˙ 2 ) · + λ 2 ϑ 2 = 0 ,
and
lim t 0 z ˙ 2 ( t ) ϑ ˙ ( t ) = 0 ; ϑ ( 1 ) = 0 ; z ( 0 ) = 1 ; z ( 1 ) = 0 .
We state now our central result as:
Theorem 2.
For the solution to (36) and (37) we have:
1. 
( z , ϑ ) gives stationary value
δ J ( z , ϑ , Λ ) = 0 ,
to the functional
J ( Z , Θ , Λ ) = 1 2 0 1 ( Z ˙ 2 Θ ˙ 2 λ Θ 2 Z Λ Θ 2 ) d t .
Here the functions ( Z , Θ ) belong to the space K defined as
K = { k : k = ( Z , Θ ) ; Z ( 0 ) = 1 ; Z ( 1 ) = 0 ; lim t 0 Z ˙ 2 ( t ) Θ ˙ ( t ) = 0 ; Θ ( 1 ) = 0 } .
2. 
The functional (38) on the solution of (36) and (37) has value equal to zero, i.e.,
J ( z , ϑ , Λ ) = 1 2 0 1 ( z ˙ 2 ϑ ˙ 2 λ ϑ 2 z Λ ϑ 2 ) d t = 0 .
3. 
The following first integrals of (36), (37) hold
3 2 z ˙ 2 ϑ ˙ 2 + 1 2 λ ϑ 2 z + 1 2 Λ ϑ 2 = C , 8 z ˙ ϑ ˙ 2 z 5 z ˙ 2 ϑ ˙ ϑ + 4 + 4 Λ λ [ 1 + 2 z ˙ ϑ ˙ 2 ] = 2 C t ,
where C = 3 2 z ˙ ( 1 ) 2 ϑ ˙ ( 1 ) 2 and is given as
C = 1 6 [ ( 4 Λ λ ) 2 + 12 ( 1 + Λ λ ) 1 ϑ ˙ ( 1 ) 2 + 4 Λ λ ] 2 ϑ ˙ ( 1 ) 2 .
Proof of Theorem 2.
Let
L = 1 2 ( Z ˙ 2 Θ ˙ 2 λ Θ 2 Z Λ Θ 2 ) .
Let I ( Z , Θ ) be a functional defined as
I ( Z , Θ ) = 0 1 L d t ,
and consider the problem of minimizing I when ( Z , Θ ) K . The necessary condition for optimality are the Euler–Lagrange equations
d d t ( L Z ˙ ) L Z = 0 ; d d t ( L Θ ˙ ) L Θ = 0 ,
which are equal to (36). Therefore
δ J ( z , ϑ , Λ ) = 0 ,
which is 1 in Theorem 1. Next by multiplying (36) 1 by ϑ and integrating, we get
0 1 ( z ˙ 2 ϑ ˙ 2 + λ ϑ 2 z + Λ ϑ 2 ) d t = 0 ,
that proves (39). Since L given by (41) does not have explicit dependance on t, we conclude that Hamiltonian H is a constant (i.e., H = L z ˙ z ˙ + L ϑ ˙ ϑ ˙ L = c o n s t . ). Thus, for (36) and (37) we have [12]
L z ˙ z ˙ + L ϑ ˙ ϑ ˙ L = 3 2 z ˙ 2 ϑ ˙ 2 + λ 2 ϑ 2 z + Λ 2 ϑ 2 = c o n s t .
or
3 2 z ˙ 2 ϑ ˙ 2 + λ 2 ϑ 2 z + Λ 2 ϑ 2 = 3 2 z ˙ 2 ( 1 ) ϑ ˙ 2 ( 1 ) .
To determine z ˙ in terms of ϑ ˙ we multiply ( 36 ) 1 by ϑ and ( 36 ) 2 by 2 z , add the result, integrate and use the boundary conditions to obtain
0 1 z ˙ 2 ϑ ˙ 2 d t = 1 Λ 0 1 ϑ 2 d t .
Also from (29) we have
0 1 ϑ 2 d t = 1 λ [ 1 + 2 z ˙ ( 1 ) ϑ ˙ ( 1 ) 2 ] .
By combining (43) and (44) and integrating (42) we obtain
3 2 z ˙ ( 1 ) 2 ϑ ˙ ( 1 ) 2 4 Λ λ z ˙ ( 1 ) ϑ ˙ ( 1 ) 2 2 ( 1 + Λ λ ) = 0 .
By solving (45) for z ˙ ( 1 ) and by using the result in (42) we obtain (40) 1 . Next we multiply (36) 1 by ϑ and (36) 2 by z and substitute in the resulting equations the term z ˙ 2 ϑ ˙ 2 from (42) to obtain
( z ˙ 2 ϑ ˙ ϑ ) · = z ˙ ( 1 ) 2 ϑ ˙ ( 1 ) 2 4 3 [ λ ϑ 2 z + Λ ϑ 2 ] , ( z z ˙ ϑ ˙ 2 ) · = z ˙ ( 1 ) 2 ϑ ˙ ( 1 ) 2 5 6 λ ϑ 2 z 1 3 Λ ϑ 2 .
Using (42) and eliminating the term λ ϑ 2 z from (46) we obtain
8 ( z z ˙ ϑ ˙ 2 ) · 5 ( z ˙ 2 ϑ ˙ ϑ ) · = 3 z ˙ ( 1 ) 2 ϑ ˙ ( 1 ) 2 + 4 Λ ϑ 2 .
Integrating the previous relation and using the boundary conditions it follows
8 z z ˙ ϑ ˙ 2 5 z ˙ 2 ϑ ˙ ϑ = 3 z ˙ ( 1 ) 2 ϑ ˙ ( 1 ) 2 t 4 + 4 Λ 0 t ϑ 2 d ξ .
Finally, by expressing 0 t ϑ 2 d ξ from (29) and recalling that 3 z ˙ ( 1 ) 2 ϑ ˙ ( 1 ) 2 = 2 C , we get
8 z ˙ ϑ ˙ 2 z 5 z ˙ 2 ϑ ˙ ϑ = 2 C t 4 4 Λ λ [ 1 + 2 z ˙ ϑ ˙ 2 ] ,
so that (40) 2 holds. ☐
Remark 2.
In the special case Λ = 0 we have C = 2 , so that first integrals take the form
3 2 z ˙ 2 ϑ ˙ 2 + 1 2 λ ϑ 2 z = 2 , 5 z ˙ 2 ϑ ˙ ϑ 8 z ˙ ϑ ˙ 2 z = 4 ( 1 t ) ,
obtained in [10]. Also, by evaluating (40) at t = 0 , it follows ϑ ( 0 ) = 2 λ .

5. Invariance Properties of the Integral (38)

Consider transformation of independent and dependent variables in (38) given by
z ¯ ( t ¯ ) = z ( t ) + ε F z ( t , z ( t ) , ϑ ( t ) ) , ϑ ¯ ( t ¯ ) = ϑ ( t ) + ε F ϑ ( t , z ( t ) , ϑ ( t ) ) , t ¯ = t + ε f ( t , z ( t ) , ϑ ( t ) ) ,
where F z , F ϑ and f are generators of the infinitesimal transformation group and ε 1 . The case of multi-time transformation and corresponding version for Noether-type first integrals is presented in [14]. Here, it is assumed that F z , F ϑ 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 z , F ϑ , f satisfy
L z F z + L ϑ F ϑ + L z ˙ ( F ˙ z z ˙ f ˙ ) + L ϑ ˙ ( F ˙ ϑ ϑ ˙ f ˙ ) + L t f + L f ˙ P ˙ = 0 ,
where L = 1 2 ( z ˙ 2 ϑ ˙ 2 λ ϑ 2 z Λ ϑ 2 ) and P = P ( t , z , ϑ ) is an arbitrary function continuously differentiable with respect to all variables (gauge function), then the system (34) has a first integral of the form
L z ˙ ( F z z ˙ f ) + L ϑ ˙ ( F ϑ ϑ ˙ f ) + L f P = c o n s t .
We apply now the Theorem 3 to the Lagrangian (41). Thus, we calculate L z , , L ϑ ˙ and substitute the result into (48). The invariance condition then becomes
1 2 λ ϑ 2 F z ( λ ϑ z + Λ ϑ ) F ϑ + z ˙ ϑ ˙ 2 ( F ˙ z z ˙ f ˙ ) + z ˙ 2 ϑ ˙ ( F ˙ ϑ ϑ ˙ f ˙ ) + 1 2 ( z ˙ 2 ϑ ˙ 2 λ ϑ 2 z Λ ϑ 2 ) f ˙ P ˙ = 0 .
Consider two special cases generators of the infinitesimal transformation group F z , F ϑ , f and gauge function P:
  • Case 1: Suppose that F z = F ϑ = P = 0 , f = A , A = c o n s t . The condition (50) is satisfied and the first integral (49) is
    3 2 z ˙ 2 ϑ ˙ 2 + 1 2 λ ϑ 2 z + 1 2 Λ ϑ 2 = C ,
    given as (40) 1 . Note that this first integral (51) may be writen as 2 z ˙ 2 ϑ ˙ 2 L = C .
  • Case 2: Suppose that F z = z , F ϑ = ϑ 2 , f = 0 . We shall determine P so that (50) is satisfied. This leads to the condition
    P ˙ = [ L + z ˙ 2 ϑ ˙ 2 λ 2 ϑ 2 z ] ,
    or
    P ˙ = C + 3 z ˙ 2 ϑ ˙ 2 1 2 λ ϑ 2 z = C 3 2 λ ϑ 2 z Λ ϑ 2 .
    By using (46) to eliminate λ ϑ 2 z we obtain
    P ˙ = C 4 + 9 8 ( z ˙ 2 ϑ ˙ ϑ ) · Λ 2 λ ( 1 + 2 z ˙ ϑ ˙ 2 ) · .
    Finally, from (29) we have ϑ 2 = ( 1 + 2 z ˙ ϑ ˙ 2 ) · λ , so that substituting in the previous equation and integrating, we obtain
    P = C 4 t + 9 8 z ˙ 2 ϑ ˙ ϑ Λ 2 λ ( 1 + 2 z ˙ ϑ ˙ 2 ) + C p ,
    where C p is a constant. The first integral (49) now becomes
    8 z ˙ ϑ ˙ 2 z 5 z ˙ 2 ϑ ˙ ϑ 2 C t + 4 Λ λ [ 1 + 2 z ˙ ϑ ˙ 2 ] = c o n s t .
    which is equal to (40) 2 when we specify the value of constant.

6. Results of Numerical Solution to (36), (37)

Using variables (16), x 4 = p 3 and x 5 = Λ the system (17), (18), (25) and (28) transforms to
x ˙ 1 = x 2 ( 2 x 2 2 1 x 4 ) 2 / 3 ; x ˙ 2 = λ x 1 x 3 x 5 x 1 ; x ˙ 3 = ( 2 x 2 2 1 x 4 ) 1 / 3 ; x ˙ 4 = λ x 1 2 ; x ˙ 5 = 0 ,
with boundary conditions
x 1 ( 1 ) = 0 ; x 2 ( 0 ) = 0 ; x 3 ( 1 ) = 0 ; x 4 ( 0 ) = 0 .
Thus, we have to choose Λ = x 5 in order to satisfy x 3 ( 1 ) = ( t 1 a ( ξ ) d ξ ) | t = 1 = 0 . Conservation laws (40) with x i , i = 1 , , 4 , become
3 2 ( x 2 ( t ) ) 2 ( 2 x 2 ( t ) 2 1 x 4 ( t ) ) 2 / 3 + λ 2 ( x 1 ( t ) ) 2 x 3 ( t ) + Λ 2 ( x 1 ( t ) ) 2 = C , 4 ( 1 x 4 ( t ) ) x 3 ( t ) 5 x 1 ( t ) x 2 ( t ) + 4 + 4 Λ λ x 4 ( t ) = 2 C t , t ( 0 , 1 ) , C = 1 6 [ ( 4 Λ λ ) 2 + 12 ( 1 + Λ λ ) 1 ϑ ˙ ( 1 ) 2 4 Λ λ ] 2 ϑ ˙ ( 1 ) 2 , ϑ ˙ ( 1 ) = x 2 ( 1 ) ( 2 x 2 2 ( 1 ) 1 x 4 ( 1 ) ) 2 / 3 .
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 Λ = x 5 is determined so that x 3 ( 0 ) = 1 , that is the volume of the column is given as w = 1 . 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 10 8 . 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 λ : λ 0 = 0 , λ 1 = 1 , λ 3 = 3 , λ = 3.893 The corresponding values of Λ are taken from Table 1. We note that to λ 0 = 0 we have Λ = 3.289868133696451 . This is a special case of optimally shaped light column, see [15] p.216, and the exact value is Λ = π 2 3 . Our numerical result given in the Table 1 agrees with the exact value up to the order 10 10 . Also for the light optimal rod, the analytical result is a ( 1 ) = 4 3 , while our numerical value is a ( 1 ) = 1.333435777857 . Another special case presented in Figure 2 corresponds to the case λ = 3.893026 and the compressive force is Λ = 0 . This shape is shown by curve λ 3.893 in Figure 2.
We present also the shape of the optimally designed column in the post-critical state. Thus, we solve the following system
X ˙ = cos Θ ; X ( 0 ) = 0 ; Y ˙ = sin Θ ; Y ( 0 ) = 0 ; Θ ˙ = M ( a * ) 2 ; Θ ( 1 ) = 0 ; M ˙ = [ λ t 1 a * ( ξ ) d ξ + Λ ] sin Θ ; M ( 0 ) = 0 ,
where a * is given by (23). The solution for X , Y of (55) is shown in Figure 3 for λ = λ c r , Λ = Λ c r + 0.5 , with ( λ c r = 1 , Λ c r = 2.433148 ) . Also, in the same Figure 3 we show the shape of the rod x , y according to linear theory.

7. Conclusions

We analysed the problem of finding the shape of the strongest compressed inverted column in a constant gravity field. The problem treated here, for the case when the column has constant cross-section, is analyzed in [27,28,29]. Our main results may be summarized as:
  • We proved that the eigenvalues of the linearized problem are bifurcation points of the nonlinear equilibrium equations for arbitrary cross-sectional area, Theorem 1.
  • We formulated optimality conditions in the form (31) and (32). The variational principle and two first integrals are determined.
  • We determined transformation group that generates first integrals via Noether’s theorem. Those transformation groups leave the functional (38), absolutely in Case 1, and gauge invariant in Case 2.
  • The shape of the column, as a function of arc-length, is determined for several values of parameters given in Table 1. The special case λ = 0 corresponds to the known analytical solution of Lagrange problem.
  • The post critical behaviour of the optimally shaped column is presented for one specific value of parameters.
The increase of the buckling force for optimally shaped column and the column with constant cross-section is between 12% for Λ = 0 , see [10], and 33%, for λ = 0 , see [3].

Author Contributions

Investigation, E.K.; Supervision, T.M.A. and C.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Serbian Academy of Arts and Sciences, grant number F-64.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lagrange, J.L. Sur la figure des colonnes. In Oeuvres de Lagrange V. 2; Serret, M.J.A., Ed.; Gauthier-Villars: Paris, France, 1868; pp. 125–170. [Google Scholar]
  2. Clausen, T. Über die form architektonischer säulen, Bull. cl, physico math. Acad. St. Pétersbourg 1851, 9, 369–380. [Google Scholar]
  3. Tadjbakhsh, I.; Keller, J.B. Strongest columns and isoperimetric inequalities for eigenvalues. J. Appl. Mech. (Trans. ASME) 1962, 29, 159–164. [Google Scholar] [CrossRef]
  4. Olhoff, N.; Rasmussen, S.H. On the single and bimodal optimum buckling loads of clamped columns. Int. J. Solids Struct. 1977, 13, 605–614. [Google Scholar] [CrossRef]
  5. Seyranian, A.P. Lagrange problem for optimal shape of a column. Dev. Mech. (Uspehi Mehaniki) 2003, 2, 45–96. [Google Scholar]
  6. Olhoff, N.; Seyranian, A.P. Bifurcation and post-buckling analysis of bimodal optimum columns. Int. J. Solids Struct. 2008, 45, 3967–3995. [Google Scholar] [CrossRef] [Green Version]
  7. Atanackovic, T.M.; Seyranian, A.P. Application of Pontryagin’s principle to bimodal optimization problems. Struct. Multidisc. Optim. 2008, 25, 393–410. [Google Scholar] [CrossRef]
  8. Atanackovic, T.M. Optimal shape of a column with own weight: Bi and single modal optimization. Meccanica 2006, 41, 173–196. [Google Scholar] [CrossRef]
  9. Atanackovic, T.M. On the optimal shape of a compressed rotating rod. Meccanica 2004, 39, 147–157. [Google Scholar] [CrossRef]
  10. Atanackovic, T.M. Optimal shape of a strongest inverted column. J. Comput. Appl. Math. 2007, 203, 209–218. [Google Scholar] [CrossRef] [Green Version]
  11. Vujanovic, B.D.; Jones, S.E. Variational Methods in Nonconservative Phenomena; Academic Press: Boston, MA, USA, 1989. [Google Scholar]
  12. Vujanovic, B.D.; Atanackovic, T.M. An Introduction to Modern Variational Techniques in Mechanics and Engineering; Birkhäuser: Boston, MA, USA, 2004. [Google Scholar]
  13. Atanackovic, T.M.; Konjik, S.; Pilipović, S.; Simic, S.S. Variational problems with fractional derivatives: Invariance conditions and Nöther’s theorem. Nonlinear Anal. 2009, 71, 1504–1517. [Google Scholar] [CrossRef] [Green Version]
  14. Treanţă, S. Noether-Type First Integrals Associated with Autonomous Second-Order Lagrangians. Symmetry 2019, 11, 1088. [Google Scholar] [CrossRef] [Green Version]
  15. Atanackovic, T.M. Stability Theory of Elastic Rods; World Scientific: Singapore, 1997. [Google Scholar]
  16. Zettl, A. Sturm—Liouville Theory; American Mathematical Society, Mathematical Surveys and Monographs: Providence, RI, USA, 2005; Volume 121. [Google Scholar]
  17. Crandall, M.G.; Rabinowitz, P.H. Bifurcation from simple eigenvalues. J. Funct. Anal. 1971, 8, 321–340. [Google Scholar] [CrossRef] [Green Version]
  18. Kielhöfer, H. Bifurcation Theory: An Introduction with Applications to PDEs; Springer: New York, NY, USA, 2004. [Google Scholar]
  19. Churchill, R. Fourier Series and Boundary Value Problems; McGraw-Hill: New York, NY, USA, 1969. [Google Scholar]
  20. Yusof, F.M.; Abdullah, F.A.; Ismail, A.I.M. Modeling and optimal control on the spread of hantavirus infection. Mathematics 2019, 7, 1192. [Google Scholar] [CrossRef] [Green Version]
  21. Mangasarian, O.L. Sufficient conditions for the optimal control of nonlinear systems. SIAM J. Control 1966, 4, 139–152. [Google Scholar] [CrossRef]
  22. Chiang, A.C. Elements of Dynamic Optimization; McGraw-Hill: New York, NY, USA, 1992. [Google Scholar]
  23. Malanowski, K.; Maurer, H.; Pickenhain, S. Second-order sufficient conditions for state-constrained optimal control Problems. J. Optim. Theory Appl. 2004, 123, 595–617. [Google Scholar] [CrossRef]
  24. Keller, J.B.; Niordson, F.I. The tallest column. J. Math. Mech. 1966, 16, 433–446. [Google Scholar] [CrossRef]
  25. McCarthy, C.M. An Investigation of the Optimal Design of the Tallest Unloaded Column. Ph.D. Thesis, Rice University, Houston, TX, USA, 1997. [Google Scholar]
  26. McCarthy, C.M. The tallest column—Optimality revisited. J. Comput. Appl. Math. 1999, 101, 27–37. [Google Scholar] [CrossRef] [Green Version]
  27. Willers, A. Das knicken schwerer geatänge. Z. Angew. Math. Mech. (ZAMM) 1941, 21, 43–51. [Google Scholar] [CrossRef]
  28. Wang, C.Y.; Drachman, B. Stability of a heavy column with an end load. J. Appl. Mech. (Trans. ASME) 1981, 48, 668–669. [Google Scholar] [CrossRef]
  29. Wang, C.M.; Wang, C.Y.; Reddy, J.N. Exact Solutions for Buckling of Structural Members; CRC Press: Boca Raton, FL, USA, 2005. [Google Scholar]
Figure 1. Coordinate system for the rod under self-weight and concentrated force.
Figure 1. Coordinate system for the rod under self-weight and concentrated force.
Mathematics 08 00334 g001
Figure 2. Cross-sectional area of the rod.
Figure 2. Cross-sectional area of the rod.
Mathematics 08 00334 g002
Figure 3. Post-critical shape of optimally shaped rod X, Y and x = t , y , where y ˙ = ϑ , and ϑ is the eigenfunction of (12) and (13).
Figure 3. Post-critical shape of optimally shaped rod X, Y and x = t , y , where y ˙ = ϑ , and ϑ is the eigenfunction of (12) and (13).
Mathematics 08 00334 g003
Table 1. Critical values of the load parameter.
Table 1. Critical values of the load parameter.
λ01233.53.893026
Λ 3.289868 2.433148 1.583948 0.743068 0.326056 0

Share and Cite

MDPI and ACS Style

Kacapor, E.; Atanackovic, T.M.; Dolicanin, C. Optimal Shape and First Integrals for Inverted Compressed Column. Mathematics 2020, 8, 334. https://doi.org/10.3390/math8030334

AMA Style

Kacapor E, Atanackovic TM, Dolicanin C. Optimal Shape and First Integrals for Inverted Compressed Column. Mathematics. 2020; 8(3):334. https://doi.org/10.3390/math8030334

Chicago/Turabian Style

Kacapor, Enes, Teodor M. Atanackovic, and Cemal Dolicanin. 2020. "Optimal Shape and First Integrals for Inverted Compressed Column" Mathematics 8, no. 3: 334. https://doi.org/10.3390/math8030334

APA Style

Kacapor, E., Atanackovic, T. M., & Dolicanin, C. (2020). Optimal Shape and First Integrals for Inverted Compressed Column. Mathematics, 8(3), 334. https://doi.org/10.3390/math8030334

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop