Diffraction-based determination of single-crystal elastic constants of polycrystalline titanium alloys

The single-crystal elastic constants of polycrystalline titanium alloys have been determined using neutron and synchrotron powder diffraction.


Introduction
Single-crystal elastic constants are essential material parameters in fundamental materials science as well as in engineering. In particular, the elastic properties of any polycrystalline bulk material are based on its single-crystal elastic constants. A variety of micromechanical models have been established to describe the relations between the elastic properties of a polycrystalline bulk material and the singlecrystal elastic constants, the most common applied approaches being those introduced by Voigt (1928), Reuss (1929), Hill (1952), Kroener (1958), de Wit (1997) and Matthies et al. (2001).
In classical stress analysis using diffraction, the residual stresses of alloys are determined by measuring the lattice strains for different sample orientations. Lattice strains and stresses are related by the diffraction elastic constants (DECs), which depend on the direction in the crystal as described by the lattice planes (hkl) of the corresponding Bragg reflections. DECs in turn are based on the single-crystal elastic constants. Thus, knowledge of the single-crystal elastic constants is essential for diffraction-based stress analysis.
The conventional stress analysis can be modified in such a way as to apply a defined external stress on the sample (by tensile load or compression) while the lattice strains are measured under various sample orientations. The singlecrystal elastic constants of any polycrystalline material may then be determined by an inversion of the usual calculations applied in the residual stress analysis and fitting by a 2 -minimization technique. The feasibility of such an ISSN 1600-5767 approach was demonstrated by Hauk & Kockelmann (1979) and later applied by various authors (Gnä upel-Herold et al., 1998;Kisi & Howard, 1998;Singh et al., 1998;Howard & Kisi, 1999;Matthies et al., 2001;Fré our et al., 2005;Stebner et al., 2013;Lunt et al., 2014;Kim et al., 2016) for alloys and ceramics by diffraction under tensile stress or compression. This technique is of particular relevance for any engineering material which cannot be obtained as a single crystal, especially multiphase or highly twinned materials.
In this contribution, we will first review the background to this technique, with particular focus on the measurement geometry and data evaluation. For the validation of our data collection and evaluation procedure, face-centred cubic (f.c.c.) and body-centred cubic (b.c.c.) phases in steels and cast iron samples were studied. The results are compared with literature data on single-crystal elastic constants and bulk properties. Approaches for the modelling of grain-to-grain interactions by Voigt (1928), Reuss (1929), Hill (1952), Kroener (1958), de Wit (1997) and Matthies et al. (2001) are applied in the calculations. Since in earlier work the influence of the texture was not investigated in detail, we systematically examine here the influence of texture weightings compared with a quasiisotropic approach. In addition, a load-partition model is included in the modelling approaches.
Despite the fact that titanium alloys are important highperformance alloys, few data are available concerning DECs or single-crystal elastic constants. Howard & Kisi (1999) obtained the single-crystal elastic constants in the near -alloy Ti-6Al-4V by diffraction on a polycrystalline specimen. In a similar way, Fré our et al. (2005) determined the elastic constants of the b.c.c. phase in the dual-phase alloy Ti-17.
Ti64 is a high-strength titanium alloy and is considered as the 'workhorse' alloy for aerospace applications. It offers excellent strength and toughness up to 673 K combined with good fabricability. As a 'near -alloy' it consists mainly of the hexagonal phase and low amounts of a cubic phase.
Ti6246 is a very high strength titanium alloy but with lower toughness and weldability than Ti64, although it can be operated at higher temperatures. It consists of two phases, one hexagonal phase and one b.c.c. phase.
In general, -alloys have high corrosion resistance and very high strength but a lower elastic modulus than -alloys.

Theory
In the following, a general derivation of the current methods is provided. This includes the single-crystal elastic constants, their relation to the crystal lattice and their significance for the DECs. For multiphase alloys, a method of quantifying the load transfer from one phase to another as described by Behnken (2003) is implemented in the evaluation process. A fitting routine for the elastic constants based on the approach by Gnä upel-Herold et al. (1998) is applied to different micromechanical models.
The well known Hooke's law (1) provides the relation between the second-rank tensors of strains and stresses for any material under elastic strain. The proportional constants are given by the fourth-rank tensors of elastic compliances A and elastic constants C.
Diffraction studies under the influence of an applied mechanical load enable investigation of the strains perpendicular to particular (hkl) lattice planes, i.e. strains in different crystallographic directions. As outlined in detail by Gnä upel-Herold et al. (1998), one needs to consider the relations between three coordinate systems. The measurement frame is defined by the scattering vector Q coinciding with the reciprocal-lattice vector h for each lattice plane (hkl). The load axis determining the applied stress defines the sample frame, while the single-crystal elastic compliances are expressed in the crystal frame. As illustrated in Fig. 1, the measurement frame L is defined in such a way that the L 3 axis is parallel to Q. In the sample frame S 3 is oriented along the load axis in the case of uniaxial tensile (or compression) experiments. As shown in Fig. 1, the orientation between the scattering vector Q and the load axis (S 3 ) is given by the angles and . S can be transformed into L via the rotation ! as given in equation (2) The transformation from the crystal frame into the measurement frame is done with the rotation (Behnken, 2003): In order to ensure the correct orientations of L 3 to the crystal frame, the rotations are defined by the Euler angles f þ =2; ; 'g and f' 2 ; 2 ; =2 À 2 g (Gnä upel- Herold et al., 1998;Brakman, 1983). From equation (1) the lattice strains in the measurement frame for a Bragg reflection (hkl) under orientations ' and are given by with the stress factor F ij given by A 0 uwij ðgÞ are the general single-crystal compliances expressed in the measurement frame. Therefore A from equation (1) is transformed from the crystal into the measurement frame by the rotation matrix in the following way: The function f(g) describes the orientation distribution of the grains (ODF) and is typically related to the sample frame (Behnken, 2003), i.e. in equation (5) the integrals are evaluated over all grains contributing to the diffraction signal. In equation (4) the strains are measured along Q, and the only observed component of the strain tensor is 33 .
Equation (7) can be transformed with the knowledge of the rotation symmetry of the average tensor and the measurement direction in the measurement frame, yielding ð'; ; hklÞ ¼ hA 0 3311 i r tr 0 ð Þ þ hA 0 3333 i r À hA 0 Equation (9) can be transformed into the sample reference frame, which leads to the general equation of stress analysis (Behnken, 2003) 33 ¼ s 1 ð 11 þ 22 þ 33 Þ þ 1 2 s 2 cos 2 ð'Þ 11 þ sin 2 ð'Þ 22 with the DECs s 1 (hkl) and 1 2 s 2 ðhklÞ s 1 ðhklÞ ¼ hA 0 3311 i r ; ð11Þ 1 2 s 2 ðhklÞ ¼ hA 0 3333 i r À hA 0 3311 i r : An analytical solution of DECs from the single-crystal elastic constants can be calculated for different crystal symmetries using a model assumption for the grain-to-grain interaction.
In order to establish relationships between bulk and singlecrystalline properties, Voigt assumed that in all grains homogeneous strains appear as a function of an external stress. In this case the general elastic compliances are calculated from the averaged elastic constants hc(g)i À1 (Voigt, 1928). Later, Reuss (1929) assumed homogeneous stress in grains; thus A depends on the elastic compliances. Kroener (1958) developed a model based on Eshelby's assumption of spherical elastic inclusions in an isotropic matrix (Eshelby, 1957), which was extended by de Wit (1997). Thus, Reuss, sðgÞ þ tðgÞ Eshelby/Kroener.

<
: ð13Þ The DECs based on the Voigt (1928) model for any crystal symmetries are given in equations (14) and (15): ; ð14Þ x ¼ c 11 þ c 22 þ c 33 ; In equation (15) the single-crystal elastic constants are represented in Voigt's notation. The solutions for Reuss's approximation have been developed by Behnken (2003) for arbitrary crystal symmetries: 2 s 2 ðhklÞ ¼ 1 2 m n ð3 p q À pq Þ A mnpq ; ð17Þ On the basis of these equations, one can derive the DECs for different crystal symmetries. In the special case of cubic crystal symmetry the DECs can written as with A 0 = A 1111 À A 1122 À 2A 1212 and the crystal-orientation parameter À: It has been shown by Hill (1952) that the assumptions made by Voigt (1928) and Reuss (1929) are only valid in specific cases, giving an upper and lower limit for the bulk elastic properties for most materials. Therefore, the arithmetic average of the bulk moduli obtained by both models was proposed for more realistic values. However, due to the tensor nature of the elastic constants, the physical relation A = C À1 is not granted in this case. Matthies et al. (2001) showed that a geometric average between Voigt (1928) and Reuss (1929) calculations obeys the physical relation A = C À1 . Kroener (1958) introduced a self-consistent model to solve the case of spherical inclusions in a matrix for cubic symmetries where the shear modulus G of the matrix is compared with the shear modulus G 0 of the inclusion: with the parameterization and DECs de Wit (1997) modified Kroener's approach, recognizing that, by applying statistical averages, the stress is averaged over all spatial directions, whereas the strain is only measured in the direction normal to the diffraction planes. This leads to a slight modification in the description of the parameters (de Wit, 1997): In a multi-phase alloy subjected to elastic strain, a load transfer towards stiffer phases is expected. Diffraction studies reveal the effective elastic constants of the individual phases due to their interactions with the other phases present. Behnken (2003) derived an approach to quantify the load transfer in multiphase systems with different elastic properties. The mean phase stress consists of a macro stress-dependent contribution and an independent part denoted by the index 0.
f is the fourth-rank tensor of the stress transition factors for phase , L is the applied load to the sample and I are the residual stresses. To calculate the transition factors, the following equation is derived for an inclusion in a homogeneous matrix (Behnken, 2003): In the fourth-rank tensor equation (27), C are the singlecrystal elastic constants of the phase (in this case the inclusions), I is the unity tensor, and C and S are the macroscopic elastic moduli and compliances, respectively, based on the following equation: where the phase fractions of phases and are denoted by p , , Young's modulus is denoted by E and is a Lame constant. w À1 depends on the shape of the inclusions and is often called the Eshelby tensor. An analytical solution has been found by Kneer (1965) for spherical and ellipsoidal inclusions: In the case of spherical inclusions, w À1 is isotropic and depends only on the Poisson ratio (Kneer, 1965). The stresses of the matrix can be calculated either via the transition factors of the matrix which are obtained via p f + p f = I or via the condition ð1 À pÞh m i + p i = , where the index m represents the matrix and i the inclusion, and p is the phase fraction of the inclusion. One needs to address the minimization problem to obtain the single-crystal elastic constants from diffraction on polycrystalline materials. We choose the method of least squares ( 2 minimization) to optimize the elastic compliances by minimizing the sum of squares between deviations of measured and calculated observables.
The first possibility is to derive the DECs from the measured data and subsequently minimize the differences between them and the computed elastic constants (Gnä upel-Herold et al., 1998), In equation (30), S is a two-dimensional quantity and the total length must be minimized.
research papers The technique for obtaining the values of S 1 (hkl) meas and 1 2 S 2 ðhklÞ meas is explained in the following section. S 1 (hkl, A ijkl ) calc and 1 2 S 2 ðhkl; A ijkl Þ calc are computed with the help of equations (16), (18), (14) or (24), depending on the model used.
The second possibility is to express 2 directly in terms of the measured strains 33 , leading to the following expression for 2 : In principle, both techniques will lead to similar results, but there are important differences between them. Equation (31) is a more specialized form of equation (32), which leads to a large reduction in the summation terms of 2 , e.g. k ) n because the S i are combined values from the different measurement directions and only depend on the scattering planes.
Additionally, in both cases the texture of the sample can be taken into account through weightings of either the DECs or the strains directly, based on the multiples of random distribution (m.r.d.) values from the ODF of the different scattering planes involved.

Experimental setup
Diffraction studies were carried out using both neutron and high-energy X-ray diffractometers. The large gauge volume in neutron experiments yields higher grain statistics compared with X-ray studies. On the other hand, high-energy (shortwavelength) X-ray experiments enable the recording of complete Debye-Scherrer rings under a continuous increase of the applied load.
A load frame especially designed for elastic anisotropy and texture measurements was used, allowing the load axis to be oriented with respect to the incident beam by a rotation as on an Eulerian cradle (Hoelzel et al., 2013).
Neutron diffraction experiments were performed on the high-flux powder diffractometer D20 at the Institut Laue-Langevin in Grenoble, France (Hansen et al., 2008) and the high-resolution powder diffractometer SPODI at the Heinz Maier-Leibnitz Zentrum (MLZ) in Garching near Munich, Germany (Hoelzel et al., 2012). Both instruments operate at a constant wavelength and offer a broad scattering angle range up to 160 in 2Â. In this setup the position of the sample frame with respect to the measurement frame is given by (tilt axis of the load frame) and (rotation axis of the sample table), as shown in Fig. 2. This setup allows a simultaneous analysis of various (hkl) reflections, while the individual orientations between the scattering vectors Q(hkl) and the load axis need to be considered in the data analysis for each reflection individually. By changing the orientation of the load axis by and , different angles and ' can be obtained for each (hkl) reflection. In addition, measurements at different angles can be used to account for possible texture effects.
The experiments were done using standard tensile samples of 6 or 8 mm in diameter under uniaxial tension; no shear or bending forces were applied, nor did they occur on the sample during the tensile experiments. The macroscopic strain was measured along the load axis using a strain gauge (clip-on extensometer Sandner EXA-15). For each sample composition, tensile tests were carried out prior to the neutron diffraction studies to determine Young's modulus and the yield stress. Diffraction patterns were collected at a minimum of three and up to five stress levels, up to a maximum value of about 70% of the yield stress. At each load step, measurements at a minimum of five and up to ten angles were performed.
To validate our method for data collection and analysis, ferritic (b.c.c.) and austentic (f.c.c) phases in steels and cast iron samples were studied on the SPODI diffractometer at a monochromator angle of 155 using Ge(551) for a wavelength of 1.5483 Å . The instrument offers a good angular resolution to separate the ferrite and austenite peaks and to observe their shifts under load at high precision. The investigations were carried out on ferritic structural steel S235JR, an austenitic stainless steel of AISI type 304 (X5CrNi 18-10), duplex steel (X2CrNiMoN 22-5-3) and an austempered ductile iron (ADI) sample consisting of ferrite, austenite and graphite.
The neutron scattering power of titanium alloys is quite poor. This is particularly the case for the phase of titanium owing to the negative scattering length of Ti and a significant quantity of alloying elements (with positive scattering lengths). Therefore, the investigations on titanium alloys  A schematic view of the rotatable tensile rig. The orientation between the sample and the incident beam is given by the angles , and '.
X-ray diffraction experiments were carried out on Ti-6Al-2Sn-4Zr-6Mo on the High Energy Materials Science beamline (HEMS) at the synchrotron facility DESY in Hamburg, Germany (Schell et al., 2014). An energy of 98.25 keV (corresponding to a wavelength of 0.12619 Å ) was used to investigate tensile samples of 6 mm diameter. Since full Debye-Scherrer rings were obtained, no rotation was necessary to cover the different orientations of the load axis, which was kept perpendicular to the incident beam, i.e. = 90 . The data were collected in 10 N steps up to 30 kN.
For texture analysis, neutron pole-figure measurements were carried out on the STRESS-SPEC instrument at the MLZ (Brokmeier et al., 2011) using a Ge(311) monochromator for a wavelength of 1.68 Å . Bragg reflection intensities were measured from 0 to 360 in 5 steps in ' in seven different orientations in .
Preparation for electron backscatter diffraction (EBSD) and energy-dispersive X-ray diffraction (EDX) was performed by mechanical cutting and polishing. After cutting, plane grinding was performed with an MD-Mezzo surface followed by a single fine grinding with a 9 mm diamond suspension abrasive on an MD-Largo surface. The first polishing was performed with an MD-Chem surface and a 0.04 mm collodial silica abrasive and a final step of ion polishing with a Hitachi IM4000PLUS cross section ion polisher operated at 5 kV and 200 mA for 60 min with sample oscillation. The sample was cut along the cylinder axis, such that the sample normal is the transverse direction of the cylinder. EBSD maps were measured on a Hitachi SU5000 field-emission scanning electron microscope equipped with an Oxford Instruments NordlysII EBSD detector and an ULTIM MAX EDX detector. The employed acceleration voltage was 20 kV. EBDS and EDX signals were collected simultaneously with the sample inclined by 70 towards the EBSD detector. Data evaluation was performed with the Oxford Instruments AZTec and CHANNEL5 software. The EDX analysis was carried out in standardless mode.

Data evaluation
For the data evaluation, a software package was developed enabling the choice of all mentioned models and, optionally, the texture implementation. For dual-phase alloys, a possible load transfer from one phase to the other can be quantified.
The hkl-dependent strains ðhklÞ = ½dðhklÞ À d 0 =d 0 were obtained from the changes in the reflection positions under load in a set of diffraction patterns. An example diffraction pattern for Ti64 measured on D20 is shown in Fig. 3. To fit all free variables, at least three or five different Bragg peaks are needed for the cubic and hexagonal symmetries, respectively. The Bragg reflections were chosen to cover the maximum range of crystal directions to track the dependencies of the compliance tensor components.
For the ferritic phase (Im3m)  A diffraction pattern for Ti64 measured on D20 at = 45 and = 90 . The artifact seen near 2Â = 80 is caused by the detector and does not influence the results, since no hkl reflection is affected.

Figure 4
The linear relationship between the strain divided by the applied stress and cos 2 ( ) is used to fit the DEC of the (220) plane of Ti38644. vector Q. An example fit of the (220) DEC of Ti38644 is given in Fig. 4. It should be emphasized that this evaluation is not affected by intergranular residual stresses, as long as they do not change during the elastic loading, because the lattice strains are given by shifts of the reflections under load rather than by absolute d values.
The The load transfer is calculated by a self-consistent scheme until the change in the transition factors after an iteration is smaller than a given limit. The initial values for the elastic constants are used in this approach, with no consideration of the load transfer. The stress factors for the inclusion are calculated with equation (27) defining the stress factors for the matrix. The parameters 0 and I in equation (26) are zero in our case, as the effects of residual stresses cancel out during the evaluation procedure as mentioned earlier: By means of equation (34) the associated stress for each Bragg reflection can be updated and the adjusted elastic constants are calculated. This iteration is repeated until the change in the stress factors is smaller than 10 À5 .

Results
All obtained data were fitted with all available models except the Voigt (1928) model. The calculations based on the approach by Voigt showed large instabilities concerning the fitting routine, led to multiple results or did not converge. Therefore the results for Voigt (1928) are only shown once, for the structural steel S235JR. The values for the single-phase alloys agree well with available literature data. The results for the ferritic phase are listed in Table 1. Finkel (2016) and Gnä upel-Herold et al.  Table 2 illustrates the results for the austenitic phase calculated for different models compared with the literature data. The data published by Ledbetter (1985) Table 1 The single-crystal elastic constants and bulk properties for S235JR measured on SPODI.
The anisotropy A is calculated after Zener (1948). For the Voigt model, the fitting routine does not converge to a single solution.
Model  204.8, 138.9 and 124.2 GPa for the texture adaptation shows good agreement in both cases.
The results for the duplex steel and ADI are shown in Tables 3 and 4, respectively. Both consist of austenitic and ferritic phases. However, ADI contains a large amount of graphite of approximately 10 vol.% in the form of nodules with Young's and shear moduli of essentially zero. For c 11 and c 12 the deviation is of the same order of magnitude as the uncertainty in the values, but for c 44 the observed deviation is higher. In tensile experiments the parameter c 44 is related to the shear stresses/strains and is therefore only indirectly accessible with the different orientations covered during the measurements. This leads to higher uncertainties for c 44 during the determination in the fitting process. The accuracy for c 44 may be improved by including torsion experiments, as suggested by Woracek et al. (2012).
The tables also reveal larger discrepancies between the models in the Zener anisotropy defined by A = 2c 44 = ðc 11 À c12Þ (Zener, 1948) and the c 12 =c 11 ratio. In a systematic study of the single-crystal elastic constants of different monocrystals of austenitic stainless steels and Fe-Cr-Ni alloys, Ledbetter (1985) found that both ratios remain nearly  Table 2 The results for AISI type 304 (X5CrNi 18-10), a single-phase stainless steel, measured on SPODI; the texture data were measured on STRESS-SPEC.
The anisotropy A is calculated after Zener (1948).
Model  Table 3 The single-crystal elastic constants for X2CrNiMoN 22-5-3 duplex steel, the dual-phase stainless steel alloy, measured on SPODI; the texture data were measured on STRESS-SPEC.
The anisotropy A is calculated after Zener (1948).
Model constant at A = 3.53 and c 12 =c 11 = 0.635 (Ledbetter, 1985). Owing to the low variance in the values of A and c 12 =c 11 in Ledbetter's study, we believe these ratios can be used to estimate the best model for the material under investigation. In all models it turns out that the texture does not influence the results above their uncertainties. This is most likely due to the elastic measurements being performed on single-crystal domains and therefore not directly affected by the texture. The only parameter affected is the average strain measured via diffraction. Matthies et al. (2001) concluded that just a couple of thousand grains are enough to ensure the statistical significance of the average strains. Therefore, the isotropic  Table 4 The single-crystal elastic constants for austempered ductile iron consisting of ferrite, austenite and graphite measured on SPODI; the texture data were measured on STRESS-SPEC.
The anisotropy A is calculated after Zener (1948

Figure 6
An overview of the single-crystal elastic constants of the h.c.p. phase in Ti64 and Ti6246.
grain orientation assumption is considered to be adequate for the investigated materials as the texture is assumed not to change significantly in the elastic regime with the applied load. Further, our diffraction data show no changes in the Bragg intensities in the elastic regime. The elastic properties of both phases in duplex steel are very similar, resulting in a load transfer of about 0.3%. Thus, the load transfer approach reveals practically the same elastic constants. A similar behaviour would be expected for ADI but was not investigated further here as the additional graphite phase could not be taken into account.
The results for the obtained elastic constants in the hexagonal phase in Ti64 and Ti6246 are shown in Fig. 6, while An overview of the single-crystal elastic constants of the b.c.c. phase in Ti38644 and Ti6246. Table 5 The single-crystal elastic constants and macroscopic values for Ti38644 measured on D20.
The anisotropy A is calculated after Zener (1948 Table 6 The single-crystal elastic constants for Ti64 measured on D20; texture data were measured on STRESS-SPEC. Model a dual-phase ( + ) titanium alloy. Tables 5, 6 and 7 reveal quite similar results for the different model assumptions. In addition, the values obtained for the + alloy Ti6246 agree quite well, especially if the load transfer is taken into account, with the corresponding data for the phase in Ti64 and the phase in Ti38644, respectively. Table 6 also includes the results of Howard & Kisi (1999) on the Ti64 alloy determined by the Reuss (1929) approach, as well as values obtained by ultrasonic studies on single crystals of pure Ti (Fisher & Renken, 1964).
Good agreement with ultrasonic data was found in the hexagonal phases for c 11 and c 44 . The largest deviation between our results and earlier work is found for c 33 , and minor deviations in the c 12 and c 13 elastic constants. In hexagonal systems the behaviour in the ab plane is completely isotropic owing to the requirement for the c 66 = ðc 11 À c 12 Þ=2 elastic constant. As we are able to determine c 11 and c 12 accurately, this also gives access to the shearing parameter c 66 . On the other hand, hexagonal systems show additional anisotropies regarding the c axis which influence the precision of parameters c 12 and c 13 .
Macroscopic Young's moduli were measured in a series of tensile tests to 106.6 (2.1) GPa for Ti64, to 88.5 (1.6) GPa for Ti38644 and to 113.6 (2.6) GPa for Ti6246. Tables 5, 6 and 7 indicate that there is very good agreement for Ti64 and Ti38644 between these and the values calculated on the basis of the diffraction studies, while somewhat higher values were obtained for Ti6246.
The phase fractions of the dual-phase Ti6246 alloys were evaluated from the obtained diffraction patterns using Rietveld analysis (MAUD and FullProf suite) (Lutterotti et al., 1997;Rodriguez-Carvajal, 1993). For the alloy, phase fractions of 78 and 22% were found for the and phases, respectively. The results of phase distribution versus phase composition obtained by EBSD and EDX are shown for the Ti6246 alloy in Fig. 8. The phase fractions determined by EBSD amount to 78 and 22%, the same as for the Rietveld analysis. The main chemical difference between the two phases resides in the Mo content and the compositions are tabulated in Table 8. While the doping elements Sn, Zr and Al are relatively homogeneously distributed, Mo shows an inhomogeneous distribution with elliptical Mo-poor regions in the 5-10 mm range (Fig. 8).
The load partitioning applied to the Ti6246 sample indicates a significant load transfer from the phase to the phase. The stress in the phase is reduced by 11.01% and transferred to the phase, increasing its stress by 3.10% (note that the smaller increase is mainly due to the ratio of phase fractions of around 1:4). The single-crystal elastic constants corrected for load transfer are shown in Table 7. They reveal a clear shift from the uncorrected values towards the corresponding values of the single-phase samples. The impact of the load transfer in the phase for the Reuss (1929), Hill (1952) and Matthies et al. (2001) models yields a change of about 10% in most cases, whereas for the Kroener (1958) and de Wit (1997)  The single-phase values of 124.9, 69.9 and 38.0 GPa for the same model almost match the 'load-transfer-corrected' ones. Similar behaviour is also observed for all other models. Owing to the lower increase in stress in the phase the changes remain rather small, but the results achieved with the loadtransfer model match the single-phase values more consistently. The comparatively small changes in the elastic constants when using the Kroener (1958) and de Wit (1997) models are supposedly due to the assumption of inclusions in a matrix, which is essentially the same assumption as used for the load transfer.
Similar elastic constants for phases in different titanium alloys can be expected owing to the low quantity of alloying elements, shown already by the comparison of Howard & Kisi (1999) with pure titanium. Significant variations in the elastic constants for phases in different alloys can be found in the literature (Hounkpati et al., 2016;Fré our et al., 2005). This may result from the fact that -alloys contain a large number of different -stabilizing elements in contrast with the phase. However, our investigations show a larger impact of the modelling between the two investigated phases.

Conclusions
Diffraction experiments are a feasible and powerful route to determine single-crystal elastic constants from polycrystalline multiphase alloys due to the phase selectivity, as shown in the case of Ti6246.
The data-evaluation procedure included several approximations to consider grain-to-grain interactions in polycrystals based on the approches by Voigt (1928), Reuss (1929), Hill (1952), Kroener (1958), de Wit (1997) and Matthies et al. (2001). In terms of computational stability, the Reuss (1929), Hill (1952) and Matthies et al. (2001) assumptions were the most stable in our data analysis. Concerning the approaches of Kroener (1958) and de Wit (1997), not all possible combinations of the elastic constants led to a solution of the selfconsistent equation (22). The Voigt equations could only be used for the evaluation with very limiting constraints. Good indicators of which models will suit best are the anisotropy values and the c 12 /c 11 ratios obtained for the different models. We systematically investigated the effects of the texture on the modelling and find its influence to be smaller than the differences observed for different grain-to-grain interaction models.
For the first time a full analysis of the load transfer has been included in the evaluation of the elastic constants from diffraction data for the example of Ti6246.
For the dual-phase alloy Ti6246, the load-transfer approach allows a direct comparison of measured 'effective' elastic constants with 'load-transfer-corrected' elastic constants. In this alloy, a significant load relocation of about 14% was found from the to the phase. By including the load transfer, the evaluations of the constants of the phase in the dual-phase alloy were significantly shifted towards those of the singlephase alloys, while the effect of the load transfer on the phase was significantly smaller. In particular, the loadtransfer-corrected elastic constants in the phase of Ti6246 show good agreement with the corresponding values for the near -alloy in Ti64. In addition, for the phase the loadtransfer-corrected constants are in excellent agreement with the corresponding single-phase result for the pure -alloy.