Hostname: page-component-745bb68f8f-grxwn Total loading time: 0 Render date: 2025-02-11T10:42:25.712Z Has data issue: false hasContentIssue false
  • English
  • Français

A software for diffraction stress factor calculations for textured materials

Published online by Cambridge University Press:  15 June 2012

Thomas Gnäupel-Herold
Thomas Gnäupel-Herold*
Affiliation:
Materials Science and Engineering, University of Maryland, Building 090, Room 2135, College Park, Maryland 20742 Center for Neutron Research, NIST, 100 Bureau Dr stop 6102, Gaithersburg, Maryland 20899-6102
*
a)Author to whom correspondence should be addressed. Electronic mail: tg-h@nist.gov
Rights & Permissions [Opens in a new window]

Abstract

A software for the calculation of diffraction elastic constants (DEC) for materials both with and without preferred orientation was developed. All grain-interaction models that can use the crystallite orientation distribution function (ODF) are incorporated, including Kröner, Hill, inverse Kröner, and Reuss. The functions of the software include: reading the ODF in common textual formats, pole figure calculation, calculation of DEC for different (hkl,φ,ψ), calculation of anisotropic bulk constants from the ODF, calculation of macro-stress from lattice strain and vice versa, as well as mixture ratios of (hkl) of overlapped reflections in textured materials.

Keywords

diffractionelastic constantstexturestrainstress
Type
Technical Articles
Information
Powder Diffraction , Volume 27 , Issue 2 , June 2012 , pp. 114 - 116
Copyright
Copyright © International Centre for Diffraction Data 2012

I. INTRODUCTION

Diffraction-based stress analysis depends critically on the use of the correct diffraction elastic constants (DEC). A survey through the literature makes it abundantly clear that in the vast majority of cases in which lattice strain needed conversion to stress-preferred grain orientation (texture) is – justifiably or not – disregarded, and isotropic diffraction elastic constants were used. The main reasons for this apparent oversimplification are the added need to quantify the degree of preferred grain orientation through pole figure measurements, and the use of pole figures to determine the orientation distribution function (ODF). Even with the ODF, further calculations are hampered by the lack of freely available software tools to perform calculations of anisotropic diffraction elastic constants. The IsoDEC (Gnäupel-Herold, Reference Gnäupel-Herold2011) software program introduced here was developed to address this need.

II. FEATURES AND FUNCTIONALITY

The main purpose of IsoDEC is the calculation of the orientation and (hkl) dependent stress factors using the ODF (in the textual form as provided by the output of freely available texture packages), the single crystal elastic constants and the bulk elastic constants as input. The latter can be calculated from the ODF. The software can also be used to calculate stress from measured lattice strains using the thus determined stress factors or, in reverse, calculate lattice strains from a given stress tensor. Some of the implemented models for DEC calculations also allow including the grain shape, expressed as a general ellipsoid, into the calculation of the DEC.

Evidently the ODF is crucial. It needs to be supplied in textual form, i.e., using the output from popLA (Kallend et al., Reference Kallend, Kocks, Rollett and Wenk1991) or MTEX (Hielscher and Schaeben, Reference Hielscher and Schaeben2008). The ODF is expected in a 5° × 5° × 5° equidistant spacing for the three Euler angles φ 1 Φ φ 2 (Bunge notation). This tabular form is internally expanded to the full interval [0…2π; 0…π; 0…2π] and made quasi-continuous by means of tri-linear interpolation. This allows the computation of pole figures (in popLA format) and orientation fibers (i.e., ODF values for a series of Euler angle combinations) through Euler space. Orientation fibers are a key element of DEC calculations that use the ODF. The models for DEC calculations are the Kröner model (Kröner, Reference Kröner1958; Bollenrath et al., Reference Bollenrath, Hauk and Müller1967; Behnken and Hauk, Reference Behnken and Hauk1986), the Reuss model (Reuss, Reference Reuss1929; Möller & Martin, Reference Möller and Martin1939), the Hill model (Voigt, Reference Voigt1928; Reuss, Reference Reuss1929; Hill, Reference Hill1952), the inverse Kröner model (Kröner, Reference Kröner1958; Gnäupel-Herold, Reference Gnäupel-Herold2012), and the modified Voigt model (Murray and Noyan, Reference Murray and Noyan1999). The two Kröner-type models allow the user to devise complex elastic interactions: first, one can approximate multi-phase composites (e.g. elastically hard particles in a soft matrix such as SiC in Al) by using appropriate set of single crystal constants and bulk constants. Second, the effects of elongated or otherwise non-spherical grains can be included by changing the grain-shape parameters in IsoDEC.

Input and output in IsoDEC can be done through the spreadsheet that allows the import and export of text files with its content as well as basic copy and paste exchanges with common spreadsheet programs such as EXCEL. Large datasets consisting of d-spacings, their uncertainties, unstressed d-spacings with uncertainties, and measurement directions (φ ψ) can be treated this way. Stress tensor components and the unstressed d-spacing can be fixed or free parameters in the stress fit, depending on boundary/equilibrium conditions.

III. OPERABILITY

IsoDEC offers two visual interfaces: one for isotropic calculations and one for calculations that include preferred orientation. The interface for isotropic calculations is shown in Figure 1. The fewest possible steps for calculating the isotropic constants s 1(hkl) and ½s 2(hkl) are loading a set of material constants from the drop-down list, entering the indices h, k, and l, and starting the calculation. Elastic constants that are not found in the database included can be entered manually, but must include the lattice parameters and, after entering the single crystal elastic constants, a calculation of the isotropic bulk constants (upper left).

Figure 1. Screenshot of IsoDEC for isotropic calculations.

For anisotropic calculations, the amount of information that has to be provided increases as shown in Figure 2. Generally the ODF has to be provided either in popLA format or in MTEX format (“Load ODF”). There is also the possibility of not using an ODF, and entering user values for the bulk constants, e.g. using a rotated single crystal constants tensor to simulate the effects of individual texture components on the stress factors. Of course, the anisotropic bulk constants can be calculated from the ODF, the average grain shape (as ellipsoid axis ratios) can be entered, and the model (Hill, Reuss, Kröner, inverse Kröner and Voigt/modified Voigt) can be chosen. The stress factors are calculated for the fields in which h, k, l, phi and psi are given. This step has to precede any stress or strain calculations that require the stress factors. A stress calculation needs the additional input of d-spacings, the d-uncertainties, the reference or unstressed d-spacing and their uncertainties, and a sorting qualifier (column 1 in the worksheet) that allows to group those data together that are to be used for the calculation of the same stress tensor (i.e., same location on the sample or similar). The reverse calculation d-spacing/strain from stress naturally requires stress values entered in their appropriate columns. The output of this calculation is in the column “d_calc”, which means that if strain is desired one has to perform the explicit calculation (d − d 0)/d 0 oneself.

Figure 2. Interface for anisotropic calculations.

IV. PLATFORM AND AVAILABILITY

IsoDEC was written in Turbo Delphi for the Windows operating system. It can also be used under MacOS X+ and x86 linux if Windows emulator software is present. Both Darwine (MacOS X and Wine) have been tested successfully. The installation package contains a help file, a database of single crystal elastic constants and several ODFs in textual form that are readable by IsoDEC for use in further calculations of DEC, pole figures, anisotropic bulk constants, and ODF intensities for given specimen directions. The program home page and the download link can be found at http://www.ncnr.nist.gov/programs/crystallography/software/isodec/.

DISCLAIMER

Certain commercial firms and trade names are identified in this report in order to specify aspects of the experimental procedure adequately. Such identification is not intended to imply recommendation or endorsement by the National Institute of Standards and Technology, nor is it intended to imply that the materials or equipment identified are necessarily the best available for the purpose.

References

Behnken, H. and Hauk, V. (1986). “Berechnung der röntgenographischen Elastizitäts-konstanten (REK) des Vielkristalls aus Einkristalldaten für beliebige Kristallsymmetrie,” Z. Met.kd 77, 620–626.Google Scholar
Bollenrath, F., Hauk, V. and Müller, E. H. (1967). “Zur Berechnung der vielkristallinen Elastizitätskonstanten aus den Werten des Einkristalls,” Z. Met.kd 58, 76–82.Google Scholar
Gnäupel-Herold, T. (2011). “IsoDEC: a software for calculating diffraction elastic constants,” J. Appl. Cryst. (Submitted).Google Scholar
Gnäupel-Herold, T. (2012). “A model for calculating diffraction elastic constants,” J. Appl. Cryst. 45(2), 197–206.CrossRefGoogle Scholar
Hielscher, R. and Schaeben, H. (2008). “A novel pole figure inversion method: specification of the MTEX algorithm,” J. Appl. Cryst. 41, 1024–1037.Google Scholar
Hill, R. (1952). “The elastic behaviour of a crystalline aggregate,” Proc. Phys. Soc. London 65, 349–354.CrossRefGoogle Scholar
Kallend, J. S., Kocks, U. F., Rollett, A. D., and Wenk, H.-R. (1991). “Operational texture analysis,” Mater. Sci. Eng. A 132, 1–11.CrossRefGoogle Scholar
Kröner, E. (1958). “Berechnung der elastischen Konstanten des Vielkristalls aus den Konstanten des Einkristalls,” Z. Phys. 151, 504–518.Google Scholar
Möller, H. and Martin, G. (1939). “Elastische Anisotropie und röntgenographische Spannungsmessung,” Mitt. Kaiser Wilhelm Inst. Eisenforsch. 21, 261–269.Google Scholar
Murray, C.E. and Noyan, I.C. (1999). “A modified Voigt method for calculation of the elastic constants of ensembles selected by diffraction methods,” Phil. Mag. A 79, 371–389.CrossRefGoogle Scholar
Reuss, A. (1929). “Berechnung der Fließgrenze von Mischkristallen aufgrund der Plastizitätsbedingung für Einkristalle,” Z. Angew. Math. Mech. 9, 49–59.CrossRefGoogle Scholar
Voigt, W. (1928). Lehrbuch der Kristallphysik (Teubner, Leipzig).Google Scholar
Figure 0

Figure 1. Screenshot of IsoDEC for isotropic calculations.

Figure 1

Figure 2. Interface for anisotropic calculations.