A generalized anisotropic quadric yield criterion and its application to bone tissue at multiple length scales

Nonlinear computational analysis of materials showing elasto-plasticity or damage relies on knowledge of their yield behavior and strengths under complex stress states. In this work, a generalized anisotropic quadric yield criterion is proposed that is homogeneous of degree one and takes a convex quadric shape with a smooth transition from ellipsoidal to cylindrical or conical surfaces. If in the case of material identification, the shape of the yield function is not known a priori, a minimization using the quadric criterion will result in the optimal shape among the convex quadrics. The convexity limits of the criterion and the transition points between the different shapes are identified. Several special cases of the criterion for distinct material symmetries such as isotropy, cubic symmetry, fabric-based orthotropy and general orthotropy are presented and discussed. The generality of the formulation is demonstrated by showing its degeneration to several classical yield surfaces like the von Mises, Drucker–Prager, Tsai–Wu, Liu, generalized Hill and classical Hill criteria under appropriate conditions. Applicability of the formulation for micromechanical analyses was shown by transformation of a criterion for porous cohesive-frictional materials by Maghous et al. In order to demonstrate the advantages of the generalized formulation, bone is chosen as an example material, since it features yield envelopes with different shapes depending on the considered length scale. A fabric- and density-based quadric criterion for the description of homogenized material behavior of trabecular bone is identified from uniaxial, multiaxial and torsional experimental data. Also, a fabric- and density-based Tsai–Wu yield criterion for homogenized trabecular bone from in silico data is converted to an equivalent quadric criterion by introduction of a transformation of the interaction parameters. Finally, a quadric yield criterion for lamellar bone at the microscale is identified from a nanoindentation study reported in the literature, thus demonstrating the applicability of the generalized formulation to the description of the yield envelope of bone at multiple length scales.

Effect of boundary conditions on yield properties of human femoral trabecular bone

Trabecular bone plays an important mechanical role in bone fractures and implant stability. Homogenized nonlinear finite element (FE) analysis of whole bones can deliver improved fracture risk and implant loosening assessment. Such simulations require the knowledge of mechanical properties such as an appropriate yield behavior and criterion for trabecular bone. Identification of a complete yield surface is extremely difficult experimentally but can be achieved in silico by using micro-FE analysis on cubical trabecular volume elements. Nevertheless, the influence of the boundary conditions (BCs), which are applied to such volume elements, on the obtained yield properties remains unknown. Therefore, this study compared homogenized yield properties along 17 load cases of 126 human femoral trabecular cubic specimens computed with classical kinematic uniform BCs (KUBCs) and a new set of mixed uniform BCs, namely periodicity-compatible mixed uniform BCs (PMUBCs). In stress space, PMUBCs lead to 7–72 % lower yield stresses compared to KUBCs. The yield surfaces obtained with both KUBCs and PMUBCs demonstrate a pressure-sensitive ellipsoidal shape. A volume fraction and fabric-based quadric yield function successfully fitted the yield surfaces of both BCs with a correlation coefficient R2≥0.93. As expected, yield strains show only a weak dependency on bone volume fraction and fabric. The role of the two BCs in homogenized FE analysis of whole bones will need to be investigated and validated with experimental results at the whole bone level in future studies.

Relation between Excitation Power Density and Er3+ Doping Yielding the Highest Absolute Upconversion Quantum Yield

The upconversion quantum yield (UCQY) is one of the most significant parameters for upconverter materials. A high UCQY is essential for a succesful integration of upconversion in many applications, such as harvesting of the solar radiation. However, little is known about which doping level of the rare-earth ions yields the highest UCQY in the different host lattices and what are the underlying causes. Here, we investigate which Er3+ doping yields the highest UCQY in the host lattices β-NaYF4 and Gd2O2S under 4I15/2 → 4I13/2 excitation. We show for both host lattices that the optimum Er3+ doping is not fixed and it actually decreases as the irradiance of the excitation increases. To find the optimum Er3+ doping for a given irradiance, we determined the peak position of the internal UCQY as a function of the average Er−Er distance. For this purpose, we used a fit on experimental data, where the average Er−Er distance was calculated from the Er3+ doping of the upconverter samples and the lattice parameters of the host materials. We observe optimum average Er−Er distances for the host lattices β-NaYF4 and Gd2O2S with differences <14% at the same irradiance levels, whereas the optimum Er3+ doping are around 2× higher for β-NaYF4 than for Gd2O2S. Estimations by extrapolation to higher irradiances indicate that the optimum average Er−Er distance converges to values around 0.88 and 0.83 nm for β-NaYF4 and Gd2O2S, respectively. Our findings point to a fundamental relationship and focusing on the average distance between the active rare-earth ions might be a very efficient way to optimize the doping of rare-earth ions with regard to the highest achievable UCQY.