REVIEW 3 major objections 4 minor 30 references
Orthotropic Viscoelastic Creep in Cellular Scaffolds
T0 review · 3 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read In hierarchical models of Norway spruce, the creep compliance tensor stays proportional to the elastic compliance tensor at every scale and in every tissue type, so the measured direction-dependence of creep must originate in nonlinear…
desk verdict The paper's main conclusion is baked into its input assumptions; the useful part is the hierarchical modeling framework, not the negative result about topology. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is a four-element Kelvin-Voigt chain whose creep compliance is $J_{c,AB}(t) = C^{-1}_{0,AB} \sum_{i=1}^{4} (1/\gamma_i^{\mathrm{ve}})(1-e^{-t/\tau_i})$, with fixed times $\tau_i = [0.1, 1, 10, 100]$ hours and unknown dimensionless factors $\gamma_i^{\mathrm{ve}}$ linking each viscous element to the elastic compliance. The same $\gamma_i^{\mathrm{ve}}$ set is assigned to every cell-wall layer, so each upscaling step—from layer to Y-shaped cell to tissue to growth ring—preserves the proportionality between elastic and viscoelastic behavior. The paper then extracts component-wise $\gamma_i$ from each model and shows that they collapse onto the input set, making the proportionality a property of the whole hierarchy rather than a single-scale assumption.
What would settle it
Run the same hierarchical upscaling with layer-specific $\gamma_i$ or $\tau_i$ (or with a single layer whose creep is not proportional to its elasticity) and check whether the normalized creep curves $J_{c,AB}(t)/C^{-1}_{0,AB}$ separate by direction. A separation would show that topology alone can create orthotropic creep, contradicting the paper's conclusion that nonlinear stress-concentration behavior must be the cause.
Extended reading notes
Core claim
The central result is that the effective viscoelastic compliance remains proportional to the effective elastic compliance at every hierarchical level of the model: cell wall layers, tissues, and the assembled growth ring. A single set of proportionality factors $\gamma_i^{\mathrm{ve}}$ with fixed characteristic times describes normalized creep in every direction, and this collapse is not broken by geometric nonlinearity or by replacing regular honeycomb cells with real disordered tissue images. The paper therefore asserts that the experimentally observed directional dependence of spruce creep does not emerge from tissue-scale topology. Instead, the missing directional dependence must come from nonlinear material behavior concentrated at stress concentration sites, where localized hinge-like deformation can develop.
Load-bearing premise
Everything rests on assigning every cell-wall layer the same creep spectrum, with the same proportionality factors $\gamma_i$ and the same characteristic times; if real layers differ in these parameters, tissue topology could in principle produce direction-dependent creep even in linear viscoelasticity.
Editorial extensions
If this is right
- Linear multiscale wood models that assign one scalar creep-to-elasticity factor will not reproduce experimental orthotropic creep, so long-term deflection predictions for radial and tangential directions need nonlinear cell-wall input.
- A single creep test in one direction can identify the whole linear viscoelastic spectrum of the growth ring, because all directions share the same proportionality factors in the model.
- Adding geometric nonlinearity or realistic disordered cell geometry does not change the collapsed factors, so neither effect can rescue the topological explanation.
- The next generation of wood creep models should place nonlinear constitutive laws at stress concentration sites such as cell corners, where hinge-like localized deformation can develop.
Reading between the lines
- If the collapse is generic rather than an artifact of the shared-spectrum assumption, then any measured directional spread in normalized creep is a direct experimental signal of nonlinearity, and the deviation of component-wise $\gamma_i$ from the universal set could quantify the intensity of stress-concentration mechanisms.
- The same proportionality-collapse should occur in other linear cellular composites provided all wall material has identical spectra; observing direction-dependent normalized creep in such a material would indicate nonlinearity.
- Because the input assumption of identical $\gamma_i$ and $\tau_i$ in every layer pre-structures the result, a natural extension is to repeat the study with layer-specific spectra to establish whether compositional differences between S1, S2, and S3 layers change the conclusion.
- A direct test is to measure creep of isolated earlywood tissue under radial versus tangential loading and compare component-wise normalized curves; non-coinciding curves would indicate that even a single tissue type carries nonlinear responses.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops a hierarchical multiscale finite-element model of Norway spruce, from individual cell wall layers (ML, P, S1, S2, S3) through regular Y-shaped unit cells and image-based disordered tissue RVEs up to a growth-ring RVE. Viscoelasticity is represented by a four-element Kelvin-Voigt chain in which each layer's creep compliance is proportional to its elastic compliance, with common proportionality factors gamma_i and fixed characteristic times tau_i = [0.1, 1, 10, 100] h. The gamma_i are identified inversely against macroscopic normalized creep measurements via a surrogate-based optimizer. The paper's central claim is that all simulated creep curves collapse onto a single universal set of proportionality factors regardless of tissue type, disorder, or geometric nonlinearity, and that therefore directional creep differences cannot originate from tissue topology but require nonlinear material behavior at stress concentration points.
Significance. If the central claim were established, the paper would justify a strong conclusion: that linear hierarchical models with a universal elastic-to-viscoelastic scaling cannot reproduce anisotropic creep of wood, so nonlinear cell wall behavior is necessary. The elastic part of the manuscript is a useful and largely convincing exercise: the predicted tissue and growth-ring elastic constants are compared with experimental and literature values, and the FEM framework with periodic boundary conditions is carefully described. However, the creep conclusion is not supported by the evidence. The collapse of the creep curves is a direct consequence of the constitutive assumption that every cell wall layer shares the same gamma_i and tau_i, and the paper does not test the alternative case in which different layers have different viscoelastic spectra. As submitted, the main scientific message about wood mechanics is therefore not established, although the multiscale modeling and identification pipeline could be valuable if the proportionality assumption is relaxed and the central question is revisited.
major comments (3)
- [Sec. 2.1.2, Eq. (4); Sec. 2.2.1] The model prescribes C_i^{-1} = C_0^{-1}/gamma_i and identical tau_i in every cell wall layer. Consequently each layer's compliance is C_0^{-1}(1 + sum_i gamma_i^{-1}(1 - exp(-t/tau_i))), i.e. the same scalar time function multiplies every layer's elastic compliance. For any linear elastic composite with fixed geometry and boundary conditions, scaling all constituent compliances by the same time function scales the effective compliance by that same function, independent of topology or disorder. The collapse shown in Figs. 7-9 is therefore a mathematical consequence of the input assumption, not an emergent numerical discovery. The statement in Sec. 2.2.1 that this assumption 'was numerically verified before this study using RVEs with elastic fibers in a viscoelastic matrix' refers to a different composite where fibers and matrix have different spectra and does not test the load-bearing condition that all layers share the same gamma_i and tau_i.
- [Sec. 4, Conclusions] The conclusion that 'directional differences do not emerge from topology, but must relate to non-linear material behavior at stress concentration points' is an overreach. The simulations only show that within the restricted model class in which every layer shares one gamma_i/tau_i spectrum, no directional creep emerges. If layers were assigned different spectra, for example different characteristic times for the S2 layer versus the compound middle lamella, topology could in principle produce direction-dependent effective creep even in linear viscoelasticity. The paper contains no simulations or analytical argument covering that case, so the central negative claim about topology is not supported by the presented evidence.
- [Eq. (9), Sec. 2.5] The objective function in Eq. (9) compares normalized experimental creep curves for all AB components with the simulated normalized creep. Since the simulated curves collapse by construction to one master curve, the optimization cannot extract directional information from the experiments; it effectively fits a single curve. The resulting 'universal set' gamma_i is the same set that was prescribed as input for every layer in Sec. 2.2.1, so the agreement is circular rather than an independent identification. The paper should state clearly that the inverse procedure only calibrates the common spectrum, and cannot validate the proportionality assumption against the experimental spread shown in Fig. 8b.
minor comments (4)
- [Sec. 2.1.1] The positive-definiteness condition for the 3x3 subdeterminant is written as 'nu_12 nu_23 nu_31 < 0.5(1 - ...) < 0.5', which appears malformed; the intended inequalities should be corrected.
- [Fig. 1] The label '12mm' in the tissue-scale schematic is likely intended to be '12 um'; please verify the units.
- [Sec. 3.2.1] The text alternates between 'RUC' and 'RCU' (e.g., in Fig. 7 and Sec. 3.2.1); the terminology should be harmonized throughout.
- [Sec. 3.2.1] The phrase 'proving the analogy between elastic and viscoelastic deformation mechanisms' is stronger than justified; the simulations illustrate the analogy under the assumed proportionality but do not prove it for wood.
Circularity Check
Universal creep collapse is an input, not an emergent result: identical KV spectra are prescribed in every cell-wall layer, so the conclusion that topology cannot cause directional creep is forced by construction.
-
fitted input called prediction
[Sec. 2.1.2 (Eq. 4); Sec. 2.2.1; Sec. 3.2.1 (Fig. 7)]
"The viscoelastic compliance tensor of each layer is assumed to be proportional to the corresponding elastic compliance tensor (see Sec. 2.1.1) and, even though C^{-1}_0 and C^{-1}_i are different, the ratios gamma^{ve}_i are identical for all layers. ... Remarkably, for each gamma^{ve}_i, the values collapse onto a single proportionality factor ... they match the input set gamma^{ve,opt}_i."
The input assumption is that every cell-wall layer and phase has the same gamma_i and the same tau_i = [0.1, 1, 10, 100] h. In the linear viscoelastic model used here, if every constituent compliance is a scalar multiple of its elastic compliance with one common time kernel, the effective compliance of any composite is the same scalar multiple of the effective elastic compliance, independent of topology or disorder. The inverse procedure (Eq. 9) fits gamma_i to the experimental normalized creep curves and then presents the fitted set as a discovered universal set. But the universal collapse is exactly the prescribed common spectrum propagated through the model.
-
self citation load bearing
[Sec. 2.2.1, parenthetical verification statement]
"Note that this common assumption was numerically verified before this study using RVEs with elastic fibers in a viscoelastic matrix."
The only justification offered in the text for the load-bearing identical-ratios premise is this appeal to a prior numerical verification, performed before the present study within the authors' own modeling framework. The verification described, elastic fibers in a viscoelastic matrix, does not establish that all cell-wall layers share one gamma_i/tau_i spectrum; it is a different composite morphology and does not test the excluded degree of freedom. Thus the central premise that produces the collapse is anchored in self-referential prior work rather than an independent, machine-checkable, or externally testable result. Relaxing that premise would invalidate the paper's key conclusion about topology and nonlinearity.
full rationale
The paper's central scientific claim is that directional creep differences do not emerge from tissue-scale topology and therefore require nonlinear material behavior at stress concentrations. That claim is loaded into the model's input. Section 2.1.2 states that the viscoelastic compliance tensor is proportional to the elastic compliance tensor, with the same gamma_i and the same tau_i = [0.1, 1, 10, 100] h chosen a priori. Section 2.2.1 then states that the ratios gamma_i are identical for all cell-wall layers. Under these assumptions, every constituent has the same normalized creep kernel, and any linear composite of such constituents has an effective normalized creep curve equal to that same kernel, independent of the tissue RUC/RVE geometry. The simulations in Sec. 3.2 confirm this by showing that the extracted gamma_i values 'match the input set gamma^{ve,opt}_i'. The universal collapse is therefore a mathematical consequence of the shared input spectrum, not an emergent property of wood structure. The fitted gamma_i values do contain some empirical content through the optimization against macroscopic creep data, but the conclusion that topology cannot cause anisotropic creep is not supported because the model class excludes the only mechanism that could produce such anisotropy in linear viscoelasticity, namely layer-dependent spectra. The elastic validation against literature values does not remedy this, since it concerns the skleronomous response only. The self-citation supporting the identical-ratios assumption is additionally load-bearing and not independently demonstrated. The paper is not entirely circular, because the inverse parameter identification and surrogate optimization are genuine numerical procedures, but the central interpretive conclusion is forced by the input ansatz.
Assumptions & free parameters
free parameters (2)
- γ_i^ve (i=1..4) =
[18.880, 3.386, 24.257, 2.537]
- τ_i =
[0.1, 1, 10, 100] hours
assumptions (5)
- domain assumption Viscoelastic compliance of each cell wall layer is proportional to its elastic compliance with identical γ_i and τ_i for all layers
- domain assumption Linear viscoelasticity with a Kelvin-Voigt chain of 4 elements captures wood creep
- domain assumption Orthotropy of wood and the local layer coordinate systems
- domain assumption Classical micromechanics (Halpin-Tsai, laminate theory) with literature constituent properties accurately represent cell wall layers
- domain assumption The experimental data from Maas and Wittel (2025a,b) are accurate and comparable to the model at 12% moisture
Cite this review
Pith. "Pith review of Orthotropic Viscoelastic Creep in Cellular Scaffolds." pith.science (2026). https://pith.science/paper/E4RGRWQS
@misc{pith2026250701071,
author = {Pith},
title = {Pith review of: Orthotropic Viscoelastic Creep in Cellular Scaffolds},
year = {2026},
howpublished = {\url{https://pith.science/paper/E4RGRWQS}},
note = {Machine review of arXiv:2507.01071}
}
read the original abstract
Recent measurements of Norway spruce have revealed stress-state-dependent normalized creep behavior, highlighting a gap in our fundamental understanding. This study examines whether the anisotropic response originates from the micro-structural, cellular nature of composite cell walls with varying tracheid types. Cell wall creep parameters are identified via surrogate-based inverse parameter identification, applied to hierarchical micro-mechanical and FEM models of increasing topological complexity up to the growth ring scale. Despite microstructural disorder, simulated creep curves converge toward a universal set of proportionality factors. The results indicate that directional creep behavior cannot be attributed solely to tissue-scale topology, and that realistic predictions require the inclusion of non-linear material responses at stress concentration sites.
Reference graph
Works this paper leans on
-
[1]
Physical Review E 109, 044139 (2024)
Amando de Barros, J.O., Wittel, F.K.: Unifying model for the rheological behavior of hygroresponsive materials. Physical Review E 109, 044139 (2024)
work page 2024
-
[2]
Mechanics of Time-Dependent Materials 28, 2069–2083 (2024) 21
Bengtsson, R., Bergeron, L., Afshar, R., Mousavi, M., Gamstedt, E.K.: Evaluating the viscoelastic shear properties of clear wood via off-axis compression testing and digital-image correlation. Mechanics of Time-Dependent Materials 28, 2069–2083 (2024) 21
work page 2024
-
[3]
Structural Multisiciplinary Optimization 17, 1–13 (1999)
Booker, A., Dennis Jr, J., Frank, P., Serafini, D., Torczon, V., Trosset, M.: A rig- orous framework for optimization of expensive functions by surrogates. Structural Multisiciplinary Optimization 17, 1–13 (1999)
work page 1999
-
[4]
Wood Science and Technology 12, 75–86 (1978)
Cave, I.: Modelling moisture-related mechanical properties of wood Part I: Properties of the wood constituents. Wood Science and Technology 12, 75–86 (1978)
work page 1978
-
[5]
Wood Science and Technology 12, 161––167 (1978)
Cousins, W.: Young’s modulus of hemicellulose as related to moisture content. Wood Science and Technology 12, 161––167 (1978)
work page 1978
-
[6]
Wood Science and Technology 10(1), 9–17 (1976)
Cousins, W.J.: Elastic modulus of lignin as related to moisture content. Wood Science and Technology 10(1), 9–17 (1976)
work page 1976
-
[7]
Wood Science and Technology 41(5), 443–460 (2007)
Donaldson, L.: Cellulose microfibril aggregates and their size variation with cell wall type. Wood Science and Technology 41(5), 443–460 (2007)
work page 2007
-
[8]
Holzforschung 66(8), 959–965 (2012)
Engelund, E.T., Salm´ en, L.: Tensile creep and recovery of norway spruce influenced by temperature and moisture. Holzforschung 66(8), 959–965 (2012)
work page 2012
Show all 30 references
-
[9]
Holzforschung 27, 1–7 (1973)
Fengel, V.D., Stoll, M.: ¨Uber die Ver¨ anderungen des Zellquerschnitts, der Dicke der Zellwand und der Wandschichten von Fichtenholz-Tracheiden innerhalb eines Jahrringes. Holzforschung 27, 1–7 (1973)
1973
-
[10]
Proceedings of the National Academy of Sciences 108(47), E1195–E1203 (2011)
Fernandes, A.N., Thomas, L.H., Altaner, C.M., Callow, P., Forsyth, V.T., Apperley, D.C., Kennedy, C.J., Jarvis, M.C.: Nanostructure of cellulose microfibrils in spruce wood. Proceedings of the National Academy of Sciences 108(47), E1195–E1203 (2011)
2011
-
[11]
Holzforschung 78, 624–630 (2024)
Ferrara, A., Wittel, F.K.: Micro-mechanical tests on tissue slices of Norway spruce: tensile and shear performance. Holzforschung 78, 624–630 (2024)
2024
-
[12]
Mechanics of Time-Dependent Materials 29(2), 36 (2025)
Ferrara, A., Wittel, F.K.: Tensile creep of norway spruce on the tissue scale. Mechanics of Time-Dependent Materials 29(2), 36 (2025)
2025
-
[13]
Wood Science and Technology 0(0), 18 (2025)
Ferrara, A., Amando de Barros, J.O., Koch, S.M., Wittel, F.K.: Morphogenetic vectors in the densification of cellulose scaffolds. Wood Science and Technology 0(0), 18 (2025)
2025
-
[14]
Mechanics of Time-Dependent Materials 13(4), 333–356 (2009)
Fortino, S., Mirianon, F., Toratti, T.: A 3d moisture-stress FEM analysis for time dependent problems in timber structures. Mechanics of Time-Dependent Materials 13(4), 333–356 (2009)
2009
-
[15]
Progress in Materials Science 52(8), 1263–1334 (2007)
Fratzl, P., Weinkamer, R.: Nature’s hierarchical materials. Progress in Materials Science 52(8), 1263–1334 (2007)
2007
-
[16]
Gibson, L.J., Ashby, M.F.: Cellular Solids: Structure and Properties. 2nd edn. Cambridge Solid State Science Series: Cambridge University Press (1997)
1997
-
[17]
Grimsel, M.: Mechanisches Verhalten von Holz: Struktur- und Parameteridentifikation eines anisotropen Werkstoffes. No. 1 in Dresdner Forschungen Maschinenwesen: w.e.b.-Univ.-Verl (1999)
1999
-
[18]
Gritzmann, P.: Grundlagen der Mathematischen Optimierung. Springer-Verlag (2013) Gr¨ onquist, P., Wood, D., Hassani, M.M., Wittel, F.K., Menges, A., R¨ uggeberg, M.: Analysis of hygroscopic self-shaping wood at large-scale for curved mass timber structures application. Science...
2013
-
[19]
Gutierrez-Lemini, D.: Engineering Viscoelasticity. 1st edn. Springer New York (2014)
2014
-
[20]
Hajikarimi, P., Moghadas Nejad, F.: Applications of Viscoelasticity. 1st edn. Elsevier (2021) 22
2021
-
[21]
Halpin, J., Kardos, J.L.: The Halpin-Tsai equations: a review. Polymer Engineering & Science 16(5), 344–352 (1976) Hanhij¨ arvi, A., Mackenzie-Helnwein, P.: Computational analysis of quality reduc- tion during drying of lumber due to irrecoverable deformation. I: Orthotropic v...
1976
-
[22]
Computer Methods in Applied Mechanics and Engineering 283, 1032–1060 (2015)
Hassani, M.M., Wittel, F.K., Hering, S., Herrmann, H.J.: Rheological model for wood. Computer Methods in Applied Mechanics and Engineering 283, 1032–1060 (2015)
2015
-
[23]
a review COST action e35 2004–2008: Wood machining – micromechanics and fracture
Hofstetter, K., Gamstedt, E.K.: Hierarchical modelling of microstructural effects on mechanical properties of wood. a review COST action e35 2004–2008: Wood machining – micromechanics and fracture. Holzforschung 63(2) (2009)
2009
-
[24]
Computers & structures 72(4), 459–480 (1999)
Holmberg, S., Persson, K., Petersson, H.: Nonlinear mechanical behaviour and analysis of wood and fibre materials. Computers & structures 72(4), 459–480 (1999)
1999
-
[25]
Kelly, A., Zweben, C.H.: Comprehensive composite materials (2000)
2000
-
[26]
Wood Science and Technology 53(2), 447–467 (2019) M¨ uller, J., Pich´ e, R.: Mixture surrogate models based on dempster-shafer theory for global optimization problems
Mora, D., Olaniran, S.O., R¨ uggeberg, M., Burgert, I., Herrmann, H.J., Wittel, F.K.: Mechanical behavior of chemically modified norway spruce: A generic hierarchical model for wood modifications. Wood Science and Technology 53(2), 447–467 (2019) M¨ uller, J., Pich´ e, R.: Mix...
2019
-
[27]
PhD thesis
Persson, K.: Micromechanical modelling of wood and fibre properties. PhD thesis. Lund University, (2000)
2000
-
[28]
Mechanics of Materials p 16 (2009)
Qing, H., Mishnaevsky, L.: 3D hierarchical computational model of wood as a cel- lular material with fibril reinforced, heterogeneous multiple layers. Mechanics of Materials p 16 (2009)
2009
-
[29]
International Journal of Solids and Structures 47(9), 1253–1267 (2010)
Qing, H., Mishnaevsky, L.: 3d multiscale micromechanical model of wood: From annual rings to microfibrils. International Journal of Solids and Structures 47(9), 1253–1267 (2010)
2010
-
[30]
Composites Science and Technology 72(6), 744–751 (2012) Villamizar Santamar´ ıa, S.: High-rise buildings made out of wood
Rafsanjani, A., Derome, D., Wittel, F.K., Carmeliet, J.: Computational up-scaling of anisotropic swelling and mechanical behavior of hierarchical cellular materials. Composites Science and Technology 72(6), 744–751 (2012) Villamizar Santamar´ ıa, S.: High-rise buildings made o...
2012
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.