Research Article | Open Access
Mariano Vázquez Espí, Jaime Cervera Bravo, Carlos Olmedo Rojas, "On Galileo’s Tallest Column", Mathematical Problems in Engineering, vol. 2015, Article ID 649341, 10 pages, 2015. https://doi.org/10.1155/2015/649341
On Galileo’s Tallest Column
The height at which an unloaded column will fail under its own weight was calculated for first time by Galileo for cylindrical columns. Galileo questioned himself if there exists a shape function for the cross section of the column with which it can attain a greater height than the cylindrical column. The problem is not solved since then, although the definition of the so named “constant maximum strength” solids seems to give an affirmative answer to Galileo’s question, in the form of shapes which seem to attain infinite height, even when loaded with a useful load at the top. The main contribution of this work is to show that Galileo’s problem is (i) an important problem for structural design theory of buildings and other structures, (ii) not solved by the time being in any sense, and (iii) an interesting problem for mathematicians involved in related but very different problems (as Euler’s tallest column). A contemporary formulation of the problem is included as a result of a research on the subject.
In 1638, in his Discorsi e Dimostrazioni Matematiche , Galileo postulated the existence of the tallest column, that is, a cylindrical column, such that it attains the maximum height once the area of its cross section and the strength of its material are prescribed. Therefore, Galileo’s tallest column is in the limit of resistance only bearing its own weight. The rationale of his proof gave rise to the square-cube law, a mathematical principle that considers the relationship between the flow through the surface of a volume and the stock into the latter, in the mechanical case, for example, the stress with the weight. This principle has been very useful and applied in a variety of scientific fields, mainly in biology [2–4].
In the following century, Euler  pointed out a very different problem, that is, to find the shape of a stable column, axially symmetric with respect to the vertical axis, such that it attains the maximum height once its volume, specific weight, and Young’s module are given, buckling due to the action of a load bear at its top.
In our view, both problems are not yet completely solved nowadays, although Galileo’s problem has received much less attention than Euler’s. Furthermore, we think that Galileo’s problem is more meaningful for a theory of design of structures subjected to small limit on strains and displacements, as buildings and other structures in civil engineering [6, 7].
Hereafter, we consider the solid continuum with the following standard assumptions:(1)The scope of the analysis is the classical theory of Elasticity.(2)The process of deformation is isothermal and quasi-static; heat or kinetic energy is not taken in consideration into energy balancing.(3)We are only interested in solutions whose strains and displacements will be very small; hence equilibrium and compatibility equations approximately hold in the geometry of the undeformed body.
Section 2 outlines important aspects of Galileo’s problem, comparing it to Euler’s and enlightening its importance and profound meaning for a theory of structural design. Our main working hypothesis is formulated there. Section 3 deals with some clues that support our working hypothesis; that is, an insurmountable size exists for a fairly large set of structural problems, as it is the case of cylindrical columns of Galileo’s first insight into the tallest column. This section covers the main aim of this paper which is to attract mathematicians to work out Galileo’s problem, because we are architects and our mathematical knowledge is, to say the least, limited, but to solve the problem is a key point to continue the development of structural design theory. Finally, Section 4 is devoted to formulate the problem formally in contemporary terms.
2. Galileo’s Problem on the Tallest Column
Proposition VII. Among heavy prisms and cylinders of similar figure, there is one and only one which under the stress of its own weight lies just on the limit between breaking and not breaking so that every larger one is unable to carry the load of its own weight and breaks while every smaller one is able to withstand some additional force tending to break it. (GALILEO, 1638)
Consider a cylindrical column of height and diameter , of a lineal elastic material defined by Young’s Modulus , allowable compressive stress , and specific weight and subjected to the sole action of its own weight.
Such a column will be unsafe in simple compression if the applied load exceeds the column strength , where is the area of the cross section—if the base is a circle, —or, which is the same, if the applied stress exceeds . This fact means that the height of such a column may not be greater than a characteristic length of the material, . We name this length “structural scope” of the material, . And we name “structural scope” of cylindrical columns to the maximum height of safe columns. In this simple case, is numerically equal to the material scope , but generally is related to but is not equal to . Therefore the first conclusion of Galileo can be expressed asLater, Galileo considers in which way this insurmountable limit can be increased. He envisaged two main ways: (a) to increase the material scope , or (b) to change the shape of the column. In the latter case, he reasons—in a funny paragraph—that if the giants exist, they would have a very different aspect and proportions compared to human beings; specifically the bones of their legs would have a greater diameter/length ratio, because otherwise their weight that increases in proportion to would be greater in proportion to their strength that increases as , and as a result the giants would suffer stresses—which increase as —greater than human beings, and considering the bone material very similar in all the mammals, the giants would be unable to perform in their life as well as human beings. A few centuries later, this result could be confirmed comparing dinosaurs of different sizes but of same suborder or family (i.e., Theropoda or Tyrannosauridae) .
As a material with infinite strength or null specific weight does not exist, it is clear that following first Galileo’s way we can only increase the insurmountable height but remaining finite. If we adopt the latter way, the main question arises: does there exist an optimal shape which has infinite height? We are looking for an answer to this question because it is a key into the theory of structural design. If the answer is “Yes,” then Galileo’s problem has a solution for any size considered. But if the answer is “No,” there exist instances of structural problems which have no solution; that is, there are unsolvable problems in structural design. Furthermore, as we will show below, near the unsurmountable size, any solution for the problem will have an unaffordable physical cost, so it would be infeasible from a practical view.
Our working hypothesis is that a finite insurmountable size exists for a fairly large set of structural problems (not only for Galileo’s problem). Moreover, the optimal shape for each problem—which maximises the finite insurmountable size—is a sound reference to measure the efficiency of all other shapes with size lesser than the one of the optimal shape [7, 9, 10].
In a first approximation, we can represent the physical cost of a structure by its self-weight, as many costs during the manufacturing, but not all, are approximately proportional to the self-weight of the structure: CO2 emissions, mineral resources consumption, and so forth. For a given structural problem, we define the structural efficiency as the ratio between the useful load and the total load (i.e., the useful load plus the self-weight) required to solve that problem in a particular structure. Galileo postulated also the relationship between the size of a structure and its ability to resist a useful load. Let us consider a cylindrical column of size . It can resist an additional useful load , the value of which is at most the weight difference between this column and the column of insurmountable height . Hence, according to the previous definition, the efficiency of such a column will beNote that (2), which we name Galileo’s rule, is exact in the case of cylindrical columns, but it is not proved that it would be a general rule. The best result we get up to date is that Galileo’s rule is a very good estimate in canonical problems like bending of beams and bridges . We define the load cost as the inverse of efficiency, hence always higher than unity, . Then, the self-weight of the column isAs a reward, Galileo’s rule, apart from the cost, gives us a sound estimate for the self-weight that it is a required datum for the final project but unknown in the preliminary phases of design of large structures.
2.1. Comparison between Euler’s and Galileo’s Problems
Remember the cylindrical column of height and diameter . As we saw, such a column will be unsafe if its height is equal to or greater than material scope . But the column may also fail by elastic buckling. According to Landau and Lifshitz’ Course , the critical height for buckling is related to the diameter byThe ratio is another characteristic length of the structural material. Whereas the scope is its specific strength, is its specific stiffness. Let us define the geometrical slenderness of the column as the ratio . ThenTherefore, the safety of a given column bearing only its own weight requires the fact that two conditions hold: (i) and (ii) . It is worth noting that it is always possible to satisfy the second condition, as for each height we can choose such that . However the first condition is an absolute one, as it only depends on the properties of material. Hence, the height of a safe, cylindrical column would be lesser than or equal to .
This limit, as noted above, only could be modified in two ways: changing material’s properties or changing the shape of the column (or both). The interesting point here is that to answer Galileo’s question we must elucidate if a finite structural scope (related to the material scope ) exists for any shape of the column. The advances in the analysis of the solution for Euler’s problem [12–14] are useless to this aim, mainly because of ignoring the limit that strength condition imposes on the shape (condition (i) for cylindrical columns).
Furthermore, as it is well-known, the classical solutions of Euler’s problems on buckling are contradicting the experimental data. “This contradiction between theory and experiments is not surprising. The ideal appearance of a phenomenon is always more or less influenced in practice by multiple causes that can deform it to the point of leaving none but a caricature. In the problems of instability, the theory considers only perfect elements, both form and structure indefinitely elastic and resistant. The test pieces, or the members actually built, are very far from perfection; the materials are inhomogeneous, and they are approximately elastic and within certain limits.”  As a consequence, in the engineering practice, there is no bifurcation between two different equilibria. On the contrary, as the slenderness approaches to its critical value , the failure changes continuously from simple flattening to bending with net compression. In fact, the so named critical load (or Euler’s load) is not a “load” at all but a stiffness of the column against lateral displacement, and the failure occurs for real loads numerically lesser than this stiffness. (Unfortunately, this stiffness can be expressed in load units (), but it is better understood with stiffness unit (), showing that it is a ratio between the bending moment and the lateral displacement.) Moreover, although the buckling of real structural members (with negligible self-weight) is a nonlineal problem in a first, mathematical view, therefore candidate to a numerical solution, it is possible to overcome the difficulties and to solve the problem by a direct albeit nonexact formulation when structural strength is taken into account .
3. On the Existence of a Finite Height for Galileo’s Tallest Column
Our epistemological situation confronting Galileo’s problem is analogous to the situation that algorithm designers are confronting when using the well-known Theory of NP-Completeness . We, the structural designers, do not know if a finite insurmountable size exists for the problem at hand; hence we cannot know in advance if our problem is solvable or not. But if we believe that this limit exists, we can manage at least a rough estimate of its value and, armed with this knowledge, take a decision about the solvability of the problem. Indeed, if we know the size limit that different types of structures can reach for our problem, we can evaluate approximately the relative merit of each type and select the most promising one for the actual size of our problem. So the existence (or not) of a finite height for Galileo’s tallest column is a key point for our everyday work.
Let us consider the two main approaches to the problem: first that such a limit does not exist because it is easy to find the corresponding shape and second that such a limit probably exists because it is very hard to find out any shape that can overcome a given finite limit on its height.
3.1. The Known Solutions with Infinite Height Are Unfeasible
Although Euler’s tallest column and related problems have received very much attention up to date, some researchers while studying that problem spent a few minutes to study, too, problems related with that of Galileo. This is the case of Karihaloo and Hemp  that study the “maximum strength design” of structural members. In their approach, all cross sections of a structural member attain the maximum allowable stress for the given material; therefore the solution is also referred to as “constant maximum strength design.”
Let us examine with some detail the constant maximum strength design of Karihaloo and Hemp for tension members; see  (Section 2.1) and Figure 1. Consider a cable of length and cross-sectional area , hung on its top edge () and when its bottom edge () is a free boundary. If the gravity and the external load are acting in the negative direction of -axis, the condition of the constant maximum strength isThe solution isas it can be checked obtaining its derivative and comparing the result with (6). If the cross section is circular, the radius gives us the contour of the member as follows:This solution can have an infinite height with constant stress and bear a useful load at the bottom edge (, ). Anyhow, its volume grows exponentially with its size (), so in practice it is an “intractable” solution—in the same meaning that the term is used in algorithm complexity theory —with a load cost as follows:
Furthermore, the solution is not feasible from the point of view of equilibrium because only the equilibrium in the direction is considered for obtaining a constant stress . Let us consider the 2D-case for the sake of simplicity; see Figure 1(b). As the border (or ) is stress-free, the tangent in any point is a principal direction (), as the orthogonal direction is (). As , we haveand hence if is the angle between the principal stress direction and coordinate axes. This value is greater than allowable stress for all . In fact, as increases, exponentially, and grows in the same way. As a consequence, the classical solution is not a feasible one for common failure criteria, that is, as Von Mises criterion. This error is common to all solutions obtained making use of Bernoulli-Euler theory or Navier hypothesis, as these models are not applicable in the case of variable cross section beams—a well-known drawback; see, for example, . These solutions do not prove in any way that an infinite size for columns or beams would be feasible. Of course, these solutions are almost exact when the size is much lesser than the material structural scope , as the exponential function grows very slowly when its argument is very small, and for small size the solutions obtained pass fairly well through experimental checks . But they are useless to answer the question pointed out by Galileo, because to this issue we must explore sizes of the same order of magnitude that material scope, that is, large structures for which the equilibrium equations must be completely fulfilled.
3.2. Trying to Refute the Existence of a Finite Height
In 2010 (unpublished work), we tried to refute the existence of a finite structural scope whatever the structural material of Galileo’s column is. Our try was naive and unsuccessful. However we think it can help others to understand the difficulties of the problem and perhaps gives clues to better searches for a complete solution. To be short, let us consider a 2D-version.
We have an elastic linear material as before. We choose Von Mises criterion as the failure one; hence in any point of a feasible body the following expression hold:Let us suppose that we are able to determine a suitable stress field for an instance of Galileo’s problem, defined over all the plane. This field implicitly defined the shape of a body for which (11) holds. The support line must intersect in some point of some solution ofLet and be the principal planes of the stress tensor. Then the curves that can form the stress-free contour of the column will be solutions of one of the following equations:depending on if we are looking for a compressive or tensile solution. Furthermore, these curves intersect in some point of solution (12) and also the solution ofIndeed, if in (13) we select the direction , we must select now in (14) and vice versa. This intersection point will be named “base-end vertex” hereafter. Selecting an appropriate set of arcs from solutions of (13), with the additional condition that they define a closed region on the plane together with the support line, we get the shape of a solution for this instance of Galileo’s problem and we can determine its height relative to the material scope; see Figure 2.
If we would be able to explore completely the set of all possible stress fields—and, of course, we are not—with the maximal height obtained we could answer—“Yes” or “No”—Galileo’s question.
3.2.1. Generating Subsets of Stress Fields
To check the possibility of our idea we recall on well-known Airy’s function. Whatever Airy’s stress function is, for which the biharmonic equation holds, it can be defined aswhere and are analytical functions in . The function satisfies both compatibility and equilibrium equations, and the displacements can be calculated without integration of the stress functions, from the complex potentials and .
The stress field isUsing the Hooke law with this field, the displacement functions arewhere is Poisson’s modulus.
The complex potentials considered were simple polynomials as follows:The coordinate origin will be at the top of a column (or at the bottom of a cable). Some boundary conditions must hold always as follows:We considered also different supports as follows:With a concrete selection of boundary conditions and support geometry, an instance of the problem is defined. We considered several Standard Problem Definitions as follows.
Mountain or Peak. A shortening of the shape can be measured with . Consider
MountainBis. A shortening of the shape can be measured with . Consider
MountainIII or Summit. A shortening of the shape can be measured with . Consider
The trial stress field is derived from and , which are defined through some constant parameters to be determined. In this way, is completely defined. Let us name to the set of parameters to be determined for the expressions of , , , and . These parameters, once determined, will define the body (or bodies) generated by the trial stress field.
3.2.2. Getting a Feasible Body
Once a trial stress field is selected in algebraic form, (19) and (20) must hold simultaneously for all . Let us represent this condition with the set of problem equations. One hasIt is worth noting that, into (24), all the equations needed for satisfying (19) or (20) for all must be included. For example, if one of (20) isthen the following three equations will be included into (24). Consider
Solving (24) for all , we obtain a set of relationships between the parameters of that we can write aswhere is the set of independent parameters for a given such that . It is worth noting that some (or all) parameters in can disappear from the stress or displacement fields, because only derivatives of are present in the expressions of these fields.
If the stress field depends on some components of , these components can be chosen freely as equilibrium or compatibility equations will hold for arbitrary values of . This fact means that the given represent a family or set of solutions, not a unique one.
Remeber that is the comparison stress, and it must be calculated to assert if the solution is feasible or not; that is, . The given or unknown line of the support (straight or curve) must intersect some solutions ofThis requirement is necessary if the body attains its maximal resistance.
Let and be the principal planes of the stress tensor. Then the curves that can form the contour free of the stress will be solutions of one of the following equations:We must choose one of the two, depending on if we look for a column (compression) or a cable (tension). The solutions would be curves of the form .
Selecting support curves with free contour curves in such a manner that they form a closed domain, a set of shapes is determined. For each shape, the safety criterion must be imposed in any interior point. In this phase, some parameters of could be dependent on others, and in this case a new reduction of the number of independent parameters results in the following:If all the parameters become determined now, the stress field corresponds to a unique shape. Otherwise, we have a family of shapes depending on arbitrary values of .
Perhaps, the main problem is to determine the base-end vertex mentioned before. If the curves that define the shape can be obtained in an explicit form (), they can be managed directly. Otherwise, the shape will be defined by inequations, whose sign can be determined for each function considering the sign of its value in , where is the size used in (20), if the origin was specified as a point of the free contour (null surface stress) in (19).
Anyhow, it is always possible simply to draw the curves of the contour of the body and the support line for an arbitrarily chosen value of the size . Generally, the last one will intersect the former and in two different points. But these two points may be the same base-end vertex of the contour; the first point mentioned is the vertex of the base where the stress-free contour ends; the second one is the vertex where the own base ends. Hence, adjusting the value of in such a manner that the two points be the same, we will find out the correct base-end vertex and the height of the shape by a simple albeit tedious trial-and-error method.
3.2.3. An Example in Detail
Bodies. There is a family dependent on parameter of . This parameter defines the abscissa of the base-end vertex.(i)Vertices: ; ; For a normal steel (see Table 1), the real domain of is approximately with roots in the extremes of the interval and in . Consider(ii)Stress-free arcs: they do not have analytical expressions. Drawing them, it is clear that there are solutions only for . The base width varies from 0,9 up to 1,4. The height of all shapes is constant, only depending on material properties; see Figure 3.(iii)Insurmountable size : .
3.2.4. Provisional Conclusion
After our research, we can claim nothing about Galileo’s question, in spite of the fact that all shapes we found have insurmountable size lesser than Galileo’s column; see Table 1. But at least we have tried to do the best, trying to refute our own conjectures. Maybe it is possible to search on all the set of complex potentials with methods of high mathematics, which of course are beyond our knowledge.
4. Formal Definitions of Galileo’s Problem
After all, we have two main hypotheses that we outline informally as follows.
Conjecture 1. A finite insurmountable size exists for a fairly large set of structural problems (not only for Galileo’s problem) when the self-weight and stress limit are taken into consideration.
The second one is suggested for the results of our search in Section 3.2, and it is stronger than the first one.
Conjecture 2. In case of original Galileo’s problem, the insurmountable size for a solid column (without holes of any kind) is equal to the material scope; that is, .
Refuting any of both hypotheses consists in showing a given problem—including support, boundary conditions, and failure criterion on stresses—and a shape family for structures that can solve it and that includes a shape of infinite size.
4.1. The 2D Galileo Problem
For the sake of simplicity let us consider a 2D-universe.
Problem 1 (see ). To find a -symmetrical body of maximal height, placed in the semiplane and supported in the line, only bearing its own weight, of a homogeneous, linear elastic material defined by Young’s Modulus , Poisson’s modulus , allowable compressive stress , and specific weight and subjected to displacement constraints and and to some criterion on stresses to be fulfilled over all the body that can be expressed asThe support line can be a surface with friction following Coulomb theory. Hence the tangential stress is subjected toThere are not additional fundamental constraints on shape, but someone can be imposed for convenience.
The original Galileo column is simply a rectangular domain of height and whatever width . The principal stresses are and . Further, the Von Mises stress is .
If anyone can envisage a general proof of our conjectures, the related problems would be directly solved. In any case else, the general formulation of the problem can be stated as to find out a shape with infinite size or height, proving in this way that our working hypothesis is false but being the shape of maximum scope determined. We think on this class of problems as good candidates for some topology or shape optimisation methods [22–26].
4.2. Other Related Problems
As the general formulation can be hard to attack with available methods, we can suggest some alternative problems which in our view could be equivalent (or at least approximately equivalent) to the original.
Let be a given volume in the 2D-universe. We can consider the problem of finding a shape with this given volume that maximises the height of the figure subject to the same stress tensor constraint.
Perhaps the stress constraint can be replaced by minimising the (maximum or mean) Von Mises stress in the volume, being the latter unbounded, and the total height of the figure fixed to a given value . With this problem it should be the case that we will get solutions with maximum absolute Von Mises stress lesser than ; hence with appropriated scaling we will get a solution higher than Galileo’s column.
Another approach arises from considering the calculus of the maximum scope of a shape as a limit case. Let us consider a useful load at a height as a uniform load along a width . The problem is now to find a shape of minimal weight in equilibrium with and its self-weight with the stress constraints as above. One additional constraint on the shape will be that it must lie into the region limited by and . If this problem can be solved, the structure scope will be the limit of when or . Obviously, a solution is Galileo’s column of constant width equal to , but is there another one? The useful load can be defined too as , where is a given constant. In this case the function can be viewed as a design variable, or its integral over can be viewed as an additional constraint on stresses.
More equivalents formulations can exist or can be proposed following these lines.
We think that a minimum compliance approach is not equivalent to problems in Galileo’s realm due to self-weight. But it could be the case that minimum compliance objective leads to useful solutions that after appropriate scaling provide that the stress constraint be fulfilled.
Galileo’s problem has theoretical interest in mathematics and very practical interest in the structural design theory. It would be a benchmark problem for topology or shape optimisation methods. Each different stress constraint or material model (e.g., plasticity) leads to new instances of the problem. In 2038 it will be 400 hundred years since Galileo formulated the problem.
We will appreciate any insight into it.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Part of the research work for this paper was carried out during the one-year stay of Mariano Vázquez Espí at CIMNE (Universitat Politécnica de Cataluña, Barcelona), financed by the Universidad Politécnica de Madrid. Thanks are due to Professor Eugenio Oñate and the CIMNE staff for their hospitality. Thanks are due too to the Academic Editor, Reza Jazar, for improving the original paper.
- G. Galilei, Discorsi e Dimostrazioni Matematiche, Elsevier, Leiden, The Netherlands, 1638.
- J. B. S. Haldane, “On being the right size,” Harper's Magazine, vol. 152, pp. 424–427, 1926.
- T. A. McMahon, “Size and shape in biology,” Science, vol. 179, no. 4079, pp. 1201–1204, 1973.
- D. W. Thompson, On Growth and Form, Cambridge University Press, Cambridge, UK, 1942.
- L. Euler, “On the strength of columns,” Memoires de L'Academic des Sciences et Belles-Lettres, vol. 13, pp. 252–282, 1759.
- J. Antuña and M. Vázquez, “Are there unsolvable structural problems? An open question,” Informes de la Construccion, vol. 64, no. 525, pp. 103–109, 2012.
- J. C. Bravo and M. V. Espí, “Galileo, Maxwell, Michell, Aroca: measuring the structural efficiency,” in Structural Milestone in Architecture and Engineering. International Conference on Research in Construction, Madrid, november 2011. IETcc-UPM, Instituto de Ciencias de la Construction (CSIC), 2011, http://oa.upm.es/9931/.
- T. A. McMahon, J. T. Bonner, and W. H. Freeman, On Size and Life, Scientific American Library, New York, NY, USA, 1983.
- J. C. Bravo, “The structures and the self-weight,” Informes de la Construction, vol. 42, no. 407, pp. 73–85, 1990.
- C. O. Rojas, M. V. Espí, and J. C. Bravo, “On the insurmountable size of truss-like structures,” in Proceedings of 3rd International Conference on Mechanical Models in Structural Engineering, pp. 477–493, University of Seville/University of Granada/CSIC, Sevilla, Spain, 2015.
- L. D. Landau and E. M. Lifshitz, Theory of Elasticity, vol. 7 of Course of Theoretical Physics, Pergamon Press, New York, NY, USA, 2nd edition, 1970.
- S. J. Cox and C. McCarthy, “The shape of the tallest column,” SIAM Journal on Mathematical Analysis, vol. 29, no. 3, pp. 517–551, 1998.
- J. I. Díaz and M. Sauvageot, “Euler's tallest column revisited,” Nonlinear Analysis: Real World Applications, vol. 11, no. 4, pp. 2731–2747, 2010.
- Y. V. Egorov, “On the tallest column,” Comptes Rendus Mecanique, vol. 338, no. 5, pp. 266–270, 2010.
- J. Dutheil, “The theory of instability through disturbance of equilibrium,” in Proceedings of the 4th Congress of IABSE, pp. 275–295, Association Internationale des Pont et Charpentes, 1952.
- J. Cervera Bravo, J. Ortiz Herrera, M. Vázquez Espí, and A. Aznar López, “Steel compression design: the buckling's weigth,” Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería, vol. 29, no. 1, pp. 79–91, 2013.
- M. R. Garey and D. S. Johnson, Computers and Intractability: A Guide to NP-Completeness, WH Freeman, New York, NY, USA, 1979.
- B. L. Karihaloo and W. S. Hemp, “Maximum strength/stiffness design of structural members in presence of self-weight,” Proceedings of The Royal Society of London, Series A: Mathematical and Physical Sciences, vol. 389, no. 1796, pp. 119–132, 1983.
- X. Tong, B. Tabarrok, and K. Y. Yeh, “Vibration analysis of Timoshenko beams with non-homogeneity and varying cross-section,” Journal of Sound and Vibration, vol. 186, no. 5, pp. 821–835, 1995.
- R. H. Plaut and L. N. Virgin, “Optimal design of cantilevered elastica for minimum tip deflection under self-weight,” Structural and Multidisciplinary Optimization, vol. 43, no. 5, pp. 657–664, 2011.
- É. Goursat, “On some equations of Monge explicitly integrable,” Annali della Scuola Normale Superiore di Pisa-Classe di Scienze, vol. 1, no. 1-2, pp. 35–59, 1932.
- H. Azegami, S. Fukumoto, and T. Aoyama, “Shape optimization of continua using NURBS as basis functions,” Structural and Multidisciplinary Optimization, vol. 47, no. 2, pp. 247–258, 2013.
- F. Navarrina, I. Muiños, I. Colominas, and M. Casteleiro, “Topology optimization of structures: a minimum weight approach with stress constraints,” Advances in Engineering Software, vol. 36, no. 9, pp. 599–606, 2005.
- J. París, F. Navarrina, I. Colominas, and M. Casteleiro, “Stress constraints sensitivity analysis in structural topology optimization,” Computer Methods in Applied Mechanics and Engineering, vol. 199, no. 33–36, pp. 2110–2122, 2010.
- M. Victoria, P. Martí, and O. M. Querin, “Topology design of two-dimensional continuum structures using isolines,” Computers and Structures, vol. 87, no. 1-2, pp. 101–109, 2009.
- B. Zheng, C. J. Chang, and H. C. Gea, “Topology optimization considering body forces,” International Journal for Simulation and Multidisciplinary Design Optimization, vol. 3, pp. 316–320, 2009.
Copyright © 2015 Mariano Vázquez Espí et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.