Numerical Simulations of Transverse Compression and Densification in Wood

Authors

  • John A. Nairn

Keywords:

Wood, compression, material point method, numerical modeling, plastic collapse, densification

Abstract

Numerical modeling, such as finite element analysis (FEA), of complex structures and complex materials is a useful tool for stress analysis and for failure modeling. Although FEA of wood as an anisotropic continuum is used, numerical modeling of realistic wood structures, including details of wood anatomy and variations in structure within specimens, has been beyond the capabilities of FEA and other methods. In contrast, the recently derived material point method (MPM) has features that make it amenable to analysis of realistic wood structures. To demonstrate the capabilities of MPM, simulations were done for wood in transverse compression. Some advantages of MPM are that it is easy to discretize micrographs of wood specimens into a numerical model, it can handle large deformations, it can model elastic-plastic cell-wall properties, and it automatically accounts for contact between cell walls. MPM simulations were run for softwood and hardwood loaded in either radial or tangential compression. The simulations reproduced many features of wood compression, gave insight into effects of wood anatomy on compression, and may be the first numerical calculations of realistic wood structures extended through to full densification without numerical difficulties.

References

Atluri, S. N., and S. P. Shen. 2002. The meshless local Petrov-Galerkin (MLPG) method. Tech. Science Press, Forsyth, GA.nBardenhagen, S. G., and E. M. Kober. 2004. The generalized interpolation material point method. Computer Modeling Eng. & Sci.B:477-496.nBardenhagen, S. G., J. U. Brackbill, and D. Sulsky. 2000. The material point method for granular materials. Computer Methods Appl. Mech. Eng.187:529-541.nBardenhagen, S. G., J. E. Guilkey, K. M. Roessig, J. U. Brackbill, W. M. Witzel, and J. C. Foster. 2001. An improved contact algorithm for the material point method and application to stress propagation in granular materials. Computer Modeling Eng. Sci.2:509-522.nBardenhagen, S. G., A. D. Brydon, and J. E. Guilkey. 2005. Insight into the physics of foam densification via numerical simulation. J. Mech. Phys. Solids.53:597-617.nBelytschko, T., Y. Y. Lu, and L. Gu. 1994. Element-free Galerkin methods. Int. J. Num. Meth. Eng.37:229-256.nBodig, J. 1963. The peculiarity of compression of conifers in radial direction. Forest Prod. J., 13:438.nBodig, J. 1965. The effect of anatomy on the initial stressstrain relationship in transverse compression. Forest Prod. J.15:197-202.nBodig, J. 1966. Stress-strain relationship for wood in transverse compression. J. Materials1:645-666.nBodig, J., and B. A. Jayne. 1982. Mechanics of wood and wood composites. Van Nostran-Reinhold Co, Inc., New York, NY.nBourke, P. 2004. BMP file format. http://astronomy.swin.edu.au/~pbourke/dataformats/bmp/'>http://astronomy.swin.edu.au/~pbourke/dataformats/bmp/nBrydon, A., and S. Bardenhagen. 2005. Personal communication.nBrydon, A., and S. Bardenhagen., E. A. Miller, and G. T. Seidler. 2005. Simulation of the densification of real open cell foam microstructures. J. Mech. Phys. Solids.53:2638-2660.nDavids, W. G., E. N. Landis, and S. Vasic. 2003. Lattice models for the prediction of load-induced failure and damage in wood. Wood Fiber Sci.35:120-134.nEasterling, K. E., R. Harrysson, L. J. Gibson, and M. F. Ashby. 1982. On the mechanics of balsa and other woods. Proc. Royal Soc. Lond. A.383:31-41.nGibson, L. J., and M. F. Ashby. 1997. Cellular solids: Structure and properties. Cambridge University Press, UK.nGibson, L. J., K. E. Easterling, and M. F. Ashby. 1981. The structure and mechanics of cork. Proc. Royal Soc. Lond. A.377:99-117.nGibson, L. J., M. F. Ashby, G. S. Schajer, and C. I. Robertson. 1982. The mechanics of two-dimensional cellular materials. Proc. Royal Soc. Lond. A.382:25-42.nGuilkey, J. E., and J. A. Weiss. 2003. Implicit time integration for the material point method: quantitative and algorithmic comparisons with the finite element method. Int. J. Num. Meth. Eng.57:1323-1338.nGuilkey, J. E., T. Harmon, B. Kashiwa, A. Xia, and P. A. McMurtry. 2003. An Eulerian-Lagrangian approach for large deformation fluid-structure interaction problems, part 1: algorithm development. Proc. 2nd Inter. Conf. on Fluid-Structure Interactions. June 24-26, 2003, Cadiz, Spain.nGuo, Y., and J. A. Nairn. 2004. Calculation of J-Integral and stress intensity factors using the material point method. Computer Modeling Eng. Sci.6:295-308.nHaygreen, J. G., and J. L. Bowyer. 1996. Forest products and wood science: An Introduction. Iowa State University Press, Ames, IA.nJones, R. M. 1975. Mechanics of composite materials. Scripta Book Company, Washington, DC.nKennedy, R. W. 1968. Wood in transverse compression. Forest Prod. J.18:36-40.nKolsky, H. 1963. Stress waves in solids. Dover, New York, NY.nKultikova, E. V. 1999. Structure and properties relationships of densified wood. M.S. Thesis, Wood Science and Forest Products, Virginia Polytechnic Institute and State University, Blacksburg, VA.nKunesh, R. H. 1968. Strength and elastic properties of wood in transverse compression. Forest Prod. J.18:36-40.nLandis, E. N., S. Vasic, W. G. Davids, and P. Parrod. 2002. Coupled experiments and simulations of microstructural damage in wood. Exp. Mech.42:389-394.nNairn, J. A. 2003. Material point method calculations with explicit cracks. Computer Modeling Eng. Sci.4:649-664.nNairn, J. A. 2005a. Source code and documentation for NairnMPM Code for Material Point Calculations http://oregonstate.edu/~nairnj/'>http://oregonstate.edu/~nairnj/nNairn, J. A. 2005b. Diffusion calculations in NairnMPM solid mechanics code. Unpublished results, see Nairn, 2005a for source code.nParker, S. G. 2002. A component-based architecture for parallel multi-physics PDE simulation. International Conference on Computational Science (ICCS2002) Workshop on PDE Software, April 21-24, 2002, Amsterdam, Netherlands.nPellicane, P. J., J. Bodig, and A. L. Mrema. 1994a. Behavior of wood in transverse compression. J. Testing Eval.22:383-387.nPellicane, P. J., J. Bodig, and A. L. Mrema. 1994b. Modeling wood in transverse compression. J. Testing Eval.22:376-382.nShiari, B., and P. M. Wild. 2004. Finite element analysis of individual wood-pulp fibers subjected to transverse compression. Wood Fiber Sci.36:135-142.nSmith, R. A., M. J. Paulus, J. M. Branning, and P. J. Phillips. 2001. X-ray computed tomography on a cellular polysiloxane under compression. J. Cellular Plastics37: 231-248.nSmith, I., E. Landis, and E. Gong. 2003. Fracture and fatigue in wood. John Wiley & Sons, Ltd., New York, NY.nSulsky, D., and H. K. Schreyer. 1996. Axisymmetric form of the material point method with applications to upsetting and Taylor impact problems. Computer Methods. Appl. Mech. Eng.139:409-429.nSulsky, D., Z. Chen, and H. L. Schreyer. 1994. A particle method for history-dependent materials. Computer Methods Appl. Mech. Eng.118:179-186.nSulsky, D., S.-J. Zhou, and H. L. Schreyer. 1995. Application of a particle-in-cell method to solid mechanics. Comput. Phys. Commun.87:236-252.nTabarsa, T. 1999. Compression perpendicular-to-grain behaviour of wood. Ph.D. Thesis, Forestry and Environmental Management, The University of New Brunswick, Canada.nTabarsa, T., and Y. H. Chui. 2000. Stress-strain response of wood under radial compression. Part I. Test method and influences of cellular properties. Wood Fiber Sci.32: 144-152.nTabarsa, T., and Y. H. Chui. 2001. Characterizing microscopic behavior of wood under transverse compression. Part II. Effect of species and loading direction. Wood Fiber Sci.33:223-232.nZhou, S. 1998. The numerical prediction of material failure based on the material point method. Ph.D. Thesis, University of Mexico, Albuquerque, NM.nZhu, H. X., N. J. Mills, and J. F. Knott. 1997. Analysis of the Higher Strain Compression of Open Cell Foams. J. Mech. Phys. Solids45:1875-1904.n

Downloads

Published

2007-06-05

Issue

Section

Research Contributions