This has been published : C.R. Acad. Sci. Paris, t. 311,série II, pp 37-43, 1990. The whole text with the equations may be may be downloaded here.

Abstract

The shape of tree branches is computed using the theory of flexure of beams. Numerical computation is used and takes into account large displacements and growth. The coupling between elastic bending and growth results in a branch shape having an inflection point connected with a permanent deformation.

Growth stresses

The study of stresses in living materials has been the subject of many papers. For example the stresses due to anisotropic growth [1], influence of stresses on growth [2], and even viscoelasticity [3] or plasticity. The influence of gravity on growth is known under the names of gravimorphism [4] or geotropism [5]. The proportions of trees have been found to be limited by elastic criteria to stand under their own weight [6]. Growth stresses and particularly those due to reaction wood play an important role in the mechanics of trees [8], but this paper is focused essentially on the physical influence of gravity. Growth stresses are predominant in the trunk of a vertical tree, where gravitational forces creating nearly hydrostatic pressures may be neglected [8]. In branches, where bending is dominant, stresses due to gravitational forces  are much larger than in the trunk. In this paper, growth stresses will be considered as a distinct phenomenon that will not be taken into account in a first approximation.

Upright growth of trees

A forest tree that has been shifted from its normal upright position during a storm will probably grow straight again. Reaction wood may force the tree to an upright position again, but bending strain was not found to affect radial xylem growth in Douglas-fir [5]. Old branches have an inflection point and their curvature remaining almost unchanged when cut, their deformation is permanent. It is often not possible to straighten these branches without breaking them.

Bending of beams and of  branches

In the absence of gravity a branch would probably grow straight in the direction of light. With growth rings of constant thickness and a constant yearly increase in length, the shape of a branch would be a cone. Trees may be considered as structures made of beams (dead branches) subjected to gravity. Classical Strength of Materials considers a beam with a given shape, applies a load and calculates the resulting shape, slightly different from the original one for thick beams. Thin beams may have large deflections resulting in non linear displacement, though still elastic. To apply elasticity to a tree branch, it is necessary to know its shape before any calculation. The theory shows that the curvature of the bent beam has a constant sign,  there is no inflection point. Elasticity is adequate to calculate the deflection of a branch when the changes in thickness and length are negligible as for a branch temporarily loaded with snow or fruits.

Irreversibility of the bend

When the load is permanent, viscous or plastic deformation may occur, but the coupling between the simultaneous increase of the load and the thickening of the branch is more important. A branch bends continuously with time, even if it thickens. Even with a constant load there would be no decrease of the deflection when the branch thickens. On the contrary, a decrease of the thickness of the branch, would also increase the deflection. The process of growth being time-dependent, the shape of a branch is a function of time.

Method to calculate the branch bending

Dividing time into infinitely small intervals, it is possible to divide each time step in a few more steps: growth, loading and bending. Growth produces an increase in length and thickness of the branch and therefore a small change in geometry. Using the new geometry, the increase in load being small, linear elasticity may be used to compute the deflection, giving another change in geometry of the branch. The growth cycle may then be repeated for the next time step and so on. Because of the changes in thickness and length, the principle of superposition does not apply directly. Therefore, the stress distribution through the branch is no more linear and a permanent deformation occurs even if the incremental stress distribution is linear, e.g. elastic behaviour of the wood.

Numerical method

Many methods have been devised by engineers to analyse the mechanical behaviour of structures. Numerous configurations have been solved with analytical formulae, but they are valid only for simple geometries. With the advent of computers, numerical methods, such as finite elements and finite differences have been developed to calculate complicated structures. None of them (at our knowledge), takes into account the growth phenomenon (occurring in the same manner in the construction of buildings as for trees). In order to calculate the shape of a branch, a simple finite difference method has been used to integrate the differential equation of flexure step by step along the branch and iterated  in a time marching process. At each time step, thickness, length and deflection of the branch are adjusted to take care of growth. Few input data are necessary: the thickness of the annual growth rings, the annual length increase of the branches, the growth angle, the density, the longitudinal elastic modulus of wood and the acceleration of gravity.

Numerical results

The results are synthesised in a picture of the whole tree: the young branches are at the top, almost straight and pointing in the direction of growth. The old branches, at the bottom of the tree, are curved and deflected towards the bottom. Figures 1 to 3 show the influence of the growth angle. On figure 3 b, after touching the ground, the branches grow horizontally, which is not realistic.  On figure 4, the branches continue to grow with a support on the ground. Figure 6 shows a real tree where the branches have touched the ground.

The shapes of the branches are characterised by an inflection point, not predictable with the elastic criteria of MacMahon [6].

Fig
Fig 1: Small growth angle. Growth direction at 80°.

Fig

Fig. 2 : Growth direction at 45°.

Fig
Fig. 3 a : 90° growth angle without ground.

Fig

Fig. 3 b : 90° growth angle with ground.

 

Fig 4 : 45° growth angle.

Fig

 Fig. 4 : Growth direction at 45°, growth speed twice slower than on the other figures. This shape is not geometrically similar to figure 2 (scale effect).

Fig

Fig. 5 : Old tree with branches growing up again after touching earth.

Fig

Fig. 6 a : Old cypress from the Bois de Boulogne in Paris

Fig_6_b

 

Fig 6 b : Another old branch having grown after touching the ground.

Araucaria_araucana_GC  Sans_titre

Sans_titre

Fig 7 : An araucaria araucana (jardin des Poètes d'Auteuil) looking like the numerical model in 2008 and 2010 (spring and autumn).

Sans_titre

Fig. 8 Same in february 2011

2009

Fig 9 Same in september 2011

2012

Fig 10 Same in march 2012 

Araucaria

Fig 11 Same in july 2012

IMG_0508

Fig 11 Same in august 2013

Fig 12 Same in april 2014

Fig 12 Same in april 2014

Here is the original paper (C.R. Acad. Sci. Paris, t. 311,série II, pp 37-43, 1990) accompanied with the listing of the software : Branche

References

[1] R. R. ARCHER, F.E. BYRNES, On the Distribution of Tree Growth Stresses - Part I: An
Anisotropic Plane Strain Theory, Wood Sci. Technol., 8 (1974) pp. 184-196.
[2] FENG-HSIANG HSU, The influences of Mechanical loads on the form of a growing
elastic body, J. Biomech. 1 (1968) p. 303-311.
[3] J.L. NOWINSKI, Mechanics of growing materials, Int. J. Mech. Sci., 20 (1978) p. 493-
504.
[4] J. CRABBE, H. LAKHOUA, Ann. Sci. Nat., 12ème série, t. 19 (1978) p. 125.
[5] R.M. KELLOG, G.L. STEUCEK, Mechanical stimulation and xylem production of
Douglas fir, CAPPI Forest Biology Wood Chemistry Conf., Madison, 1977, pp 151-157.
[6] T. MACMAHON, Size and Shape in Biology, Science, 179 (1973), pp. 1201-1204.
[7] M. FOURNIER, Déformations de maturation, contraintes "de croissance" dans l'arbre sur
pied, réorientation et stabilité des tiges, Sém. "Développement architectural, apparition de bois
de réaction et mécanique de l'arbre sur pied", Montpellier, 20/1/1989, Ed. LMGMC, USTL
Montpellier.
[8] P.P. GILLIS, Theory of Growth Stresses, Holzforschung, Bd. 27, Heft 6 (1973) pp. 197-
207.
[9] S.P. TIMOSHENKO, Théorie de la stabilité élastique, Dunod, Paris, 1966
[10] M. FOURNIER, Mécanique de l'arbre sur pied: maturation, poids propre, contraintes
climatiques dans la tige standard, Thèse, INPL, Nancy, 1989.
[11] B. SCHAEFFER, Rhéologie des propergols en cours de polymérisation, Industrie
Minérale, N° spécial Rhéologie, T. IV (1977) N° 5, pp. 225-234.