Load transfer from the growing fibre into the growing medium: application to plant leaf growth/koormuse ulekanne kasvavalt kiult kasvavale keskkonnale: rakendus taime lehe kasvu kirjeldamiseks.1. INTRODUCTION The modelling of the biological growth is a promising bias of the contemporary mechanics of continuous media, which is connected with tissue engineering and manufacturing the novel types of composites for biomedical and technical applications. Biological tissues can be considered as multiphase continuous media composed of a solid porous skeleton (cellular walls, membranes, extracellular substances, conducting vessels) filled with liquids (intracellular, interstitial and delivering fluids). The delivering fluids contain dissolved substances that can be absorbed by a solid phase providing the mass increase of a growing body due to increase in the volume (surface or volumetric growth) and decrease in porosity (inner growth or remodelling) [1,2]. In growing plants cell expansion is produced by creating a driving force for water uptake by decreasing hydrostatic pressure (turgor) through stress relaxation in the cell wall [3]. Orientation of the newformed cell walls, conducting elements, and extracellular fibres is determined by the principals of the stress tensor. Cells are sensitive to the stress, and the threshold is [[sigma].sup.*] = 0.030.05 MPa [4]. In that way the biomaterials with optimal mechanical properties (maximal strength for the given external load at minimal total weight) may be produced by the growth at certain conditions. Numerous experimental observations and measurements of animal and plant tissues have revealed that tension stimulates growth in the corresponding direction while compression depresses it, and there is a linear dependence between the growth rate and the stress tensors. The Maxwell model of the viscoelastic fluid has been used for modelling the growth of the spine [5], plant roots [6] and leaves [7;8], cellular monolayers [9], and laminated tissues [10] at different loads. Comparison of the numerical calculations and experimental measurements of the growth kinematics at different conditions poses an inverse problem for determination of the rheological parameters of the Maxwell model. Many measurements for the stressfree grown plant leaves of different size, shape, and venation type show that the own growth rates [e.sup.g] ~ 0.1, depending on the growth stage and availability of water and nutrients [11]. Image analysis of dailytaken digital pictures of separate leaves have revealed some regularities between the diameters and branching angles of the conducting elements (veins) in the junctions, diameters, and lengths of the separate veins and their influence domains (the areas of the leaf blade supplied by water and nutrients by the corresponding vein) [11] that correspond to the model of the optimal transportation network of the rigid tubes with permeable walls [12]. Since the veins reinforce the leaf blade, whose mechanical strength and stability are determined by orientation, thickness, and branching angles of the main veins, the problem of optimal reinforcement of the biological tissues can be solved and the Natureinspired composite materials may be constructed on the basis of our measurement data [13]. In this paper, the important problem for the growing fibrereinforced composite, connected with load transfer from the elongating fibre into the growing surrounding tissues, is considered. 2. MATERIAL AND METHODS Plants of Abutilon hybridum, Plectrantus L., and Cissus rhombifolia were cultured in a phytotron. Dailytaken digital pictures were analysed using the computer software described in [11]. A set of Lagrange coordinates ([[xi].sup.1.sub.j], [[xi].sup.2.sub.j]) of fixed marks (vein junctions) on a growing leaf blade (j = 1 ... N, N = 200300 junctions depending on the leaf size) were chosen and their coordinates measured on the consecutive pictures using the image analysis software. The Euler coordinate system ([x.sup.1], [x.sup.2]) was connected with the direction of the main vein and the point of attachment of the leaf to the petiole (see Fig. 1). The corresponding transition law [x.sup.1,2] = [X.sup.1,2](t, [[xi].sup.1], [[xi].sup.2]) was determined in a usual way. The time series of the Euler coordinates [x.sup.1.sub.j], [x.sup.2.sub.j] of the marks j = 1 ... N was used for numerical calculations of the growth velocity field [??] and components of the strain rate tensor [e.sub.ik] = ([[nabla].sub.i][v.sub.k] + [[nabla].sub.k][v.sub.i])/2. [FIGURE 1 OMITTED] It was established that for the Abutilon hybridum leaf the azimuthal component of the growth rate was practically zero [v.sub.[theta]] ~ 0 and the radial component [v.sub.r] = a(t)r, where the function a(t) depends on the growth regime (light, nutrition, and water availability) only (Fig. 1a). Hence growth deformations of the leaf plate are the similarity transformations. In that case the anisotropic plate exhibits isotropic growth. Anisotropy of the leaf is defined by the orientation of the main, secondary, and minor veins whose growth is anisotropic and connected with the extension of the veins in the longitudinal direction without any additional cell divisions and appearance of new cell walls. Some typical examples of anisotropy in fibre orientation and distribution are presented in [1,11]. Growth of the main leaf tissue (mesophile) is determined by cell divisions in different directions and production of the isomorphic cells, which may be slightly elongated and oriented along the corresponding conducting element when they are near the vein. It means that the cell orientation and stimulation of the cell divisions in the main tissue can be determined by the load transfer from the growing vein into the mesophile. Generalization of Melan's problem of load transfer from the growing stringer into viscoelastic media with different growth rheology can be obtained in [14]. Another growth kinematics was observed for Plectrantus and Cissus leaves, where [v.sub.r] = [a.sub.r](t)r, [v.sub.[theta]] = [a.sub.[theta]](t)r[theta], and [a.sub.r] > [a.sub.[theta]] for Plectrantus leaves (growth elongation of a leaf blade, see Fig. 1b) and [a.sub.r] < [a.sub.[theta]] for Cissus leaves (growth expansion that is especially noticeable for the leaf base, see Fig. 1c). The straight lines in Fig. 1b,c correspond to the similarity transformations of the leaf blade, while the marks "plus", "rectangle", and "multiplication sign" indicate the positions of the corresponding Lagrange points on the real growing leaf surface. In that case anisotropic growth (elongation) of the veins is combined with angledependent growth of the mesophile. Theoretical aspects of anisotropic volumetric growth are discussed in [15,16]. As far as orientation of the growing veins in the leaf cannot be only determined by genetic factors at variable natural conditions (wind and rain load, gravity) and fixed in space, variation of the angles between the main veins may be determined by the dependence of the growth of the main tissues on plant sap delivery through the corresponding vessel, while the flow rate of the sap is defined by the diameter of the vein [12,13] and load transfer between the tissues or layers with different own growth rates [10,14]. In this paper the problem of the coordinated growth of the fibre and the main tissue (mesophile) by means of the stress field is investigated. 3. PROBLEM FORMULATION AND SOLUTION We consider a 2D problem for the growing leaf with pinnate venation as a rectangular plate composed of the main vein [S.sub.1] [member of] {x [member of] [0,L], y [member of] [a, a]} and the main tissue [S.sub.2] [member of] {x [member of] [0,L], y [member of] [h  a, a] [union] [a, a + h]} with different rheological properties. Due to symmetry, the upper part of the leaf is only considered (Fig. 2a). The strain tensor [??] is supposed to be a sum of the elastic [[??].sup.e] and growth strain tensors [[??].sup.g]. The tensor [[??].sup.e] obeys Hook's law [??] = [??][[??].sup.e], where [??] is the tensor of elastic moduli and [??] is the stress tensor. In agreement with experimental data D/Dt([[??].sup.g]) = [??] + [??][??], where D/Dt is time derivative, [??] is the tensor of own growth rates (at zero stress growth), and [??] is the tensor of the growth viscosities [1]. Finally we obtain the Maxwell model of the viscoelastic body with an additional source term. The momentum equation is considered in the static form due to large characteristic growth times [T.sup.*] = 12 days. Hence the problem formulation is presented by the momentum equations and constitutive equations for the two media m = 1, 2 with different rheological parameters: [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1) The boundary conditions are symmetry conditions at y = 0; zero normal stress at the perimeter; growth rate and stress continuity conditions at the interfaces y = [+ or ]a and the attachment condition to the leaf petiole at x = 0, namely: [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2) Solution of the problem (1)(2) will be found as Fourier expansion on the basis of the results for the laminated growing composites [10]. Let us introduce the Airy functions [[THETA].sup.m](x, y) = [[summation].sup.inf.sub.n=0][[theta].sup.m.sub.n](y)cos([[gamma].sub.n]x). Then the boundary conditions at x = L are satisfied and the components of the stress tensors are [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3) [FIGURE 2 OMITTED] Substitution of (3) into (1) gives the fourthorder differential equations for determination of the functions [[theta].sup.m.sub.n](y). Solution of the problem can be easily obtained in the following form: [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] Substitution of (3) into (2) gives a set of algebraic equations for determination of the eight unknown coefficients [C.sup.m.sub.14,n], m = 12, which is not presented here for brevity. Solution of the obtained system was computed for the following ranges of the material and geometrical parameters that correspond to measurement data for different plant leaves: L = 130 cm, a = 12 mm, H = 0.25 cm, [A.sub.11] ~ [A.sub.22] = 0.050.2, [A.sub.12] = [A.sub.21] = 0, [B.sub.11] ~ [B.sub.22] = 2.5 x [10.sup.6], [B.sub.12] ~ [B.sub.21] = (0.10.5)[B.sub.11], [E.sub.ik] ~ [10.sup.6][10.sup.8] Pa. 4. RESULTS AND DISCUSSIONS When [A.sub.ik] = [a.sub.ik], the problem (1)(2) has a trivial solution [[sigma].sup.1.sub.ik] = [[sigma].sup.2.sub.ik] = 0. A typical growth rate field is presented in Fig. 2b. In this case the size of the plate depends on the relation between [a.sub.11] and [a.sub.22] only. For instance, the kinematically different growth types obtained for the growing leaves and presented in Fig. 1ac can be achieved here at [a.sub.11] = [a.sub.22], [a.sub.11] > [a.sub.22], and [a.sub.11] < [a.sub.22], respectively. The relations may be connected rather with features of the growth at a cellular level (fibre orientation and distribution in the cellular wall and concentration of the phytohormones) than with the mechanical factors [3]. When the fibre exhibits the anisotropic growth and [A.sub.22] = 0 while [a.sub.11] = [a.sub.22] < [A.sub.11], the stress concentration is observed at the interface y = a only. For clarity the distribution of the stress intensity [[sigma].sup.[??]] is given in Fig. 3a at [a.sub.11] = [a.sub.22] = 0.1, [A.sub.11] = 0.2. Computer simulations for the mentioned ranges of the material parameters revealed that the region with [[sigma].sup.[??]] [not equal to] 0 is placed along the interface and is quite narrow (d [less than or equal to] a). As an illustration the growth rate field for the regions with different rheological parameters [A.sub.ik] and [B.sub.ik] is plotted in Fig. 3b. Far from the interface the growth is not affected by the stress field and is governed by the own growth rate [a.sub.ik] only. The same result has been obtained by FEM calculations for the plates with [a.sub.11] = [a.sub.22] >> [A.sub.11] (growth restriction along the fibre). In that way the stress field provides a local growth coordination within several cellular layers. When the parameters [A.sub.ik] and [a.sub.ik] differ significantly, the rectangular shape of the form is not maintained and the problem can be solved by the FEM. As far as in plant leaves the initial shape (for instance, the considered initial rectangular shape) is kept during the growth, it means that there are no significant differences in [A.sub.ik] and [a.sub.ik], though growth mechanisms may be different in veins (elongation) and in the attached mesophile (cell divisions). When a ~ h, the differences in [A.sub.11] and [a.sub.11] may lead to growth instabilities and shape variations far from the allometric law as discussed in [6,7]. This will make the subject of further investigation. [FIGURE 3 OMITTED] ACKNOWLEDGEMENTS I thank the organizers of the Euromech Colloquium 478, Prof. Juri Engelbrecht and Prof. Arkadi Berezovski, the participants in the Colloquium, and especially Prof. Gerard Maugin for stimulating discussions on the subject and a promising new insight into the problem consideration and data analysis. Received 15 January 2007 REFERENCES [1.] Belousov, L. V. and Stein, A. A. (eds). Mechanics of Growth and Morphogenesis. Modern Problems of Biomechanics, Vol. 10. Moscow University Press, Moscow, 2002. [2.] Cowin, S. C. Strain or deformation rate dependent finite growth in soft tissues. J. Biomech., 1996, 29, 647649. [3.] Schopfer, P. Biomechanics of plant growth. Amer. J. Botany, 2006, 93, 14151425. [4.] Cosgrove, D. J. Wall relaxation and the driving forces for cell expansive growth. Plant Physiol., 1987, 84, 561564. [5.] Entov, V. M. Mechanical model of scoliosis. J. Mech. Solids, 1983, 18, 199206. [6.] Stein, A. A. Deformation of the rod from the growing biological material at longitudinal compression. Appl. Math. Mech., 1998, 59, 149157. [7.] Kantor, B. Ya. and Kizilova, N. N. Mechanics of growing biological continuum. Proc. Nat. Acad. Sci. Ukraine, 2003, 2, 5660. [8.] Kantor, B. Ya. and Kizilova, N. N. Stressstrain state investigation in bidimensional growing biological material at growth restrictions. Kharkov Nat. Univ. Vestnik, 2003, 582, 107120. [9.] Stein, A. A. and Logvenkov, S. A. Space selforganization of a thin layer of the biological material growing on the substrate. Rep. Russian Acad. Sci., 1993, 328, 443446. [10.] Kizilova, N. N. and Egorova, E. S. Modelling of laminated growing biological materials. J. Mech. Eng., 2005, 56, 330342. [11.] Kizilova, N. N. Computational approach to optimal transport network construction in biomechanics. Lect. Notes in Computer Sci., 2004, 3044, 476485. [12.] Kizilova, N. N. Hydraulic properties of branching pipelines with permeable walls. Int. J. Fluid Mech. Res., 2005, 32, 98109. [13.] Kizilova, N. N. Construction principles and control over transport systems organization in biological tissues. In Physics and Control International Conference Proceedings. IEEE Computer Society, Washington, 2003, 1, 303308. [14.] Kizilova, N. N. and Kravchenko, H. P. Investigation of stressstrain state of the 2D growing media. Mech. Solid Body, 2003, 33, 158168. [15.] Maugin, G. A. and Imatani, S. Anisotropic growth of materials. J. Physique IV (Proceedings), 2003, 105, 365372. [16.] Menzel, A. Modelling of anisotropic growth in biological tissues. Biomechan. Model. Mechanobiol., 2005, 3, 147171. Natalya Kizilova Department of Theoretical Mechanics, Kharkov National University, Svobody Sq. 4, 61077 Kharkov, Ukraine; nnk_@bk.ru 

Reader Opinion