Shell Stability Related to Pattern Formation in Plants

[+] Author and Article Information
C. R. Steele

Division of Mechanics and Computation, Department of Mechanical Engineering, Stanford University, Stanford, CA 94305

J. Appl. Mech 67(2), 237-247 (Dec 10, 1999) (11 pages) doi:10.1115/1.1305333 History: Received December 10, 1999
Copyright © 2000 by ASME
Your Session has timed out. Please sign back in to continue.


Jean, R. V., 1994, Phyllotaxis: A Systemic Study in Plant Morphogenesis, Cambridge University Press, Cambridge, UK.
Jean, R. V., and Barabé, D., eds., 1998, Symmetry in Plants, World Scientific, Singapore.
Turing,  A. M., 1952, “The Chemical Basis of Morphogenesis,” Philos. Trans. R. Soc. London, Ser. B, B237, pp. 37–72.
Harrison, L. G., 1993, Kinetic Theory of Living Pattern, Cambridge University Press, Cambridge, UK.
Schwendener, S., 1874, Das mechanische Prinzip im anatomischen Bau der Monocotyledonen., Engelmann, Liepzig.
Karam,  G. N., and Gibson,  L. J., 1995, “Elastic Buckling of Cylindrical Shells With Elastic Cores—I, Analysis,” Int. J. Solids Struct., 32, Nos. 8–9, pp. 1259–1283.
Karam,  G. N., and Gibson,  L. J., 1995, “Elastic Buckling of Cylindrical Shells With Elastic Cores—II, Experiments,” Int. J. Solids Struct., 32, Nos. 8–9, pp. 1285–1306.
Martynov,  L. A., 1975, “A Morphogenetic Mechanism Involving Instability of Initial Form,” J. Theor. Biol., 52, pp. 471–480.
Harrison,  L. G., Snell,  J., Verdi,  R., Vogt,  D. E., Zeiss,  G. D., and Green,  B. R., 1981, “Hair Morphogenesis in Acetabularia mediterranea: Temperature-Dependent Spacing and Models of Morphogen Waves,” Protoplasma, 105, pp. 211–221.
Harrison,  L. G., and Hillier,  N. A., 1985, “Quantitative Control of Acetabularia Morphogenesis by Extracellular Calcium: A Test of Kinetic Theory,” J. Theor. Biol., 114, pp. 177–192.
Dumais,  J., and Harrison,  L. G., 2000, “Whorl Morphogenesis in the Dasycladalean Algae: The Pattern Formation Viewpoint,” Philos. Trans. R. Soc. London, Ser. B, 355, pp. 281–305.
Hernández,  L. F., 1991, “Morphometry and Surface Growth Dynamics of the Sunflower (Helianthus annuus L.) Receptacle, Its Importance in the Determination of Yield,” Suelo y Planta, 1, pp. 91–103.
Wu, C.-H., 1993, “Fourier Series for Stability of Shallow Shells With One Large Element: Application to Plant Morphogenesis,” Ph.D. thesis, Stanford University, Stanford, CA.
Dumais, J., 1999, personal communication.
Green,  P. B., 1992, “Pattern Formation in Shoots: A Likely Role for Minimal Energy Configurations of the Tunica,” Int. J. Plant Sci., 153, No. 3, pp. S59–S75.
Szilard, R., 1974, Theory and Analysis of Plates: Classical and Numerical Methods, Prentice-Hall, Engelwood Cliffs, NJ.
Green,  P. B., 1996, “Transductions to Generate Plant Form and Pattern: An Essay on Cause and Effect,” Ann. Bot. (London), 78, pp. 279–281.
Green, P. B., Steele, C. R., and Rennich, S. C., 1998, “How Plants Produce Pattern: A Review and a Proposal That Undulating Field Behavior is the Mechanism,” Symmetry in Plants, Jean, R. V., and Barabé, D., eds., World Scientific, Singapore, pp. 359–392.
Carpenter,  R., Copsey,  L., Vincent,  C., Doyle,  S., Magrath,  R., and Coen,  E., 1995, “Control of Flower Development and Phyllotaxy by Meristem Identity Genes in Antirrhinum,” Plant Cell, 7, pp. 2001–2011.
Kirchoff,  B. K., and Rutishauser,  R., 1990, “The Phyllotaxy of Costus (Costaceae),” Botanical Gaz, 151, No. 1, pp. 88–105.
Hernández,  L. F., and Palmer,  J. H., 1988, “Regeneration of the Sunflower Capitulum After Cylindrical Wounding of the Receptacle,” Am. J. Botany, 75, No. 9, pp. 1253–1261.
Kwiatkowska,  D., and Florek-Marwitz,  J., 1999, “Ontogenetic Variation of Phyllotaxis and Apex Geometry in Vegetative Shoots of Sedum Maximum (L.) Hoffm.,” Acta Soc. Botan. Polon., 68, No. 2, pp. 85–95.
Green,  P. B., Steele,  C. R., and Rennich,  S. C., 1996, “Phyllotactic Patterns: A Biophysical Mechanism for Their Origin,” Ann. Bot. (London), 77, pp. 515–527.
Green,  P. B., 1999, “Expression of Pattern in Plants: Combining Molecular and Calculus-Based Biophysical Paradigms,” Am. J. Botany, 86, No. 8, pp. 1059–1076.
Chen,  J. C. W., 1973, “The Kinetics of Tip Growth in the Nitella Rhizoid,” Plant Cell Physiol., 14, pp. 631–640.
Hejnowicz,  Z., Heinenann,  B., and Sievers,  A., 1971, “Tip Growth: Patterns of Growth Rate and Stress in the Chara Rhizoids,” Z. Pflanzenphysiol., 81, pp. 409–424.
Green,  P., 1968, “Growth Physics in Nitella: A Method for Continuous in vivo Analysis of Extensibility Based on a Micro-Manometer Technique for Turgor Pressure,” Plant Physiol., 43, No. 8, pp. 1169–1184.


Grahic Jump Location
(a) Tip growth in the unicellular alga Acetabularia. The growth takes place in the end cap, which is nearly spherical (elliptical with a/b=1.2 for this example). The cylindrical portion remains of constant diameter, equal to about 50 μm. The load-bearing wall, composed mostly of mannan polymers, can be seen as a thin transparent layer surrounding the tip. The cytoplasm (dark granular region) and the central vacuole exert a pressure of 7–10 atmospheres on the wall and thus provide the driving force for elongation of the cell. (Photograph from Dumais 14.) (b) In Acetabularia, the elongation stage stops at regular intervals and the end cap changes from nearly hemispherical to ellipsoidal (a/b=1.9 for this example). Here only the wall is shown. When the ratio of radial to axial semi-major axes of the ellipse reaches a value near 1.5, significant circumferential compression occurs, which causes buckling of the surface. (From Dumais and Harrison, 11.) This is just as in a standard thin-walled pressure vessel. (c) Axial view of Acetabularia after buckling. It appears that the compressive circumferential stress causes a buckling pattern that initiates the development of an equally spaced array of lateral hairs. (Such a symmetric pattern is called a whorl.) Subsequently, each lateral hair elongates as in (a). (From Dumais and Harrison 11.)
Grahic Jump Location
(a) Sunflower during the pattern development. The older florets are at the outer region, while the new florets are generated at the rim of the inner smooth dome. Each new floret is at the “golden section” between those of the older generation, which generates the spirals. (See Jean, 1.) At this stage, the sunflower head has a diameter of around 3–4 mm. (b) Shell model for the sunflower. The outer layer of cells (tunica) have substantially thicker walls to withstand the internal pressure, in comparison with the inner cells (corpus). The region of the floret initiation has negative Gaussian curvature of the tunica. (Geometry from Hernández 12.) (c) Stress calculated for the shell model of the sunflower. The region of floret initiation has substantial compressive stress in the circumferential direction shown by the shaded area. (From Wu 13.) (d) Effects of cuts across the sunflower capitulum. The tension regions gape open, but the region of the floret initiation is pressed together because of the compressive circumferential stress. (From Dumais 14.)
Grahic Jump Location
(a) Comparison of the actual sunflower (left) and computations (right). The computation is with the use of the von Kármán equations for a plate on an elastic foundation. The plate is compressed by a uniform radial edge force. The plate is initially flat, but edge conditions on the rotation are prescribed that have the spiral pattern. For subsequent buckling, the spiral pattern propagates toward the center, as shown. Typically, the distance between calculated buckles is about twice the natural wavelength. (Calculation by S. Rennich, color enhancement of figure by J. Dumais, from the cover of the American Journal of Botany, July, 1999.) (b) Calculation of plate post-buckling, in which spiral pattern degenerates into ridges. The distance between ridges is about equal to the natural wavelength. (Calculation by S. Rennich.)
Grahic Jump Location
Representation of a plant by an elastic layer, consisting of a sandwich plate representing the tunica, attached to an elastic half-space, representing the interior cells (corpus). The compressive force N in the tunica causes a buckling deformation, indicated by the dashed line, with the wavelength λ. The thickness of the tunica is t, the thickness of the walls of the tunica is tf, the thickness of the interior cell walls is tw and the cell diameter is L.
Grahic Jump Location
Relation of observed peak to peak distance between primordia and the tunica thickness. The dotted line is the best fit, with a slope 14.1. The dashed lines show the natural wavelength for a plate on an elastic substrate with relative volume fractions of f=10 and 100. Generally, there seems to be little difference in the whorl and spiral patterns. (Data collected from literature by C. Schmid and J. Dumais.)
Grahic Jump Location
Buckling and organ shape. (A) Wild-type stamens (st) in snapdragon. Creases traverse the formative zone, making a ring of humps. (B) The deficiens mutant of snapdragon. The formative region undulates in a plane as a ribbon, forming five cup-like primordia. (C) Out-of-plane buckling mode for a thick ring. (D) In-plane buckling of an elastic ring constrained by a wall at the outer margin. (From Green 24.)
Grahic Jump Location
Thermal strain to cause in-plane buckling as a function of the harmonic index. For a given value of geometry Rg/R, the minimum gives the critical condition. For the thicker rings, the critical condition stabilizes at the harmonic n=5.
Grahic Jump Location
(a) Initial hemispherical shell of radius a. The angle Φ is measured from the apex and the arc length S is measured from the equator. (b) Current configuration consisting of hemispherical shell with a cylindrical extension added.
Grahic Jump Location
(a) Perturbation of hemispherical shell by the radial displacement of magnitude η. (b) Elliptical end cap with semi-major and minor axes a and b.
Grahic Jump Location
Growth rate as a function of local mean stress for a simple, purely mechanical model. Growth occurs between a threshold value T and a saturation value of 1.5. For the spherical and elliptic caps, the stress is near N0, and the slope of the curve has the positive value γ. For the cylindrical region, the effective stress is higher by a factor of 1.5, and growth does not occur. The value T could represent an initial frictional resistance, while γ is the inverse of the effective viscosity.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In