Xin Yang1, Wenjun Ma1, Wenping Gu2, Zhaoyang Zhang1, Ben Wang1, Yan Wang3, Shifeng Liu3. 1. College of Materials Science and Engineering, Xi'an University of Technology Xi'an 710048 China yangx@xaut.edu.cn. 2. Electronic Science and Technology, Chang'an University Xi'an 710061 China. 3. School of Metallurgical and Engineering, Xi'an University of Architecture & Technology Xi'an 710055 China Liushifeng66@xauat.edu.cn.
Titanium alloy has advantages of high specific strength, low elastic modulus, high energy absorption and strong corrosion resistance. It has been widely used in the field of medical implantation especially for its good biocompatibility, among which Ti6Al4V (Ti64) alloy accounts for more than 50%.[1] However, the elastic modulus of the bulk Ti64 alloy ranges from 90 GPa to 110 GPa[2] which is much higher than that of cortical bones such as the femur and tibia (0–35 GPa),[3] and such large elastic modulus mismatch between the implant and the surrounding bone tissue results in stress shielding around the implant, leading to failure of the regeneration of bone cells.[4] Therefore, it is of great importance to reduce modulus mismatch in the biomedical field.[5] The Ti64 lattice structure has a low elastic modulus (0–40 GPa) that matches well with that of the femur and tibia.[6] Additionally, the space between the structures promotes the integration of bone cells making it the preferred choice for medical implantation.[7] In recent years, selective laser melting (SLM) technology is capable of realizing complex lattice patient-specific structures from medical grade metallic powder materials.[8]At present, the research on SLM-Ti64 lattice structure mainly focuses on the preparation of Ti64 lattice structure and its compression performance. M. Mazur et al.[9] took the Ti64 spherical powder with an average particle size D50 of 40 μm as the raw material, and studied the manufacturability of the Ti64 alloy lattice strut using SLM technology. The results show that the minimum manufacturing strut diameter is 0.2 mm, and the minimum inclination angle is 20°. S. Gangireddy et al.[10] used Ti64 spherical powder with a particle size distribution of 15–45 μm as the raw material, and prepared the octahedral lattice structure using SLM technology. The results show that total energy absorption and compressive strength are positively correlated with relative density. When the relative density is 0.354, the compressive strength is 75 MPa. Other studies also showed that the Ti64 lattice structure prepared by SLM usually has higher compressive strength, lower elastic modulus and good biocompatibility than traditional structure.[11]Clinical results indicate that most cortical bone fractures are caused by tensile or bending strains.[12] Although only a limited number of studies have so far been conducted on the tensile properties of the SLM-ed Ti64 lattice structure, it is however possible to come up with a few. Among them, S. Raghavendra[13]et al. used SLM technology to prepare a random porous lattice structure of Ti64 alloy. When the relative density is equivalent to that of cortical bone, the tensile strength can reach 250 MPa, but its disordered microstructure is more prone to failure under tensile load than under compression. S. M. Ahmadi et al.[14] treated the diamond Ti64 lattice structure prepared by SLM technology at 1050 °C for 2 hours. The results show that the microstructure changes from fine needle-shaped α′ martensite to α + β mixed structure after heat treatment as previously observed,[15] the α + β interface is an effective barrier for crack propagation. The growing α-laths make its tensile strength drop sharply. In order to be used safely on the human body for a long time, it is urgent to obtain high tensile strength while improving its microstructure. In this study, F2CCZ (face-centered cubic unit cell with vertical struts) structure type tensile parts were prepared by SLM. Double solid solution aging heat treatment is used to homogenize the alloy microstructure while obtaining a controllable mixed structure, improving the microstructure while obtaining high tensile strength. On this basis, the tensile performance of the samples both with and without heat treatment was studied and discussed.
Experimental
A face center cubic unit cell with vertical struts (F2CCZ) lattice structure was designed by Solid Works 2013. The size of the unit cell was magnified to 1 mm × 1 mm × 1 mm. The number of unit cells was set at 3 × 3 × 20. Struts diameters of lattice structures were set from 0.2 mm to 0.33 mm. It is named F-X according to its length–diameter ratio. For example, a unit cell with a length of 1 mm and a diameter of 0.2 mm is named F-5, and the specific strut dimensions are shown in Table 1. The designed tensile sample with lattice structure and a unit cell are shown in Fig. 1.
Dimension of the struts
Number
l/mm
d/mm
l/d
F-5
1
0.2
5
F-4
1
0.25
4
F-3
1
0.33
3
Fig. 1
Schematic diagram of the tensile sample with lattice structure and a unit cell.
The lattice structures were manufactured by a selective laser melting equipment (Mlab cusing R, Concept Laser, Germany). The parameters for manufacturing samples are shown in Table 2. Argon atmosphere was used to prevent oxidation of Ti64 alloys during melting. The Ti64 raw powder morphology and particle size distribution are shown in Fig. 2a and b. The scanning strategy is a 45° checkerboard and the scanning direction is rotated by 90° after each layer, as shown in Fig. 2c. The powder exhibits high sphericity and a small amount of satellite powder. The particle size is 15–53 μm, of which D50 is 45 μm. Subsequently, a double solution-aging heat treatment process was carried out: firstly, the sample was kept at 1030 °C for 120 min, then at 1000 °C for 60 min, and finally at 600 °C for 360 min. All alloys were air-cooled.
The process parameters for SLM
Item
Parameter
Laser power
100 W
Scanning speed
600 mm s−1
Layer thickness
25 μm
Scanning space
80 μm
Fig. 2
(a) Morphology of Ti64 raw powder, (b) particle size distribution and (c) checkerboard scan strategy.
The samples were sanded with different mesh sandpaper and polished with 0.5 μm SiO2 suspension. After polishing, the samples were etched with a mixed solution (92 ml H2O + 6 ml HNO3 + 2 ml HF) for 90 s. The phase compositions of the raw powder and the shaped sample were analyzed by an XRD X-ray diffractometer (XRD-7000, Japan), and the XRD test was performed using a Cu target in a step of 0.02°. The surface morphology of the lattice structure was observed by a laser scanning confocal microscope (LEXT OLS4000, Japan). Samples for TEM (Themis Z, America) observations were prepared by electrolytic polishing using a mixture of 5% perchloric acid, 25% normal butanol and 70% methanol at 243 K and 28 V. The fracture morphology and microstructure were observed with scanning electron microscopy (Gemini SEM 300, Germany).Tensile tests were carried out on the specimens in a universal testing machine (INSTRON 5982, United States) with a 100 kN force cell. Tensile tests were conducted at a constant displacement (0.2 mm min−1). The stress is calculated by the following formula according to ISO 13314:2011, Sdesign represents the design area perpendicular to the loading direction.Since the prepared lattice structure has a regular pore structure, the density of the porous tensile specimen can be calculated by the following formula:where V is the volume of the specimen (cm3), and M is the mass of the specimen (g). Both do not consider the clamping part of the lattice tensile sample. The average porosity of the lattice sample can be calculated by the mass–volume method, as shown in the following formula:where ρlattice is the density of the lattice tensile sample (g cm−3), ρsolid is the density of fully dense Ti64 alloy (4.43 g cm−3), V is the volume of the specimen (cm3) and M is the mass of the specimen (g).
Results
Profiles of lattice structures
The profiles of lattice structures with different aspect rations are presented in Fig. 3. Some hanging slags were found on the inclined struts of F2CCZ structures. This phenomenon is often observed in the lattice structures due to the force of gravity on the melting pools.[16] The analysis of the strut diameters of lattice structures with different aspect ratios can be seen in Table 3. As the designed strut diameter increases from 0.2 mm to 0.33 mm, its actual dimensional accuracy improves. The diameter of inclined-strut is much larger than the design model, and the diameter of build direction-strut is close to the design model. The spatter particles and spherical particles on the surface of lattice structure influenced the dimensional accuracy.[17] The path planning and spot compensation of laser also affect the shape of lattice structures.[18]
Fig. 3
Profiles of lattice structure: (a) SLM-ed, (b) F-5-struts, (c) F-4-struts and (d) F-3-struts.
Analysis of dimensional accuracy
Number
Inclined-measure (mm)
Build direction-measure (mm)
Struts-design (mm)
F-5
0.31 ± 0.031
0.19 ± 0.024
0.2
F-4
0.35 ± 0.022
0.27 ± 0.003
0.25
F-3
0.36 ± 0.013
0.33 ± 0.005
0.33
XRD analyses
The XRD patterns of SLM-ed and HT-ed Ti64 alloys are shown in Fig. 4. XRD patterns of HT-ed contain a mixture of α phase and β phase. However, the XRD patterns of SLM-ed only contain α′ phase. When the cooling rate exceeds 410 °C s−1, complete α′ martensite will be produced, and this transformation is very common in the Ti64 alloy manufactured by SLM.[19] Because the cooling speed is very fast (103−105 °C s−1) during the SLM processing, most of β structures transformed to the martensitic structures. Moreover, these peaks of HT-ed show higher intensity of α characteristic peaks than SLM-ed. It implies that the size of α phase grain grew up during the heat treatment. The quantitative analyses of XRD patterns are shown in Table 4. It shows that the sample contains 5.1% β phase after heat treatment, while the SLM-ed alloy has almost no β-phase. The similar results were found by Ahmadi et al.[14] and Yuan et al.[20] Different from the SLM-ed bulk Ti64 alloy, no β phase is observed in the SLM-ed lattice structure. This may be resulted from the faster cooling rate in the SLM-ed lattice material than in the SLM-ed bulk material.
Fig. 4
XRD patterns of Ti6Al4V alloy in different states.
Quantitative analysis of XRD patterns
SLM-ed
HT-ed
α-Phase (%)
Close to 100
94.9
β-Phase (%)
Close to 0
5.1
Microstructural observation
Microstructures of SLM-ed and HT-ed alloys are shown in Fig. 5. Fig. 5a and b shows the longitudinal section, with columnar grains parallel to the building direction. These elongated structures arose from the prior β columnar grains formed during the solidification stage of the liquid and were retained to room temperature.[21] Moreover, the cooling rate in the SLM process is as high as 105 °C s−1.[21] At this cooling rate, β phase will directly transform to α′ martensite phase with an acicular appearance, and the acicular α′ martensite was distributed within the columnar structure. The acicular α′ martensitic structure occupies most of the SLM-ed, which is different from the ASTM Grade 5 alloy (α + β) and the SLM bulk Ti64 alloy (α'+ β). Comparing SLM-ed with HT-ed (Fig. 5c and d), it is found that the α laths in HT-ed became coarser, which was also reported in ref. 22. The size of α lath increases with temperature. In addition, some nanosized particles are dispersedly distributed on the α laths. It was demonstrated that these nano-sized particles are the nano-sized β particles in the research.[19] Some pores are observed in this topography, and these pores will affect the mechanical properties[23] and corrosion performance.[24]
Fig. 5
OM images (a) and (c) of SLM-ed and HT-ed, respectively, and the corresponding SEM images (b) and (d).
Fig. 6a and b show the bright field images and corresponding selected area electron diffraction (SAED) patterns of SLM-ed alloy. It is seen that the SLM-ed alloy is mainly composed of α′ phase, which is caused by the non-diffusion phase transition with high cooling rate. There are high dislocation density and many crystal defects in the martensite grains, indicating that the stress relaxation in the transition from β to α′ phase is mainly realized by dislocation formation. No residual β-phase diffraction pattern was detected in the SLM-ed alloy, which is consistent with the XRD results. The similar result was observed by Wu et al.[25]
Fig. 6
TEM bright field images (a and d) of SLM-ed and HT-ed, respectively, and the corresponding diffraction pattern (b and e), and HRTEM images (c and f).
Fig. 6d and e show the bright field image and corresponding SAED patterns of the HT-ed alloy, in which the α and β phases are observed. When solid solution at high temperature, α phase transforms to β phase. Diffusion transition occurs during air cooling, and β phase forms α lamellar and remaining β phase. The similar result was found by Zhang et al.[6] Different from the SLM-ed alloy, no α′ phase is observed in the HT-ed alloy. This resulted from the slower cooling rate in the HT-ed alloy than in the SLM-ed alloy. The orientation relationship between the β and α phases is presented by: [113]β//[1−210]α, as reported by Zhang et al.[26] and Wan et al.[27] The high-resolution TEM image presented the interface between a BCC and HCP structures (Fig. 6f), further confirmed the orientation relationship. According to the requirement for implants materials specified in ISO 20160, a homogeneous microstructure of high-density biomaterials is necessary to provide material integrity and stability.[28] According to Fig. 5 and 6, the width of acicular martensite grains of the SLM-ed alloy is in the range of 30–450 nm. The heat treatment not only eliminates the high dislocation density of the SLM-ed alloys, but also obtains martensite laths with multi-level scales. The width of the primary α produced by the first solid solution is greater than 2 μm, the width of the secondary α produced by the second solid solution is in the range of 1–2 μm, and the width of α produced by the final aging is about 450 nm. N. Hrabe et al.[29] showed that the width of α lath is directly related to its mechanical properties. The existence of multi-scale α-laths can obtain better comprehensive mechanical properties.To further identify the effect of heat treatment on the SLM-ed Ti64 lattice structure microstructure, the samples with and without heat treatment were analyzed by electron backscatter diffraction (EBSD). The results are summarized in Fig. 7. As shown in Fig. 7a–d, the grain diameters of the samples without heat treatment (SLM-ed) are all less than 10 μm, and after heat treatment, the grain diameters show a cross-length scale distribution (small size aging α and coarse size primary α exists at the same time). Further confirmed the results of SEM and TEM. With heat treatment, the maximum MUD value is reduced from 23.67 to 14.35 (Fig. 7e and f), substantially weakening the texture of the α phase. A special heat treatment system is used to homogenize the microstructure while obtaining a controllable mixed structure, thereby obtaining excellent tensile properties, which lays the foundation for the further extensive application of lattice structures.
Fig. 7
The change of grain size and texture of Ti–6Al–4V produced by SLM by dual solution-aging heat treatment: Inverse pole figure maps along the build direction (z) for the α phase (measured by EBSD) in samples without (a) and with (d) heat treatment, histograms of the grain size of samples without (b) and with (e) heat treatment and {0001} contoured pole figures (in MUD: multiples of uniform distribution) of the measured α phase in samples without (c) and with (f) heat treatment.
Mechanical properties
Tensile curves of the lattice structures are shown in Fig. 8. The modulus and tensile strength of the alloy after heat treatment both increased compared with the SLM-ed state. Additionally, the tensile curves of HT-ed alloys show smoother profiles than those of the SLM-ed alloys. Due to the existence of residual stress and non-uniform microstructure, the tensile properties of the SLM-ed alloys are dispersive, especially the elongation shows larger deviations. The microstructure of the HT-ed alloy is uniform, and the tensile properties are more stable and predictable. As shown in the inserted picture, the fracture of the SLM-ed alloy occurred at the junction of the clamping part and the gauge length. The simulation results of SLM-ed alloys by ANSYS 19.2 (this version has an additive manufacturing material module) also show that there is serious stress concentration at the transition between lattice and solid parts (Fig. 8a). Because the high density of dislocations and residual stress cause severe stress concentration at the connection. The fracture of the HT-ed alloy occurred at the center of the gauge length, which means that the heat treatment eliminated the stress concentration at the junction of the lattice and the solid parts (Fig. 8b). Studies have shown that the high-density dislocations in the SLM-ed bulk sample have a certain strengthening effect.[30] Unlike the SLM-ed bulk sample, the high density of dislocations in the SLM-ed lattice structure causes serious stress concentration and causes premature failure. Multi-level α-laths and nano-β particles make HT-ed alloys have excellent tensile properties (Fig. 8b). Z. Wang et al.[31] found that the use of SLM to form Co–Cr–Mo alloys also produces microstructures across different length scales. The multi-scale microstructure gives it excellent comprehensive mechanical properties.
Fig. 8
Tensile curves and fracture morphologies of Ti64 lattice structures with different states: (a and c) of SLM-ed, (b and d) of HT-ed.
The fracture morphologies of lattice structures are shown in Fig. 8c and d. The fracture morphology of SLM-ed consists of a small amount of shallow dimples and large smooth plane (Fig. 8c), which means that the main fracture mechanism is brittle fracture. The fracture morphology of HT-ed (Fig. 8d) shows similar fracture characteristics to the SLM-ed alloy, but the depth of dimples became deeper and there are a lot of cleavage planes. These results demonstrate that the fracture mechanism of HT-ed is quasi-cleavage fracture. This also proves that the multi-scale α-laths and nano-β particles make the HT-ed alloy have excellent comprehensive mechanical properties (high tensile strength and appropriate ductility).
Discussion
The results in present study indicate that the size of the SLM-ed lattice structure and heat treatment process affects the microstructure and phase composition of the Ti64 alloy, and play an important role in the mechanical properties of the lattice structure. Different from SLM-ed bulk alloys, the SLM-ed lattice structure samples are almost all acicular α′ martensite, and there is no residual β phase. This is due to the smaller size of the lattice structure resulted from faster cooling (>106 K s−1). The double solution-aging heat treatment reduces the dislocations and crystal defects of the SLM-ed alloy while obtaining a controllable mixed structure, ensuring the integrity and stability of the implant.As shown in Fig. 9a and b, the elastic modulus and the tensile strength of the HT-ed alloy and the SLM-ed alloy decreases with the increase of porosity. The heat treatment makes α′ martensite transform to α phase along the α′ phase grain boundary, V atom diffuses from the inside of α to the grain boundary, and β phase is generated along the V-rich region. Therefore, α + β lamellar structure is formed after solution treatment. After the aging treatment, the remaining β phase is divided into small strip-shaped aging α. The heat treatment not only eliminates the residual stress of the SLM-ed alloy, but also precipitates nano-β particles and small strips of aging α. The nano-sized β particles also were found in the Ti64 alloy fabricated by spark plasma sintering after aging at 450 °C.[32] These factors make the tensile strength and elastic modulus of HT-ed alloy higher than that of the SLM-ed one. In addition, the numerical fluctuation of Young's modulus and tensile strength of HT-ed alloys is smaller, showing that the HT-ed microstructure is homogenized and the residual stress is eliminated.
Fig. 9
(a) Tensile strength and (b) Young's modulus of the SLM-ed and HT-ed versus the density; (c) tensile strength and (d) Young's modulus versus density for various lattice structure, including our work, cortical bone, and bending-dominated lattice material.[33–36]
Fig. 9c and d summarizes the tensile strength of human bone,[33] bending-dominated lattice structure (foam),[34] and stretch-dominated lattice structure (F2CCZ). Similar to the compression behavior of the lattice structure,[9] the tensile strength and Young's modulus of the stretch-dominant lattice structure are higher than those of the bending-dominant lattice structure. This is due to the higher load-bearing capacity and higher structural efficiency. Compared with human bone, the tensile strength of SLM-ed lattice structure Ti64 is much higher than that of human bone at the same density, and the Young's modulus is also close to that of human bone. The results show that the special heat treatment can homogenize the microstructure and reduce the residual stress while also improving its tensile strength.A number of predictive models have been proposed to relate void topology and material property parameters to mechanical properties. The most noted one is the Gibson–Ashby model,[37] described by the following equation for different cellular solids.where Elattice, ρlattice and σlattice are the Young's modulus, density and tensile strength of lattice structure, and Esolid, ρsolid and σsolid are the Young's modulus, density and tensile strength of fully dense Ti64. In this study, Esolid is equal to 110 GPa,[21]σsolid is equal to 1100 MPa.[38] The constants of C1, C2, n1 and n2 are depended on given lattice structure. The value of index n is closely related to the lattice structure. The better the structural stability and the higher the material utilization, the smaller the value of n. C and n were calculated using Gibson–Ashby equation and the average porosity is 77.3%, 67.5% and 48.5%, respectively.Fig. 10 plots the ultimate strength vs. density and modulus vs. density on the logarithmic scale for both SLM-ed and HT-ed alloys using the data listed in Fig. 9. The mechanical properties of the SLM-ed alloys are relatively dispersed, and the R2 value (0.536–0.847) is small. After heat treatment, the sample obtained a high R2 value (0.893–0.975), indicating that the Gibson–Ashby model can be used to describe the SLM Ti64 lattice after the residual stress is removed. This observation also shows that the value of R2 after heat treatment is close to 1, which means the data is more stable and reliable.
Fig. 10
Log–log plots on UTS and modulus vs. density for SLM-ed (a and c) and HT-ed (b and d).
Table 5 summarizes the regression analysis of reported lattice structure experimental data. According to Gibson–Ashby theory, the C value should be in the range of 0.1–4.0. In this study, the values of C are 0.17 and 0.39 for the SLM-ed and HT-ed alloys, respectively. The increase of C value after heat treatment also represents the increase of tensile strength and Young's modulus. On the other hand, the smaller the value of n, the higher the stability and utilization of the lattice material. Compared with other bending-dominant lattice structures, the stretch-dominant lattice structure has a smaller n value and higher load-bearing capacity. Correlations were generally higher for strength than they were for modulus, though only marginally. Comparisons between the derived exponents and those predicted by the Gibson–Ashby model for bending and stretch-dominated structures provides insight into the general behavior of these topologies. For example, the exponent for the strength of FCCZ (n = 1.59) is very close to that predicted for stretch-dominated structures (n = 1.5), suggesting this topology behaves more like a stretch-dominated structure despite being bending-dominated according to the Maxwell criterion. In addition, in Table 5, it is found that when the simulated n value is less than 1, its R2 value drops rapidly, and can even be as low as 0.094. However, the values of n in this work are all lower than 1, and the minimum R2 is 0.536 and the maximum is 0.975. This shows that the structure of this study is more stable and reliable.
Regression analysis of reported and present lattice structure experimental data
Unit cell topology
Tensile strength
Young's modulus
Ref.
n1
C1
Correlation (R2)
n2
C2
Correlation (R2)
F2CCZ-SLM-ed
0.69
0.17
0.847
0.35
0.19
0.536
This work
F2CCZ-HT-ed
0.79
0.39
0.975
0.48
0.37
0.893
This work
Diamond
0.55
0.08
0.094
0.86
0.04
0.220
39
FCC
1.75
1.18
0.998
1.94
0.37
0.999
40
BCCZ
0.48
0.06
0.388
0.22
0.01
0.124
41
Dodecahedron
2.51
2.21
0.990
2.64
0.63
0.999
41
Schwartz diamond
1.36
0.77
0.997
0.76
0.05
0.966
42
FCCZ
1.59
1.32
0.994
1.50
0.35
0.919
41
FBCCZ
1.72
1.40
0.953
0.54
0.03
0.466
40
Conclusions
In this work, Ti64 lattice structure was designed and fabricated by SLM, and the influence of double solid solution-aging treatment on Ti64 lattice structures, including the microstructure and tensile properties were investigated. The following conclusions were obtained.(1) The size of the F2CCZ struts along the build direction is close to the design value, and the error does not exceed 30 μm. However, the dimensional error of the inclined struts at 45° to the build direction can be up to 140 μm.(2) SLM-ed specimens have high dislocation density and residual stress, causing serious stress concentration at the junction of lattice and solid under tensile load. The double solution aging treatment not only eliminates the residual stress and dislocation, but also produces multi-scale microstructures (multi-scale α laths and β nanoparticles). The F-4 sample has the largest increase in tensile strength is from 140 ± 18 MPa for the SLM-ed alloy to 229 ± 5.1 MPa for the HT-ed alloy.(3) According to the Gibson–Ashby model, the n values of SLM-ed and HT-ed alloys are both less than 1, which means that the F2CCZ structure has high stability and material utilization. However, the correlation coefficient (R2) of the HT-ed alloy is closer to 1, indicating that the sample is more suitable for the Gibson–Ashby model.
Authors: Harun H Bayraktar; Elise F Morgan; Glen L Niebur; Grayson E Morris; Eric K Wong; Tony M Keaveny Journal: J Biomech Date: 2004-01 Impact factor: 2.712
Authors: S M Ahmadi; G Campoli; S Amin Yavari; B Sajadi; R Wauthle; J Schrooten; H Weinans; A A Zadpoor Journal: J Mech Behav Biomed Mater Date: 2014-02-08
Authors: Duyao Zhang; Dong Qiu; Mark A Gibson; Yufeng Zheng; Hamish L Fraser; David H StJohn; Mark A Easton Journal: Nature Date: 2019-12-04 Impact factor: 49.962
Authors: S Amin Yavari; R Wauthle; J van der Stok; A C Riemslag; M Janssen; M Mulier; J P Kruth; J Schrooten; H Weinans; A A Zadpoor Journal: Mater Sci Eng C Mater Biol Appl Date: 2013-08-13 Impact factor: 7.328
Authors: C J Todaro; M A Easton; D Qiu; D Zhang; M J Bermingham; E W Lui; M Brandt; D H StJohn; M Qian Journal: Nat Commun Date: 2020-01-09 Impact factor: 14.919
Authors: Roland Chapurlat; Minh Bui; Elisabeth Sornay-Rendu; Roger Zebaze; Pierre D Delmas; Danny Liew; Eric Lespessailles; Ego Seeman Journal: J Bone Miner Res Date: 2019-12-10 Impact factor: 6.390