Surgical interventions close to vulnerable structures, such as nerves, require precise handling of surgical instruments and tools. These tools not only pose the risk of mechanical damage to soft tissues, but they also generate heat, which can lead to thermal necrosis of bone or soft tissues. Researchers and engineers are trying to improve those tools through experimentation and simulations. To simulate temperature distributions in anatomical structures, reliable material constants are needed. Therefore, this study aimed at investigating the thermal conductivity of cortical and cancellous bone. Accordingly, a custom-made steady-state experimental setup was designed and validated. 6 bovine and 3 human cortical bone samples, as well as 32 bovine cancellous bone samples, with variable bone volume fraction were tested. The cancellous bone samples were scanned by micro-computed tomography (µCT) and micro-finite element (µFE) voxel models were created to calculate iteratively the thermal conductivity of the bone marrow. The experimental results provided 0.64 ± 0.04 W/ mK for bovine cortical bone and 0.68 ± 0.01 W/mK for human cortical bone. A linear dependency of thermal conductivity on bone volume fraction was found for cancellous bone [R-square (R 2 ) = 0.8096, standard error of the estimates (SEE) = 0.0355 W/mK]. The thermal conductivity of the bone marrow was estimated to be 0.42 ± 0.05 W/mK. These results will help to improve thermal finite element simulations of the human skeleton and aid the development of new surgical tools or procedures.Keywords: Thermal conductivity of compact and trabecular bone, specific heat of bone, thermal bone necrosis, temperature of cutting or drilling of bone.

Micro-finite element ([Formula: see text]FE) analyses are often used to determine the apparent mechanical properties of trabecular bone volumes. Yet, these apparent properties depend strongly on the applied boundary conditions (BCs) for the limited size of volumes that can be obtained from human bones. To attenuate the influence of the BCs, we computed the yield properties of samples loaded via a surrounding layer of trabecular bone ("embedded configuration"). Thirteen cubic volumes (10.6 mm side length) were collected from [Formula: see text]CT reconstructions of human vertebrae and femora and converted into [Formula: see text]FE models. An isotropic elasto-plastic material model was chosen for bone tissue, and nonlinear [Formula: see text]FE analyses of six uniaxial, shear, and multi-axial load cases were simulated to determine the yield properties of a subregion (5.3 mm side length) of each volume. Three BCs were tested. Kinematic uniform BCs (KUBCs: each boundary node is constrained with uniform displacements) and periodicity-compatible mixed uniform BCs (PMUBCs: each boundary node is constrained with a uniform combination of displacements and tractions mimicking the periodic BCs for an orthotropic material) were directly applied to the subregions, while the embedded configuration was achieved by applying PMUBCs on the larger volumes instead. Yield stresses and strains, and element damage at yield were finally compared across BCs. Our findings indicate that yield strains do not depend on the BCs. However, KUBCs significantly overestimate yield stresses obtained in the embedded configuration (+43.1 ± 27.9%). PMUBCs underestimate (-10.0 ± 11.2%), but not significantly, yield stresses in the embedded situation. Similarly, KUBCs lead to higher damage levels than PMUBCs (+51.0 ± 16.9%) and embedded configurations (+48.4 ± 15.0%). PMUBCs are better suited for reproducing the loading conditions in subregions of the trabecular bone and deliver a fair estimation of their effective (asymptotic) yield properties.