Closed Form Solution in the Buckling Optimization Problem of Twisted Shafts

: The counterpart for Euler’s buckling problem is Greenhill’s problem, which studies the forming of a loop in an elastic beam under torsion. In the context of twisted shafts, the optimal shape of the beam along its axis is searched. A priori form of the cross-section remains unknown. For the solution of the actual problem, the stability equations take into account all possible convex and simply connected shapes of the cross-section. The cross-sections are similar geometric ﬁgures related by a homothetic transformation with respect to a homothetic center on the axis of the beam and vary along its axis. The distribution of material along the length of a twisted shaft is optimized so that the beam is of the constant volume and will support the maximal moment without spatial buckling. The applications of the variational method for stability problems are illustrated in this manuscript.


Introduction
Torsion shafts are tremendously significant mechanical structural elements, which are used in machinery for the transmission of power.Shafts are used broadly in machinery and are often one of the most critical parts of the machine.They are used to transfer mechanical power from the combustion engine or electromotor to a moving component, such as wheels, pulleys, or sprockets.They are used in combination with further mechanical components, such as gears and sprockets, to provide power transmission within a mechanical system.A torsion shaft characteristically consists of a cylindrical rod or bar, but can also be made in other shapes, such as triangle, squares, or hexagonal bars.Shafts are subjected to loads, including bending moments and torsional stresses, which can cause shafts to fail if they are not designed correctly.Generally, shafts are cylindrical bars that can be stepped or have the same diameter from one end to the other.Shafts that measure only one diameter across the whole shaft are much easier and cheaper to produce.Stepped shafts are made to deal with various stress levels that differ along the shaft in a mechanical system.The primary application of the optimal solutions are the moment-transmitting shafts of the engines or of the industrial machinery.Optimally stepped shafts transmit the maximal allowable moment for the predefined mass of the member.The actual manuscript provides the estimation of mass, bucking moment, stored elastic energy and stiffness of the optimally stepped shafts.This information could help the designer to estimate properly the effects of the optimization and allows for the motivated decision making for the specific design of the shafts.
"Rod" means a straight, round stick, shaft, bar, cane, or staff.The term "rod" is frequently used for the members in torsion [1][2][3][4].In structural optimization, the terms "column" or "beam" are frequently used for the designation of the member in bending state [5][6][7].A "beam" signifies a structural element that primarily resists loads applied laterally to the beam's axis.An element designed to carry the primarily axial force would be a "strut" or "column".The term "column" is applied for the heavy, vertically compressed beams."Bar" indicates a solid metal object with a uniform (round, square, hexagonal, octagonal, or rectangular) cross-section.However, the term "bar" is hardly ever used in structural mechanics.Greenhill [8] uses the word "shafting" for the designation of the twisted elastic beam.According to [9], there are no principal difference between rods, columns, and shafts in the context of stability problems.In this manuscript, the term "shaft" will be chiefly used for the description of the straight, thin, and elastic "twisted beam".
In the actual manuscript, a thin elastic shaft with an isotropic cross-section, and twisted by couples applied at its ends only [10], will be considered.The optimization for stability of the simply supported shaft were studied numerically in [1].The shaft possessed the similarly shaped cross-sections with the varying cross-sectional area.The more advanced problem of a twisted and compressed beam was studied numerically in [6].The optimization problem consisted in finding the distribution of cross-sectional areas that assigns the largest value to the critical moment.The critical moment causes the loss of stability.The constraint on the volume of material and the constraints on the admissible thickness of the rod were satisfied in the cited article.
The problem of determining the compressed and twisted beam of maximal efficiency was studied in [3].Following the structural optimization theory for a given load, the shape of the column of minimal weight was determined numerically.
The problem of the optimal design of rods under combined compression and torsion was investigated in [11].A cross-sectional area varying along the axis of the column, which led to the maximal critical loading, was sought.The varying cross-section was approximated by a function with free parameters.Alternatively, the varying cross-section was determined numerically using Pontryagin's maximum principle.
Generally speaking, for non-conservative problems, the dynamic approach for stability investigation must be applied [12].However, Euler's static approach has proved to be valid for the considered boundary value problem with conditions of fixed axes of rotation, as observed in cases 1 and 2 ([9], §5.4,Table 5.1, Section 1.2, page 127).The twisting couple Λ retains its initial direction during buckling and the boundary conditions of fixed axes of rotation are satisfied ([9], cases 1 and 2 in Table 1.1, Section 1.2, page 5).Euler's static approach is shown to be valid also for an infinite and periodically simply supported shaft.The corresponding boundary-value problems are neither self-adjoint nor conservative in the classical sense.Nevertheless, it is a conservative system of the second kind [13].Consider a thin elastic shaft with an isotropic cross-section, twisted by couples applied at its ends only (Greenhill's problem) [2,8].According to Euler's theory, the magnitude of the critical moment of the shaft is determined by the smallest positive eigenvalue.
The optimization problem for a column, loaded solely by compression forces is studied in [14].The vector of the applied forces was parallel to the axis of the column.The critical values of buckling are equal among all competitive designs of the columns.The dimensional analysis eases the mathematical technique for the optimization problem.The moments of inertia of the rod are the powers of the cross-sectional area with the exponent α.The dimensional analysis introduces two dimensionless factors, one for the total material volume and one for the total stiffness of the columns.With the method of dimensional analysis, the solution of the nonlinear algebraic equations for the Lagrange multiplier will be superfluous.The closed-form solutions for the Sturm-Liouville and mixed type boundary conditions were derived for an arbitrary positive exponent α.The solutions are expressed in terms of the higher transcendental functions.The principal results are the closed-form solution in terms of the hypergeometric and elliptic functions, the analysis of single-and bimodal regimes, and the exact bounds for the masses of the optimal columns.The isoperimetric inequality was formulated as the strict inequality sign, because the optimal solution could not be attained for any finite design parameter.The additional restriction on the minimal area of the cross-section regularized the optimization problem.
The optimal distribution of bending flexure along the axis of the simultaneously compressed and twisted rod was found in [15].The cross-section that delivers the maximum or the minimum for the critical eigenvalue must be determined among all convex and simply connected domains.The distribution of material along the length of a twisted and compressed rod was optimized so that the rod must support the maximal moment without spatial buckling, presuming its volume remains constant among all admissible rods.The static Euler's approach is applicable for the simply supported rod (hinged), twisted by the conservative moment and axial compressing force.Notably, the solution of the optimization problem for stability of simultaneously twisted and compressed rods could be determined only for three certain values of the exponent factor α. Namely, the solution in the above cited article was determined solely for α = 1, α = 2, and α = 3.In contrast, the above-cited solution [14] of the optimization of the axially compressed column allows the closed form expressions for an arbitrary positive exponent.There arises a question as to whether the optimal design problem for the shaft allows for the solutions for an arbitrary positive or a real value of the exponent factor α.
The principal task of the actual manuscript is the closed-form analytical solutions for the optimal shape of the shaft twisted along its axis.The solutions are sought for real values of the exponent factor α .

Isoperimetric Inequality for a Twisted Shaft with an Arbitrary Convex and Simply-Connected Cross-Section
In dealing with columns, rods, shafts, etc., it is assumed that they are homogeneous, elastic, and obey Hooke's law.According to Euler's theory, the magnitude of critical moment M is determined by the smallest positive eigenvalue Λ.From the beginning, the validity of the static Euler's approach for the simply supported shaft (hinged), twisted by the conservative moment, is presented.The distribution of the material along the length of the shaft is optimized so that the shaft is of the constant volume V and will support the maximal moment without spatial buckling (Figure 1).The cross section that delivers the maximum or the minimum for the critical eigenvalue must be determined among all convex and simply connected domains.A priori form of the cross-section remains unknown.For the solution of the actual problem, the stability equations take into account all possible convex and simply connected shapes of the cross-section.Consequently, the requirement of the equality of the principal moments of inertia for the cross-section could be released.
The cross-sections are similar geometric figures related by a homothetic transformation with respect to a homothetic center on the axis of the shaft and vary along its axis.The area of some reference cross-section Ω is A 0 .The area of the cross-section with the coordinate x is A(x), −L ≤ x ≤ L. The scaling dimensionless function is positive: The function s(x) is the locally integrable function over the length of the shaft (−L ≤ x ≤ L).The volume of the shaft is equal to 2V: For an arbitrary cross-section, the moments of inertia J yy , J yz , J zz of the transverse crosssectional area are the power of s(x).The moments of inertia with respect to lines, passing through a point on the neutral axis of bending and parallel to the axes y and z, read:  The constants I yy , I yz , and I zz denote the moments of inertia of the reference crosssection Ω.The moments of inertia I yy , I yz , and I zz are the components of tensor of the second rank.For an arbitrary angle Ψ between the principal axis of inertia and y axis, the moments of inertia are I yy = I 1 cos 2 Ψ + I 2 sin 2 Ψ, I yz = (I 1 − I 2 )sinΨcosΨ, I zz = I 1 sin 2 Ψ + I 2 cos 2 Ψ. (4)   The angle Ψ remains constant along the axis of the shaft.As declared above, the shaft is twisted by conservative couples, such that the torque is constant over the axis of the beam.The bending moments, in terms of curvatures of the shaft are: M z = EI yy y + EI yz z , M y = EI yz y + EI zz z . ( The axes y and z are chosen in the direction of principal axes of inertia through the center of gravity (Ψ = 0).The curvatures of the axis k z , k y in the course of buckling are assumed to be small, such that the geometrically linear equations could be applied for the solutions.The bending moments in terms of curvatures k z , k y of the shaft in the principal axes are: At first, it will be assumed that both moments of inertia are proportional to one function j(x): The constants are positive values η 1 > 0, η 2 > 0. The quantity s α is proportional to the flexural rigidity, and the exponent factor α is positive.In the simplest of cases, the exponent factor takes the values of 1, 2, and 3 (Banichuk, 1990).The case α = 2 corresponds to a congruent change in the form of the cross-section.The optimal convex shape of the simply connected cross-section with the topological genus null was determined [16].The topological genus of a surface (or Euler characteristic) is in essence the number of its "holes".The solution grounds on the following isoperimetric equation.Of all convex domains with area A 4 arbitrary angle Ψ between the principal axis of inertia and y axis, the are emains constant along the axis of the shaft.ove, the shaft is twisted by conservative couples, such that the torque axis of the beam.The bending moments, in terms of curvatures of the d  are chosen in the direction of principal axes of inertia through the = 0).The curvatures of the axis  ,  in the course of buckling are l, such that the geometrically linear equations could be applied for the ing moments in terms of curvatures  ,  of the shaft in the principal be assumed that both moments of inertia are proportional to one func- are positive values η > 0, η > 0. is proportional to the flexural rigidity, and the exponent factor  implest of cases, the exponent factor takes the values of 1, 2, and 3 he case α = 2 corresponds to a congruent change in the form of the ptimal convex shape of the simply connected cross-section with the ull was determined [16].The topological genus of a surface (or Euler ssence the number of its "holes."The solution grounds on the followuation.Of all convex domains with area  ⊿ , the equilateral triangle of the product   : with the cross-section in the form of an equilateral triangle delivers e critical eigenvalue for all convex and simply-connected domains of onal area  ⊿ .For the cross-section in the form of the equilateral trian-.9622. r cross-section, the constant is  = (4) ≈ 0.07957.The shaft with connected cross-section delivers correspondingly the minimum for ue: es describe the situations in which the form of the transverse crosshe transformation, such that one of the geometric dimensions of the es.For the technically important case of the thin-walled tubes with the f wall () and the mean diameter of tube (), the second moments − ( − /2) =  +  ≈  for  ≪  corresponds to the adjustable wall thickness (as the design variable) diameter of the tube.The analogous dependence arises in the case of wich beams [7].The case α = 3 corresponds the inconstant mean diar shaft (as the design variable) and constant wall thickness.
, the equilateral triangle yields the maximum of the product I 1 I 2 : second rank.For an arbitrary angle Ψ between the principal axis of inertia and y axis, the moments of inertia are The angle Ψ remains constant along the axis of the shaft.
As declared above, the shaft is twisted by conservative couples, such that the torque is constant over the axis of the beam.The bending moments, in terms of curvatures of the shaft are: The axes  and  are chosen in the direction of principal axes of inertia through the center of gravity (Ψ = 0).The curvatures of the axis  ,  in the course of buckling are assumed to be small, such that the geometrically linear equations could be applied for the solutions.The bending moments in terms of curvatures  ,  of the shaft in the principal axes are: At first, it will be assumed that both moments of inertia are proportional to one function (): The constants are positive values η > 0, η > 0.
The quantity  is proportional to the flexural rigidity, and the exponent factor  is positive.In the simplest of cases, the exponent factor takes the values of 1, 2, and 3 (Banichuk, 1990).The case α = 2 corresponds to a congruent change in the form of the cross-section.The optimal convex shape of the simply connected cross-section with the topological genus null was determined [16].The topological genus of a surface (or Euler characteristic) is in essence the number of its "holes."The solution grounds on the following isoperimetric equation.Of all convex domains with area  ⊿ , the equilateral triangle yields the maximum of the product   : Thus, the shaft with the cross-section in the form of an equilateral triangle delivers the maximum for the critical eigenvalue for all convex and simply-connected domains of the same cross-sectional area  ⊿ .For the cross-section in the form of the equilateral triangle  = √3 18 ⁄ ≈ 0.9622.
For the circular cross-section, the constant is  = (4) ≈ 0.07957.The shaft with the circular, simply connected cross-section delivers correspondingly the minimum for the critical eigenvalue: Two other cases describe the situations in which the form of the transverse crosssection undergoes the transformation, such that one of the geometric dimensions of the cross-section changes.For the technically important case of the thin-walled tubes with the variable thickness of wall () and the mean diameter of tube (), the second moments of: The case α = 1 corresponds to the adjustable wall thickness (as the design variable) and constant mean diameter of the tube.The analogous dependence arises in the case of the bending of sandwich beams [7].The case α = 3 corresponds the inconstant mean diameter of the tubular shaft (as the design variable) and constant wall thickness.
second rank.For an arbitrary angle Ψ between the principal axis of inertia and y axis, the moments of inertia are The angle Ψ remains constant along the axis of the shaft.
As declared above, the shaft is twisted by conservative couples, such that the torque is constant over the axis of the beam.The bending moments, in terms of curvatures of the shaft are: The axes  and  are chosen in the direction of principal axes of inertia through the center of gravity (Ψ = 0).The curvatures of the axis  ,  in the course of buckling are assumed to be small, such that the geometrically linear equations could be applied for the solutions.The bending moments in terms of curvatures  ,  of the shaft in the principal axes are: At first, it will be assumed that both moments of inertia are proportional to one function (): The constants are positive values η > 0, η > 0.
The quantity  is proportional to the flexural rigidity, and the exponent factor  is positive.In the simplest of cases, the exponent factor takes the values of 1, 2, and 3 (Banichuk, 1990).The case α = 2 corresponds to a congruent change in the form of the cross-section.The optimal convex shape of the simply connected cross-section with the topological genus null was determined [16].The topological genus of a surface (or Euler characteristic) is in essence the number of its "holes."The solution grounds on the following isoperimetric equation.Of all convex domains with area  ⊿ , the equilateral triangle yields the maximum of the product   : Thus, the shaft with the cross-section in the form of an equilateral triangle delivers the maximum for the critical eigenvalue for all convex and simply-connected domains of the same cross-sectional area  ⊿ .For the cross-section in the form of the equilateral triangle  = √3 18 ⁄ ≈ 0.9622.
For the circular cross-section, the constant is  = (4) ≈ 0.07957.The shaft with the circular, simply connected cross-section delivers correspondingly the minimum for the critical eigenvalue: Two other cases describe the situations in which the form of the transverse crosssection undergoes the transformation, such that one of the geometric dimensions of the cross-section changes.For the technically important case of the thin-walled tubes with the variable thickness of wall () and the mean diameter of tube (), the second moments of: The case α = 1 corresponds to the adjustable wall thickness (as the design variable) and constant mean diameter of the tube.The analogous dependence arises in the case of the bending of sandwich beams [7].The case α = 3 corresponds the inconstant mean diameter of the tubular shaft (as the design variable) and constant wall thickness.
Thus, the shaft with the cross-section in the form of an equilateral triangle delivers the maximum for the critical eigenvalue for all convex and simply-connected domains of the same cross-sectional area A IEW 4 second rank.For an arbitrary angle Ψ between the principal axis of inertia and y axis, the moments of inertia are The angle Ψ remains constant along the axis of the shaft.
As declared above, the shaft is twisted by conservative couples, such that the torque is constant over the axis of the beam.The bending moments, in terms of curvatures of the shaft are: The axes  and  are chosen in the direction of principal axes of inertia through the center of gravity (Ψ = 0).The curvatures of the axis  ,  in the course of buckling are assumed to be small, such that the geometrically linear equations could be applied for the solutions.The bending moments in terms of curvatures  ,  of the shaft in the principal axes are: At first, it will be assumed that both moments of inertia are proportional to one function ():  () =  (),  () =  (), () =   .
The quantity  is proportional to the flexural rigidity, and the exponent factor  is positive.In the simplest of cases, the exponent factor takes the values of 1, 2, and 3 (Banichuk, 1990).The case α = 2 corresponds to a congruent change in the form of the cross-section.The optimal convex shape of the simply connected cross-section with the topological genus null was determined [16].The topological genus of a surface (or Euler characteristic) is in essence the number of its "holes."The solution grounds on the following isoperimetric equation.Of all convex domains with area  ⊿ , the equilateral triangle yields the maximum of the product   : Thus, the shaft with the cross-section in the form of an equilateral triangle delivers the maximum for the critical eigenvalue for all convex and simply-connected domains of the same cross-sectional area  ⊿ .For the cross-section in the form of the equilateral triangle  = √3 18 ⁄ ≈ 0.9622.
For the circular cross-section, the constant is  = (4) ≈ 0.07957.The shaft with the circular, simply connected cross-section delivers correspondingly the minimum for the critical eigenvalue: Two other cases describe the situations in which the form of the transverse crosssection undergoes the transformation, such that one of the geometric dimensions of the cross-section changes.For the technically important case of the thin-walled tubes with the variable thickness of wall () and the mean diameter of tube (), the second moments of: . For the cross-section in the form of the equilateral triangle

Optimization Problem and Isoperimetric Inequality for Stability
Based on the above solution, the distribution of material along the length of a twiste beam will be optimized.The optimal distribution of the cross-sectional area along the spa of the beam is denoted by the capital letter ().The optimization problem examines th shaft of the constant volume 2 that supports the maximal moment without spati buckling.If the density of the material is , the mass 2 of the shaft is given by:  = ,   = ()  () > 0. ( The new variables will be introduced: With the new variables, the formulation of the optimization problem is the followin Particularly, if both principal moments of the second order are equal  () =  ( the factors are:  =  = 1.For example, this setting is valid for the shaft with the cros section in form of an equilateral triangle or for the shaft with a circular cross-section. The distribution of material along the length of a shaft is optimized so that the beam is of constant volume and provides the maximal moment without spatial buckling.Th variational methods will be applied for the derivation of the necessary optimality cond tions.The augmented Lagrangian is related to, but not identical to the method of Lagrang multipliers [10].Augmented Lagrangian methods are a certain class of algorithms or an alytic techniques for solving constrained optimization problems.They have parallels t penalty methods.Both methods replace a constrained optimization problem by a series o unconstrained problems and add a penalty term to the objective.The difference betwee both methods consists in their penalty terms.The augmented Lagrangian method add the term with the Lagrange multiplier.The augmented functional for the optimizatio problem ( 10) and ( 11) reads: The derivation of the necessary optimality conditions is based on the concept of fun tionals and the functional derivative.A mathematically precise discussion of the fun

Optimization Problem and Isoper
Based on the above solution, the beam will be optimized.The optimal d of the beam is denoted by the capital shaft of the constant volume 2 th buckling.If the density of the materia  = ,   = () The new variables will be introd  =   , With the new variables, the form Particularly, if both principal mo the factors are:  =  = 1.For exam section in form of an equilateral trian The distribution of material alon is of constant volume and provides t variational methods will be applied f tions.The augmented Lagrangian is re multipliers [10].Augmented Lagrang alytic techniques for solving constrai penalty methods.Both methods repla unconstrained problems and add a p both methods consists in their penal the term with the Lagrange multipli problem ( 10) and ( 11) reads: ℒ The derivation of the necessary o 2 = (4π) −1 ≈ 0.07957.The shaft with the circular, simply connected cross-section delivers correspondingly the minimum for the critical eigenvalue: Two other cases describe the situations in which the form of the transverse cross-section undergoes the transformation, such that one of the geometric dimensions of the cross-section changes.For the technically important case of the thin-walled tubes with the variable thickness of wall t(x) and the mean diameter of tube D(x), the second moments of: tD 3 for t D The case α = 1 corresponds to the adjustable wall thickness (as the design variable) and constant mean diameter of the tube.The analogous dependence arises in the case of the bending of sandwich beams [7].The case α = 3 corresponds the inconstant mean diameter of the tubular shaft (as the design variable) and constant wall thickness.
Consider the shaft hinged on two supports x = L and x = −L: The buckling equations of the shaft read: For the boundary value problem ( 7) and ( 8), the actual curvatures and displacements y * , z * minimize the quotient: Particularly, if both pri the factors are:  =  = 1.

section in form of an equilat
The distribution of mat is of constant volume and p variational methods will be tions.The augmented Lagran multipliers [10].Augmented alytic techniques for solving penalty methods.Both meth unconstrained problems and both methods consists in th the term with the Lagrange problem ( 10) and ( 11) reads: The derivation of the ne tionals and the functional d tional derivative can be foun from the first functional deri

Optimization Problem and Isoperimetric Inequality for Stability
Based on the above solution, the distribution of material along the length of a twisted beam will be optimized.The optimal distribution of the cross-sectional area along the span of the beam is denoted by the capital letter S(x).The optimization problem examines the shaft of the constant volume 2V 0 that supports the maximal moment without spatial buckling.If the density of the material is ρ, the mass 2M of the shaft is given by: The new variables will be introduced: With the new variables, the formulation of the optimization problem is the following:

Optimization Problem and Isoperimetric Inequality for Stability
Based on the above solution, the distribution of material along the length of a beam will be optimized.The optimal distribution of the cross-sectional area along of the beam is denoted by the capital letter ().The optimization problem exam shaft of the constant volume 2 that supports the maximal moment withou buckling.If the density of the material is , the mass 2 of the shaft is given by:  = ,   = ()  () > 0.
The new variables will be introduced:  =   ,  =   ,  = With the new variables, the formulation of the optimization problem is the fo Particularly, if both principal moments of the second order are equal  () the factors are:  =  = 1.For example, this setting is valid for the shaft with th section in form of an equilateral triangle or for the shaft with a circular cross-sect The distribution of material along the length of a shaft is optimized so that t is of constant volume and provides the maximal moment without spatial buckl variational methods will be applied for the derivation of the necessary optimalit tions.The augmented Lagrangian is related to, but not identical to the method of L multipliers [10].Augmented Lagrangian methods are a certain class of algorithm alytic techniques for solving constrained optimization problems.They have par penalty methods.Both methods replace a constrained optimization problem by a unconstrained problems and add a penalty term to the objective.The difference both methods consists in their penalty terms.The augmented Lagrangian meth the term with the Lagrange multiplier.The augmented functional for the optim problem ( 10) and ( 11) reads: Particularly, if both principal moments of the second order are equal I 1 (x) = I 2 (x), the factors are: η 1 = η 2 = 1.For example, this setting is valid for the shaft with the cross-section in form of an equilateral triangle or for the shaft with a circular cross-section.
The distribution of material along the length of a shaft is optimized so that the beam is of constant volume and provides the maximal moment without spatial buckling.The variational methods will be applied for the derivation of the necessary optimality conditions.The augmented Lagrangian is related to, but not identical to the method of Lagrange multipliers [10].Augmented Lagrangian methods are a certain class of algorithms or analytic techniques for solving constrained optimization problems.They have parallels to penalty methods.Both methods replace a constrained optimization problem by a series of unconstrained problems and add a penalty term to the objective.The difference between both methods consists in their penalty terms.The augmented Lagrangian method adds the term with the Lagrange multiplier.The augmented functional for the optimization problem (10) and ( 11) reads: For the boundary value problem (7  * ,  * minimize the quotient: ,  , ,

Optimization Problem and Isoperim
Based on the above solution, the di beam will be optimized.The optimal dis of the beam is denoted by the capital le shaft of the constant volume 2 tha buckling.If the density of the material The new variables will be introduc With the new variables, the formul Particularly, if both principal mom the factors are:  =  = 1.For examp section in form of an equilateral triangl The distribution of material along is of constant volume and provides the variational methods will be applied for tions.The augmented Lagrangian is rela multipliers [10].Augmented Lagrangia alytic techniques for solving constraine penalty methods.Both methods replace unconstrained problems and add a pen both methods consists in their penalty the term with the Lagrange multiplier problem (10) and (11) reads: The derivation of the necessary opt tionals and the functional derivative.tional derivative can be found in [10] o from the first functional derivative of th The derivation of the necessary optimality conditions is based on the concept of functionals and the functional derivative.A mathematically precise discussion of the functional derivative can be found in [10] or [17].The necessary optimality condition follows from the first functional derivative of the augmented Lagrangian (12): In Equation ( 13), c an auxiliary constant, that is proportional to Lagrange multiplier λ.The spatial curvature is signed in Equation ( 13) as K, where the curvatures of the deformed axis of the shaft in the planes 0xy and 0xz are: The application of Fermat's principle for the optimization problems leads to the necessary optimality condition: The necessary optimality condition ( 14) is the requirement that the augmented Lagrangian has a stationary value.The optimal second moment of the cross-section follows from ( 14) as:

Optimization Problem and Isoperimetric Inequality for
Based on the above solution, the distribution of materia beam will be optimized.The optimal distribution of the cross of the beam is denoted by the capital letter ().The optimi shaft of the constant volume 2 that supports the maxim buckling.If the density of the material is , the mass 2 of  = ,   = ()  () > 0.
The new variables will be introduced: With the new variables, the formulation of the optimiza Particularly, if both principal moments of the second o the factors are:  =  = 1.For example, this setting is valid section in form of an equilateral triangle or for the shaft with The distribution of material along the length of a shaft i is of constant volume and provides the maximal moment w variational methods will be applied for the derivation of the tions.The augmented Lagrangian is related to, but not identic multipliers [10].Augmented Lagrangian methods are a cert alytic techniques for solving constrained optimization prob penalty methods.Both methods replace a constrained optim unconstrained problems and add a penalty term to the objec both methods consists in their penalty terms.The augment the term with the Lagrange multiplier.The augmented fun problem (10) and (11) reads: The derivation of the necessary optimality conditions is tionals and the functional derivative.A mathematically pr tional derivative can be found in [10] or [17].The necessary from the first functional derivative of the augmented Lagran The applied method of scaling allows the arbitrary selection of the constant c.For the briefness of the governing equations, the constant c is introduced as the positive solution of the equation: Eα For the boundary value problem ( 7) and ( 8), the  * ,  * minimize the quotient: ,  , ,  ≡    ,

Optimization Problem and Isoperimetric Inequ
Based on the above solution, the distribution of beam will be optimized.The optimal distribution of of the beam is denoted by the capital letter ().Th shaft of the constant volume 2 that supports th buckling.If the density of the material is , the mas  = ,   = ()  () The new variables will be introduced: With the new variables, the formulation of the o Particularly, if both principal moments of the the factors are:  =  = 1.For example, this settin section in form of an equilateral triangle or for the s The distribution of material along the length of is of constant volume and provides the maximal m variational methods will be applied for the derivati tions.The augmented Lagrangian is related to, but n multipliers [10].Augmented Lagrangian methods a alytic techniques for solving constrained optimizat penalty methods.Both methods replace a constraine unconstrained problems and add a penalty term to both methods consists in their penalty terms.The the term with the Lagrange multiplier.The augme problem ( 10) and ( 11) reads: The bending stiffness reduces with this choice to:

Closed-Form Solution of the Governing Equations
The order of both Equations ( 8) could be reduced by one using the boundary value conditions (7).The buckling equations transform with Equation ( 16) to: Two Equations ( 17) with the corresponding boundary conditions for Z and Y are simultaneous nonlinear ordinary differential equations of the second order.The standard methods for the closed form solutions of such equations were apparently not reported in the literature.The principal task of the actual manuscript is to obtain a closed form solution of Equation ( 17).As the first step for the solution of two simultaneous Equations (17), the unknowns will be represented over two new functions K(x), θ(x): The substitution of Equation (18) in Equation ( 17) reduces it to two simultaneous equations for the new unknowns K and θ: The area of the optimal cross-sections S vanishes in the ends of the shaft, where the bending moment and spatial curvature K disappear.This requirement plays the role of the boundary conditions for K and S at x = ±L.The function K(x) = K(−x) must be even, and the function is odd θ(x) = −θ(−x).Thus, from the symmetry considerations are as follows θ(0) = 0.
It is principally possible to solve Equation ( 19) for K, θ, but it is preferable to find the solution for the functions S, θ.With this choice, the solution delivers the optimal cross-section S directly.For this purpose, K in Equation ( 19) is replaced by S using Equation ( 14): For a given S(x), the solution to Equation (21) with respect to θ(x) reads: This choice of the lower integration limit satisfies the symmetry condition: θ(0) = 0. Substitution of θ(x) into the Equation (20) leads to an equation in terms of S only: The Equation ( 23) is an ordinary nonlinear differential equation of the second order.For the solution of Equation ( 23), the dependent and independent variables must be exchanged: The integral (31) evaluates in terms of the beta-function [19]: For volume evaluation, the integrand in Equation (31) must be multiplied by area [18]: For estimation of elastic energy E , the integrand in Equation ( 28) must be multiplied by area t α : The normalized coordinate X = x/L is introduced for the briefness of formulas.The shape S(X) of the right side is given by the implicit function of the normalized coordinate −1 ≤ X ≤ 1 [18,19]: The expression for the volume of the half of the shaft with the unit length is: The following formula is valid for an assessment of elastic energy E in the normalized coordinates: For practically interesting cases, the shape reduces to the elementary functions, as shown in Table 1.
Table 1.Normalized cross-sectional area along the span of the shaft, the normalized volume, and normalized elastic energy.

Exponent Factor Cross-Sectional Area Normalized Volume Normalized Elastic Energy
Ratio Elastic Energy to Volume The buckling shape of the twisted column with the simply supported (hinged) ends is demonstrated on Figure 1 for the exponent factor α = 2.For the other values of parameter α, the buckling shapes are similar.Figure 2 displays the areas of cross-sections of the optimal twisted columns for the exponent factors α = 1, 2, 3. Second moments of inertia of the cross-sections of the optimal twisted columns are shown for the exponent factors α = 1, 2, 3 on Figure 3.

Effectiveness
The influence of the exponent  influences the estimations for the optimization effects.Another argument for the introduction of the invariant optimization factors is methodical.In variational calculus, it is common to obtain one factor as the optimization objective and the others as the a priori given constraints.To convert it into an unconstrained problem, the method of Lagrange multipliers is commonly used.The resulting unconstrained problem with Lagrange multipliers increases the number of variables.The new number of unknown variables is the original number of variables plus the original number of constraints.The constraints are usually solved for some of the variables in terms of the others, and the former can be substituted out of the objective function, leaving an unconstrained problem in a smaller number of variables.This method of solution leads to the nonlinear algebraic equations for Lagrange multipliers.These nonlinear equations in most cases do not possess the closed analytical solutions and are solvable only numerically.

Effectiveness
The influence of the exponent α influences the estimations for the optimization effects.Another argument for the introduction of the invariant optimization factors is methodical.In variational calculus, it is common to obtain one factor as the optimization objective and the others as the a priori given constraints.To convert it into an unconstrained problem, the method of Lagrange multipliers is commonly used.The resulting unconstrained problem with Lagrange multipliers increases the number of variables.The new number of unknown variables is the original number of variables plus the original number of constraints.The constraints are usually solved for some of the variables in terms of the others, and the former can be substituted out of the objective function, leaving an unconstrained problem in a smaller number of variables.This method of solution leads to the nonlinear algebraic equations for Lagrange multipliers.These nonlinear equations in most cases do not possess the closed analytical solutions and are solvable only numerically.
The common method for the solution of the optimization problems with stability constraints uses Lagrange multipliers.In the present manuscript, the method of dimensionless factors will be applied for the optimization analysis.This method delivers different lengths and volumes of the optimal shafts.Instead of seeking for the twisted beams of the fixed length and volume, the method of dimensionless factors allows for the comparison of the shafts with the different lengths and cross-sections.
Consider the shafts with the same form of cross-sections.The affine transformation of the shaft is the product of two elementary transformations, namely homothety and scaling.The homothety of ratio ζ multiplies lengths by ζ.Thus, ζ is the ratio of magnification, or dilation factor, or scale factor, or similitude ratio.The cross-section function S(ξ) scales by another factor , such that for the affine transformed shaft, the cross-section function will be S(ξ).Apparently, the eigenvalue Λ alters in the course of the affine transformation of the shaft.
For each fixed value of α, the two dimensionless factors will be studied: The factors F V , F E (38) will be used for the comparisons of different designs.For arbitrary powers p 1 , p 2 , p 3 , p 4 , the factors will alter for any affine transformation of the beam.The critical buckling moment M = k α EΛ inherits the factor k α and is proportional to this value.Evidently, that the ratios of the buckling loads for different designs with the same form of the cross-sections do not depend on the constants k α .For different cross-sections, the actual value of k α has to be used.
With the methods of dimensional analysis, the characteristic choice of powers is the following: Using the expressions (38) and ( 39), the factors result to: With the choice of the powers (38), the factors F V , F E do not alter for any affine transformation of the shaft.In other words, the factors F V , F E are the invariants to the affine transformation of the shaft and provide a natural basis for the comparison of different designs.
With the above factor, the estimation of the effect of mass optimization turns out to be trivial.The reference design will be selected with the constant cross-section along the span.The invariant factors for the reference design is F V .The factor is the same for all exponents α and for the boundary conditions with both hinged ends.The factor F V reads for the reference design with the constant cross-section as: The greater the factor, the higher the buckling moment for the given length and volume of the shaft.For example, the buckling force of the reference clamped shaft is four times the buckling force of the reference shaft with the hinged ends.
The dual formulations are typical of the optimization of buckling as well.For the dual formulations, the masses of the shafts for the fixed lengths and fixed buckling moments are compared.The volumes and masses of the optimal and reference shafts relate to each other as the inverse roots of the order α of the factors F V : Specifically, the shaft with the higher value of the factor F V possesses the lower mass.Consequently, the optimality is expressed for the shaft in the form of the isoperimetric inequalities.
The results of the evaluation of the dimensionless volume factors and volumes for the fixed critical eigenvalue Λ = π/2 and the half-length L = 1 are presented in Table 2.The half-volumes of the optimal columns V and the volume optimization factor F V are shown as the functions α of on Figure 4. Figure 5 demonstrates the half-volumes V and elastic energy E of the optimal columns.Table 2 displays the dimensionless volume factor and volume for the fixed critical eigenvalue Λ = π/2 and the half-length L = 1 for different exponent factors α.The reference serves the shaft of the constant cross section with the same critical eigenvalue Λ = π/2.The area of its cross-section is one, such that the half volume of the reference shaft is 1.
Specifically, the shaft with the higher value of the factor  possesses the low mass.Consequently, the optimality is expressed for the shaft in the form of th isoperimetric inequalities.
The results of the evaluation of the dimensionless volume factors and volumes fo the fixed critical eigenvalue  = /2 and the half-length  = 1 are presented in Table The half-volumes of the optimal columns  and the volume optimization factor  a shown as the functions  of on Figure 4. Figure 5 demonstrates the half-volumes  an elastic energy ℰ of the optimal columns.Table 2 displays the dimensionless volum factor and volume for the fixed critical eigenvalue  = /2 and the half-length  = 1 fo different exponent factors .The reference serves the shaft of the constant cross sectio with the same critical eigenvalue  = /2.The area of its cross-section is one, such th the half volume of the reference shaft is 1.

Conclusions
For the shaft, the optimal shape along the axis was searched.The cross section tha delivers the maximum or the minimum for the critical eigenvalue was to be determined among all convex and simply connected domains.At the beginning, the validity of th static Euler's approach for the simply supported beam (hinged), twisted by th conservative moment, was demonstrated.The applied method for integration of th optimization criteria delivers different length and volumes of the optimal twisted beams Instead of seeking for the twisted shafts of the fixed length and volume, the twisted beam with the different lengths and cross-sections using the invariant factors could be directl compared.The moments of inertia of the rod are the powers of the cross-sectional are with the certain exponent.The exponents are the positive rational numbers, which ar equal to or greater than one.For all admissible exponents, the optimal distribution of th cross-sectional area for the shaft was given in the closed form.In the general case, th optimal distribution is stated in terms of the higher transcendental functions.Particularly for exponents 1,2, and 3, the solution reduces to certain elementary functions.The fina formulas involve the length of the beam, its volume, and critical torque.Remarkably, i the torsion stability problem, the optimal shape of the beam is roughly parabolic along it length.In the torsion stability problem, the optimal shape of the simply connected, conve cross-section is the equilateral triangle.Using the methods of isoperimetric inequalitie [20], the exact boundaries for the buckling eigenvalues for the fixed volume and length o the hinged twisted shaft could be stated.

Conclusions
For the shaft, the optimal shape along the axis was searched.The cross section that delivers the maximum or the minimum for the critical eigenvalue was to be determined among all convex and simply connected domains.At the beginning, the validity of the static Euler's approach for the simply supported beam (hinged), twisted by the conservative moment, was demonstrated.The applied method for integration of the optimization criteria delivers different length and volumes of the optimal twisted beams.Instead of seeking for the twisted shafts of the fixed length and volume, the twisted beams with the different lengths and cross-sections using the invariant factors could be directly compared.The moments of inertia of the rod are the powers of the cross-sectional area with the certain exponent.The exponents are the positive rational numbers, which are equal to or greater than one.For all admissible exponents, the optimal distribution of the cross-sectional area for the shaft was given in the closed form.In the general case, the optimal distribution is stated in terms of the higher transcendental functions.Particularly, for exponents 1, 2, and 3, the solution reduces to certain elementary functions.The final formulas involve the length of the beam, its volume, and critical torque.Remarkably, in the torsion stability problem, the optimal shape of the beam is roughly parabolic along its length.In the torsion stability problem, the optimal shape of the simply connected, convex cross-section is the equilateral triangle.Using the methods of isoperimetric inequalities [20], the exact boundaries for the buckling eigenvalues for the fixed volume and length of the hinged twisted shaft could be stated.
From the presented closed-form solution arises the open questions for the future work.It is well known that different types of boundary conditions on both ends could essentially influence the character of the optimal solution.If the Sturm-Liouville conditions are not satisfied, multimodal buckling regimes are principally possible.The study of the optimization problems for the multimodal buckling could be investigated for the pure twisted and twisted-compressed beams with the presented analytical solution methods.Another task could be the study of the multiplication of supported beams of finite length with the intermediate supports along the span of the beam.Special attention is required for the optimization problem for an infinite, periodically supported shaft.
The applied Euler approach is applicable for the conservative systems and for special classes of non-conservative systems.If the loads are non-conservative, the optimization problems require the application of dynamic stability methods.For example, if the moment direction follows the direction of the deformed axis, the loads be non-conservative.Seemingly, the closed form solutions of the nonconservative optimization problems are unknown.An interesting question is whether any type of the nonconservative torsion load acknowledges the closed-form solution of the optimization problem.

Figure 1 .
Figure 1.The shaft with the simply supported (hinged) ends and its buckling curve.The projections (), () of the spatial buckling curve in Cartesian coordinates are shown in the graphs below.

Figure 1 .
Figure 1.The shaft with the simply supported (hinged) ends and its buckling curve.The projections y(x), z(x) of the spatial buckling curve in Cartesian coordinates are shown in the graphs below.

2 = √ 3 /
18 ≈ 0.9622.For the circular cross-section, the constant is Appl.Mech.2023, 4, FOR PEER REVIEW Consider the shaft hinged on two () = 0,   () = 0 () (−) = 0,   (−) = 0,  The buckling equations of the sh (   ) = −  , For the boundary value problem  * ,  * minimize the quotient: Based on the above solu beam will be optimized.The of the beam is denoted by th shaft of the constant volum buckling.If the density of th = ,   =The new variables will =   ,With the new variables,

Figure 3 .
Figure 3. Second moments of inertia of the cross-sections of the optimal shafts for  = 1,2,3.

Figure 3 .
Figure 3. Second moments of inertia of the cross-sections of the optimal shafts for α = 1, 2, 3.

Table 2 .
Dimensionless volume factor and volume for the fixed critical eigenvalue Λ = π/2 and the half-length L = 1.

Figure 4 .Figure 4 .
Figure 4. Half-volumes of the optimal shafts  and the volume optimization factor  .

Figure 5 .
Figure 5. Half-volumes  and elastic energy ℰ of the optimal shafts.

5 .
Half-volumes V and elastic energy E of the optimal shafts.