Next Article in Journal
Heteroepitaxy of Cerium Oxide Thin Films on Cu(111)
Previous Article in Journal
Aspherical Lens Design Using Genetic Algorithm for Reducing Aberrations in Multifocal Artificial Intraocular Lens
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Delineation of First-Order Elastic Property Closures for Hexagonal Metals Using Fast Fourier Transforms

Department of Mechanical Engineering, University of New Hampshire, Durham, NH 03824, USA
*
Author to whom correspondence should be addressed.
Materials 2015, 8(9), 6326-6345; https://0-doi-org.brum.beds.ac.uk/10.3390/ma8095303
Submission received: 25 July 2015 / Revised: 2 September 2015 / Accepted: 14 September 2015 / Published: 17 September 2015

Abstract

:
Property closures are envelopes representing the complete set of theoretically feasible macroscopic property combinations for a given material system. In this paper, we present a computational procedure based on fast Fourier transforms (FFTs) for delineation of elastic property closures for hexagonal close packed (HCP) metals. The procedure consists of building a database of non-zero Fourier transforms for each component of the elastic stiffness tensor, calculating the Fourier transforms of orientation distribution functions (ODFs), and calculating the ODF-to-elastic property bounds in the Fourier space. In earlier studies, HCP closures were computed using the generalized spherical harmonics (GSH) representation and an assumption of orthotropic sample symmetry; here, the FFT approach allowed us to successfully calculate the closures for a range of HCP metals without invoking any sample symmetry assumption. The methodology presented here facilitates for the first time computation of property closures involving normal-shear coupling stiffness coefficients. We found that the representation of these property linkages using FFTs need more terms compared to GSH representations. However, the use of FFT representations reduces the computational time involved in producing the property closures due to the use of fast FFT algorithms. Moreover, FFT algorithms are readily available as opposed to GSH codes.

Graphical Abstract

1. Introduction

Property closures are defined as the complete set of the single-crystal and polycrystalline homogenized anisotropic property combinations for a given material. These constructs are important for microstructure based material design. Within the applied mathematics community, the property closures are referred to as the G-closures. In the past, the G-closures have been calculated for properties such as effective conductivity and elastic stiffness for a limited set of two-dimensional microstructures comprised of isotropic phases [1,2,3,4]. More recently, a mathematical framework called microstructure sensitive design (MSD) has been conceived to treat microstructure as a continuous design variable in engineering design and optimization. MSD features invertible or bi-directional linkages between the statistical description of a microstructure and its effective properties [5,6]. These linkages facilitate the identification of microstructures satisfying targeted properties and component performance criteria through the use of the spectral representation of microstructure and material properties. This representation of microstructure facilitates building a space of all theoretically feasible microstructures. When statistical orientation distribution functions (ODFs) are the descriptor of microstructure, the complete set of theoretically feasible ODFs is called a texture hull [5,6,7,8,9] and can be quantitatively described by Fourier coefficients in the multidimensional Fourier space. Mapping the hull into a space of material properties using appropriate ODF-to-property relationships defines the property closures.
In prior work, the spectral representation necessary for MSD was implemented most commonly with generalized spherical harmonics (GSH) basis functions [10]. The GSH bases offer the most compact representations for describing the dependence of anisotropic material properties on crystal orientation because they can be customized to reflect the crystal symmetry as well as the sample symmetry [10,11]. These bases were used to delineate a number of property closures for both cubic and hexagonal metals including combinations of the elastic [7,12,13,14] and plastic [15,16] properties as well as calculations of functional properties [17]. The elastic-plastic property closures presented in these works were based on the first-order bounding theories [18,19,20]. To calculate the bounds of properties, material databases of microstructure invariant Fourier coefficients have to be built. With the databases, calculation of property bounds reduces to multiplication of the spectral coefficients for properties and those for texture (a point in the texture hull) and their summation. These calculations are essentially instantaneous, which is highly desirable for the fast computation of property closures and microstructure based material design. After calculating bounds for selected property combinations of interest by mapping the hull into property space, the problem at hand is to find textures corresponding to the boundary points of the property closure. To this end, several computational methodologies were formulated [5,12,21,22]. In the core of these methodologies is an appropriately formulated optimization procedure aimed at finding textures corresponding to the boundary points of the property closure. Interestingly, for all cubic metals, textures at the boundary points of the property closures correspond to the same set of textures [13].
Several examples exist in the current literature on using property closures to theoretically design and optimize microstructure to improve performances of components. These examples identified microstructures and associated properties maximizing the deflection in a compliant beam [5], the energy storage of a flywheel [23], and the in-plane load carrying capacity of a thin plate with a central circular hole [6] and minimizing the elastic driving force for crack extension in rotating disks [24] and internally pressurized thin-walled vessels [25]. The main microstructural feature governing the elastic and plastic properties in these studies has been the ODF, although in MSD, some limited consideration has been given to compositional variations in microstructures [26] and to fiber reinforced composites [9,27]. A design strategy termed the topology optimization has been successful in designing composites of extreme functional and mechanical properties [28,29,30,31].
In crystal plasticity [32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48], calculations of average Taylor factors and Mandel spins for all possible strain modes have also been accomplished using the GSH bases [49]. In addition, the GSH bases were used to build databases of precomputed crystal plasticity solutions for the stress, strain hardening and texture evolution function for cubic metals [50,51]. Evaluation of these functions reduced to multiplying the orientation invariant spectral coefficient with the orientation dependent GSH basis and summing as opposed to iteratively solving them using Newton’s solvers was found to accelerate crystal plasticity models. Recently, in place of the GSH representations, fast Fourier Transforms (FFTs) have been used in the ODF representation [52] as well as the crystal plasticity framework [16,53,54,55,56,57,58]. The main motivation for using FFTs is the availability of much more computationally efficient algorithms than those using GSH coefficients. The methods based on GSH are significantly less effective when applied to lower symmetry metals because the relevant number of dimensions in the Fourier space is considerably higher than that for cubic metals [7]. Additionally, FFT libraries are readily accessible while GSH codes are not as easily found. Delineation of property closures for cubic metals using FFTs has recently been accomplished [52].
In this paper, we have successfully formulated a new approach using FFTs to compute elastic property closures for HCP polycrystals. To this end, we have built a database of the orientation invariant Fourier transforms for components of the elastic stiffness and developed an optimization scheme for finding the boundary points of the convex property envelopes. This new approach is introduced in this paper and demonstrated with examples. The selected materials for this study were titanium, zinc, zirconium, beryllium, magnesium, and cobalt. These materials were chosen to validate the FFT based closures against reported property closures based on the GSH representation for orthotropic HCP metals [7]. It was shown that the closures obtained using the FFT representation are identical to those based on the GSH representation for the orthotropic terms of the stiffness tensor. The FFT representation allows delineation of first-order closures of elastic properties involving the normal–shear coupling stiffness/compliance tensor components. Although polycrystalline elastic properties were calculated from a given texture in the most general formalism of triclinic sample symmetry using the Voigt-Reuss-Hill approximation [11], this is the first report of this type of closures for HCP metals. Finally, we have shown that it is possible to build texture hulls using FFTs for HCP polycrystals.

2. Representation of ODF Using FFTs and Texture Hulls for HCP Metals

The Orientation Distribution Function (ODF), f ( g ) , is the probability density associated with the occurrence of the crystallographic orientation, g, in the sample of a polycrystalline material. An ODF can be mathematically expressed in its continuous or discrete form. While the former is required for the GSH representation, the latter is necessary for the FFT representation. In its discrete form, an ODF is expressed as [59]:
f b Δ g = N g ± Δ g / 2 N ,   b F Z f b Δ g = 1
where N is the total number of orientations in the sample; N g   ±   Δ g / 2 is the number of orientations that lie within a bin Δ g centered about g. FZ and b will be defined shortly. The orientation, g, can be described using an ordered set of three Bunge-Euler rotation angles ( ϕ 1 , Φ , ϕ 2 ) that bring into coincidence the crystal axis with the sample reference frame. The Bunge-Euler representation was employed because of the inherent periodicity of the space when defined as ( ϕ 1 [ 0 , 2 π ) ,  Φ [ 0 , 2 π ) , ϕ 2 [ 0 , 2 π ) ) . On this domain, any function of crystal orientations is guaranteed to be periodic and circumvent the Gibbs phenomenon [60], eliminating the need for additional high frequency terms in the FFT representations. The Bunge-Euler space is defined by the rotation angles that make their appearance in structure-property relationships of interest in the form of integer powers of sines and cosines of the rotation angles. The main disadvantages of the Bunge-Euler space are the numerous redundancies in the representation of crystal orientations (i.e., several locations in this space may correspond to a single crystal orientation) and the inherent distortion caused by the fact that the invariant measure in this space is proportional to Δ g = s i n ΦΔ ϕ 1 ΔΦΔ ϕ 2 [10]. In spite of the disadvantages of the Bunge-Euler space, the development of the spectral representations has been most successful for this space, largely because of the implicit periodicity in the material property functions of interest when defined in this space.
Knowing the function values on a uniform grid, FFTs can be computed [61,62,63,64,65]. In order to represent the ODF using FFTs, we sample the Bunge-Euler space using this grid, which defines the bins. The FFT representation of an ODF is expressed as:
f b 1 , b 2 , b 3 = 1 B 1 B 2 B 3 k 1 = 0 B 1 1 k 2 = 0 B 2 1 k 3 = 0 B 3 1 F k 1 k 2 k 3 e 2 π i b 1 k 1 B 1 e 2 π i b 2 k 2 B 2 e 2 π i b 3 k 3 B 3
where f b 1 , b 2 , b 3 denotes the value of the ODF at the grid point identified by b, and b 1 , b 2 , and b 3 enumerate points within uniformly discretized domain; F k 1 k 2 k 3 are the FFTs of the ODF indexed by k, while k 1 , k 2 , and k 3 are the indices of the FFTs; and B 1 , B 2 , and B 3 are the total number of grid points in the periodic domain. Since the values of the ODF are the real half of the transforms, the other half are complex conjugates i.e., F k 1 , k 2 , k 3 = F B 1 k 1 , B 2 k 2 , B 3 k 3 * . The superscript, *, denotes the complex conjugate transforms.
For simplicity of notation, Equation (2) will be expressed in a shortened notation as:
f b = 1 B k = 0 B 1 F k e 2 π i b k B .
The shortened notation will be used hereafter. For given values of the ODF on the grid over B, the FFTs are:
F k = k = 0 B 1 f b e 2 π i b k B .
Computing the ODF transforms using the FFT methods was over two orders of magnitude faster compared to that obtained using the GSH methods. Benefits of the FFT method over the GSH method increase with the number of ODFs considered in the calculations of property closures because the computational time involved in ODF-to-elastic property calculations scales linearly with the number of ODFs.
Due to hexagonal crystal symmetry, there are 24 physically indistinctive orientations over the discretized Bunge-Euler space defined above, B. FZ is the fundamental zone within B containing the complete set of all physically distinct orientations that can occur in the sample. These redundancies within B can be exploited in the computations. To this end, we calculate function values in FZ and translate these throughout the rest of the space B while preserving the uniform grid. The translations are done as follows. For crystals of any symmetry, the definitions of the Bunge-Euler angles require that locations ( ϕ 1 + π ,   2 π Φ ,   ϕ 2 + π ) correspond to the exact same crystal lattice orientation. Table 1 shows twelve equivalent crystal orientations within the rotational point symmetry group associated with HCP crystals.
Table 1. Twelve symmetry operations within hexagonal close packed (HCP) rotational point symmetry group.
Table 1. Twelve symmetry operations within hexagonal close packed (HCP) rotational point symmetry group.
Symmetries 1–4Symmetries 5–8Symmetries 9–12
( ϕ 1 , Φ , ϕ 2 ) ( ϕ 1 , Φ , ϕ 2 + π 3 ) ( ϕ 1 , Φ , ϕ 2 + 2 π 3 )
( ϕ 1 , Φ , ϕ 2 + π ) ( ϕ 1 , Φ , ϕ 2 + 4 π 3 ) ( ϕ 1 , Φ , ϕ 2 + 5 π 3 )
( ϕ 1 , Φ + π , ϕ 2 ) ( ϕ 1 , Φ + π , ϕ 2 + π 3 ) ( ϕ 1 , Φ + π , ϕ 2 + 2 π 3 )
( ϕ 1 , Φ + π , ϕ 2 + π ) ( ϕ 1 , Φ + π , ϕ 2 + 4 π 3 ) ( ϕ 1 , Φ + π , ϕ 2 + 5 π 3 )
As a result of all of the considerations described above, the FZ space for HCP is F Z = { g = ( ϕ 1 , Φ , ϕ 2 )   |   0 ϕ 1 2 π , 0 Φ π 2 , 0 ϕ 2 π 3 } . Values of any desired function can be calculated on a uniform grid in FZ and subsequently extended over the entire periodic space using the above defined equivalencies preserving the uniform grid over B.
It should be noted that the smallest space defined earlier for cubic metals was defined as FZ3 [52], which is three times the actual FZ for cubic crystals. The reason for using FZ3 instead of FZ is that a uniform grid is only possible in FZ3. FZ upon being subjected to three-fold <111> symmetry operation, does not result in a uniform grid that can be used to fill out the remainder of the B space [8]. Therefore, the values of all desired functions on a uniform grid span the FZ3 space for cubic crystals while they conveniently span the FZ domain for HCP crystals.
The FFT representation allows defining the texture hull for HCP crystals. Let F k n define the FFTs of single crystals that lie within FZ, which are computed using Equation (4). Let N enumerate the single crystals that are distributed over in FZ. It is then possible to define a convex and compact texture hull [5] as:
M = { F k   |   F k = n = 1 N α n F k n ,   α n 0 ,   n = 1 N α n = 1 }
As mentioned earlier, M represents the complete set of all theoretically feasible ODFs including all of their linear combinations. M is defined in n-dimensional FFT space, where the number of dimensions in this space is equal to the discretization B. Selected projections of M are presented in Figure 1. The FFT based hulls are analogous to those presented previously using the GSH representations [7]. In the next section, we describe the methodology for mapping the texture hull in the material property space and construction of HCP property closures.
Figure 1. Examples of the fast Fourier transforms (FFT) texture hull for HCP-triclinic materials. Two arbitrarily selected normalized projections are: (a) ReF004, ImF004, ImF041 and (b) ReF444, F200, ImF444.
Figure 1. Examples of the fast Fourier transforms (FFT) texture hull for HCP-triclinic materials. Two arbitrarily selected normalized projections are: (a) ReF004, ImF004, ImF041 and (b) ReF444, F200, ImF444.
Materials 08 05303 g001

3. Property Closures

In previous work [7], it was demonstrated that a spectral approach based on the GSH basis can be used to construct property closures for HCP metals. Here, we show that the elastic property closures can be computed much more efficiently using the more readily accessible FFT representation. The primary focus here continues to be on the crystallographic texture in the sample as primary description of the microstructure. All of the closures reported in prior studies on HCP employed the orthotropic GSH functions, which implies orthotropic sample symmetry. The closures presented here using the FFT representations circumvent this symmetry assumption.

3.1. Elastic Stiffness for HCP Metals

The components of the stiffness tensor for an HCP single crystal as a function of the local crystal orientation and the five fundamental elastic constants C 11 ,   C 12 , C 13 , C 33 , and C 44 expressed in the sample frame are:
C a b c d = C 12 δ a b δ c d + C 44 ( δ a d δ b c + δ a c δ b d )        + ( C 11 C 12 2 C 44 )   ( t = 1 3 g a t g b t g c t g d t        + 1 2 ( g a 1 g b 2 g c 1 g d 2 + g a 1 g b 2 g c 2 g d 1 + g a 2 g b 1 g c 1 g d 2 + g a 2 g b 1 g c 2 g d 1 ) )        + ( C 13 C 12 ) t = 1 2 ( g a t g b t g c 3 g d 3 + g a 3 g b 3 g c t g d t ) + ( C 33 C 11 ) g a 3 g b 3 g c 3 g d 3
where δ a b represents the Kronecker symbol and g a b is:
g a b = [ cos ϕ 1 cos ϕ 2 sin ϕ 1 cos Φ sin ϕ 2 cos ϕ 1 sin ϕ 2 sin ϕ 1 cos Φ sin ϕ 2 sin ϕ 1 sin Φ sin ϕ 1 cos ϕ 2 + cos ϕ 1 cos Φ sin ϕ 2 sin ϕ 1 sin ϕ 2 + cos ϕ 1 cos Φ cos ϕ 2 cos ϕ 1 sin Φ sin Φ sin ϕ 2 sin Φ cos ϕ 2 cos Φ ]
We use bold letters to denote tensors. It is convenient to separate the orientation dependence of the tensor [7] as:
A ( g ) = A a b c d ( g ) = t = 1 3 g a t g b t g c t g d t + 1 2 ( g a 1 g b 2 g c 1 g d 2 + g a 1 g b 2 g c 2 g d 1 + g a 2 g b 1 g c 1 g d 2 + g a 2 g b 1 g c 2 g d 1 )
B ( g ) = B a b c d ( g ) = t = 1 2 ( g a t g b t g c 3 g d 3 + g a 3 g b 3 g c t g d t )
D ( g ) = D a b c d ( g ) =   g a 3 g b 3 g c 3 g d 3
Equation (6) then becomes:
C = C 12 δ a b δ c d + C 44 ( δ a d δ b c + δ a c δ b d ) + ( C 11 C 12 2 C 44 ) A ( g )   + ( C 13 C 12 ) B ( g ) + ( C 33 C 11 ) D ( g )

3.2. Representation of the Elastic Stiffness for HCP Metals Using FFTs

Orientation dependent functions A a b c d ( g ) , B a b c d ( g ) , and D a b c d ( g ) can be evaluated on a uniform grid over B and represented using FFTs as:
P ˜ k A = k = 0 B 1 A b s i n Φ e 2 π i b k B , P ˜ k B = k = 0 B 1 B b s i n Φ e 2 π i b k B , P ˜ k D = k = 0 B 1 D b s i n Φ e 2 π i b k B
To evaluate the property transforms, we had to evaluate the functions (Equation (8)) on a regular grid over B. The grid was refined to capture all of the non-zero frequencies present in the functions. A grid with one degree grid point spacing placed over B was found to satisfy this requirement, meaning that all of the important frequencies were captured by the transform. The non-zero Fourier transforms for elastic stiffness components are presented in Appendix A. Note that these transforms are orientation invariant and valid for all HCP metals. The database contains only half of the transforms since P ˜ k A = P ˜ * B k A , P ˜ k B = P ˜ * B k B ,   and   P ˜ k D = P ˜ * B k D . Note that for example P ˜ k A is a compact way of writing P ˜ k A a b c d = P ˜ k A . The tilde indicates that the functions have been pre multiplied by s i n Φ term. The term will be clarified shortly. As can be seen in the database, the number of transforms varies for the different stiffness components. Comparing the GSH Fourier coefficients presented in [7] for the same functions, we note that the FFT representation has a few more terms but the FFTs are tremendously faster and readily available.

3.3. First-Order Elastic Stiffness Bounds

The FFT representations in Equation (10) permit the efficient computation of the property bounds for any given ODF from the HCP texture hull. The first-order upper and lower bounds of effective elastic stiffness for the diagonal components are expressed as [12,18,19,20]:
( S ¯ 1 ) a b a b C a b a b * C ¯ a b a b
and for the off-diagonal components, the bounds are expressed as:
max ( C ¯ a b c d , ( S ¯ 1 ) a b c d ) Δ a b a b Δ c d c d C a b c d * min ( C ¯ a b c d , ( S ¯ 1 ) a b c d ) + Δ a b a b Δ c d c d
Δ a b c d = C ¯ a b c d ( S ¯ 1 ) a b c d
In Equations (11) and (12), no implicit summation on repeated indices is used (The Einstein indicial notation of implicit summation on repeated indices is employed in this paper, except when explicitly noted otherwise). The bars on top of a field quantity denote the volume averaged value over constituent crystals in a polycrystal: ( S ¯ 1 ) a b a b is the polycrystal compliance matrix with an inverse taken after homogenization and C ¯ a b a b is the polycrystal stiffness matrix. This paper is the first work to compute off-diagonal term HCP property closures.

3.4. Homogenization of the Elastic Properties in Fourier Space

As mentioned in the previous section, the first-order homogenization theory requires calculations of the volume average quantities. As an example, the upper bound for the components of the homogenized elastic stiffness tensor can be expressed as the following (no implicit summation on repeated indices):
C ¯ = C ¯ a b c d = 1 V V C a b c d ( x ) d x = F Z 3 f ( g ) C a b c d ( g ) s i n Φ d ϕ 1 d Φ d ϕ 2 = C 12 δ a b δ c d + C 44 ( δ a d δ b c + δ a c δ b d ) + 1 B k = 0 B 1 [ ( C 11 C 12 2 C 44 ) P ˜ k A + ( C 13 C 12 ) P ˜ k B + ( C 33 C 11 ) P ˜ k D ] F k
where V is the physical volume of a sample and x is a location in the sample. Alternatively, the volume average functions A, B and D can be calculated by multiplication of the precomputed P ˜ k A , P ˜ k B , and P ˜ k D with the ODF transforms Fk and their summation as: A ¯ a b c d = 1 B k = 0 B 1 P ˜ k A F k , B ¯ a b c d = 1 B k = 0 B 1 P ˜ k B F k , D ¯ a b c d = 1 B k = 0 B 1 P ˜ k D F k . These quantities can be used in Equation (9) to calculate the homogenized elastic stiffness tensor. Similar expressions exist for all tensors involved in Equations (11) and (12). Note that we have exploited the orthogonal properties of the spectral representation in the last part of the Equation (13). Since the regular grid is used in the summations in each direction incrementally ( Δ ϕ 1 , Δ Φ , Δ ϕ 2 ) , the s i n Φ accounting for the distortion of the Bunge-Euler space must multiply either f ( g ) or C a b c d ( g ) (Equation (13)). We have chosen to multiply the functions A, B and D by the s i n Φ term, therefore we have added tilde symbol to the property transforms.
Equation (13) can be used to calculate the volume average elastic stiffness as well as the elastic stiffness of a single crystal. It should be noted that the number of calculations (multiplications and summations) for stiffness components is the same for a single crystal and a polycrystal comprised of any number of crystals in this spectral representation. The number of calculations is determined by the number of non-zero transforms (Appendix A). Figure 2 shows contour plots of the C 1111 ( g ) elastic stiffness component over FZ for HCP metals calculated using Equations (6) and (13). The single crystal elastic constants of Zn were used. It can be clearly seen that the results are identical, which validates the FFT representation of the ODF-to-elastic stiffness relationship. The same checks were conducted for all components of the local elastic stiffness tensor.
Figure 2. Contour plots of the C 1111 elastic stiffness component in the HCP fundamental zone of the Bunge-Euler space for Zn: (a) computed using the non-zero FFTs (Equation (13)) and (b) computed directly using Equation (6). The maximum difference between corresponding locations in the two plots is of the order 10−6 GPa.
Figure 2. Contour plots of the C 1111 elastic stiffness component in the HCP fundamental zone of the Bunge-Euler space for Zn: (a) computed using the non-zero FFTs (Equation (13)) and (b) computed directly using Equation (6). The maximum difference between corresponding locations in the two plots is of the order 10−6 GPa.
Materials 08 05303 g002

3.5. Computation of Property Closures for HCP Metals

The methodology used here for building first-order elastic stiffness closures starts with a consideration of a set of points in the texture hull that correspond to “eigen textures” [66]. The set of eigen textures is selected while ensuring an adequate coverage of the HCP FZ. The property bounds (Equations (11) and (12)) are first evaluated using the FFT representation of the ODF-to-elastic stiffness linkages (Equation (13)) for these eigen textures. A finite number of eigen textures corresponding to the boundary of a given property closure of stiffness components was then selected. Subsequently, the property combinations were evaluated for weighted combinations of these textures, taking one pair of textures at a time. The weighted combinations of pairs of the selected textures were incremented by a 0.1 weighted fraction i.e., from (0.1, 0.9) pair to (0.9, 0.1) pair, which resulted with nine weighted combinations. As expected, these computations involving combinations of textures expand the property closures. Next, a new set of textures corresponding to the new boundary of the expanded closure were selected (this time these were a mixture of eigen textures and non-eigen textures) and the property combinations corresponding to the weighted combinations of these were evaluated. This process was repeated until the closure expansion saturated. The method of delineation of closures follows the main ideas underlying genetic algorithms, where good solutions are pre-selected and used in further calculations. This approach was applied successfully to delineate closures of plastic properties involving uniform ductility and ultimate tensile strength in cubic metals in [15].

3.6. Atlases of Property Closures for HCP Metals

The FFT representation of the ODF-to-elastic stiffness bounds is used in this study to obtain the first-order property closures, which identify the complete set of theoretically realizable combinations of selected macroscale properties in a given material system through a consideration of the complete set of textures (i.e., the texture hull). The shaded areas inside the closures represent the possible property combinations of the selected elastic stiffness components that can be obtained according to the first-order bounding theories for the particular material as a function of ODFs. We provide three examples of atlases of property closures corresponding to selection of different pairs of effective elastic properties. Figure 3 depicts examples of the property closures produced in this work using the FFT representations for a range of HCP metals. This particular property combination of the effective modulus in unaxial strain and the effective shear modulus in the sample play an important role in the design of components subjected simultaneously to axial loads and twisting moments. The single crystal elastic constants for these metals were taken from literature [67]. Comparison with the closures previously reported in [7], which were computed using the GSH representations reveals excellent agreement. It should be noted that this agreement occurs in spite of the fact that the previously computed closures assumed orthotropic sample symmetry, whereas the computations reported here did not invoke any sample symmetry assumption. This is because the elastic stiffness components involved in these two closures exhibit orthotropic sample symmetry. Figure 4 depicts another set of closures for the same set of metals that involve the coupling of the normal-shear components. When orthotropic sample symmetry is assumed, C 1112 is zero. The closures presented here in Figure 4 show that the values of C 1112 can take a wide range, which can play an important role in the macroscale mechanical response of components (note that these are often ignored in practice by assuming an orthotropic sample symmetry that the material may not exhibit).
As a final comment, we note that the FFT framework developed here for the efficient calculation of elastic properties can also be used to accelerate the calculations of the elasto-viscoplastic response of polycrystalline HCP metals that is under development. This framework has been recently reported for cubic metals in [68]. The implementation presented in [68] takes the advantage of calculations of the elastic properties based on the spectral representation reported in [52]. For HCP, the components of the local elastic stiffness tensor as a function of the crystal orientation can be calculated using Equation (13) instead of Equation (6). Therefore, the elastic stiffness solutions obtained based on the FFT representation can be used within the elasto-viscoplastic crystal plasticity constitutive laws operating at every integration point within implicit finite elements [69,70,71,72,73,74,75].
Figure 3. Atlases of property closures for a range of HCP metals computed using the FFT method: (a) C1111 vs. C1212 and (b) C1111 vs. C1122.
Figure 3. Atlases of property closures for a range of HCP metals computed using the FFT method: (a) C1111 vs. C1212 and (b) C1111 vs. C1122.
Materials 08 05303 g003
Figure 4. Atlases of the normal-shear coupling stiffness closures for a range of HCP metals computed using the FFT method.
Figure 4. Atlases of the normal-shear coupling stiffness closures for a range of HCP metals computed using the FFT method.
Materials 08 05303 g004

4. Conclusions

In this paper, we have extended the recently developed first-order ODF-to-elastic stiffness linkages using FFTs for cubic metals to more computationally challenging HCP polycrystals. The benefit of using FFTs to represent these linkages is the substantial reduction in the computational time involved in the property calculations due to the use of fast and readily available FFT algorithms. First, we have shown that it is possible to construct texture hulls using FFTs for HCP polycrystals, as was done earlier for cubic polycrystals. Next, we have calculated elastic stiffness FFTs and presented the actual values in Appendix A. Consistent with the GSH representation, we found that the inherently more anisotropic HCP structure demands a larger number of FFTs than FFTs used for the cubic structure to represent these properties. However, the use of FFT representations reduced the computational time involved in producing the property closures due to the use of efficient FFT algorithms. Finally, we have successfully computed elastic property closures for HCP metals without invoking any sample symmetry assumption. To this end, the methodology presented here facilitated for the first time the delineation of the property closures involving the normal–shear coupling stiffness coefficients.

Acknowledgments

This work is based upon project supported by the National Science Foundation under grant No. CMMI-1541918. The authors gratefully acknowledge this support.

Author Contributions

Marko Knezevic defined the research described in this paper and advised Nicholas W. Landry who performed the work. The authors wrote the paper through a typical student-advisor relationship.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

P ˜ k 1 , k 2 , k 3 A
P ˜ k 1 , k 2 , k 3 A
k1 k2 k3C1111k1 k2 k3C1112k1 k2 k3C1113
0 0 03.645 × 107 + 0j0 0 01.346 × 10−3 + 0j0 0 03.474 × 10−4 + 0j
0 0 22.916 × 106 − 8.017 × 10−4j0 0 2−1.024 × 10−3 − 1.458 × 106j0 2 11.458 × 106 − 1.752 × 10−3j
0 0 42.187 × 106 + 1.490 × 10−3j0 0 4−3.287 × 10−4 − 2.187 × 106j0 2 31.458 × 106 − 9.201 × 10−4j
0 2 02.916 × 106 − 8.017 × 10−4j0 2 45.437 × 10−4 + 1.458 × 106j0 2 357−1.458 × 106 + 8.970 × 10−4j
0 2 4−1.458 × 106 − 2.353 × 10−4j0 2 3563.103 × 10−4 − 1.458 × 106j0 2 359−1.458 × 106 + 1.470 × 10−3j
0 2 356−1.458 × 106 − 1.324 × 10−3j0 4 25.819 × 10−4 + 7.290 × 105j0 4 12.187 × 106 − 3.252 × 10−3j
0 4 02.187 × 106 + 1.490 × 10−3j0 4 4−1.197 × 10−3 − 3.645 × 105j0 4 3−7.290 × 105 + 7.335 × 10−4j
0 4 2−1.458 × 106 − 2.353 × 10−4j0 4 356−3.065 × 10−4 + 3.645 × 105j0 4 3577.290 × 105 + 1.229 × 10−3j
0 4 43.645 × 105 − 4.225 × 10−4j0 4 358−1.105 × 10−3 − 7.290 × 105j0 4 359−2.187 × 106 + 1.318 × 10−3j
0 4 3563.645 × 105 − 2.071 × 10−12j
0 4 358−1.458 × 106 + 1.324 × 10−3j
k1 k2 k3C1122k1 k2 k3C1123k1 k2 k3C1133
0 0 04.374 × 106 + 0j0 0 0−8.664 × 10−4 + 0j0 0 05.832 × 106 + 0j
0 0 4−2.187 × 106 + 7.111 × 10−4j0 2 1−3.505 × 10−3 + 1.458 × 106j0 0 2−2.916 × 106 + 1.366 × 10−5j
0 2 0−2.916 × 106 + 3.868 × 10−3j0 2 32.036 × 10−3 − 1.458 × 106j0 4 0−2.916 × 106 + 3.465 × 10−3j
0 2 41.458 × 106 − 2.976 × 10−3j0 2 357−1.196 × 10−3 − 1.458 × 106j0 4 21.458 × 106 − 2.873 × 10−3j
0 2 3561.458 × 106 − 2.763 × 10−3j0 2 3592.458 × 10−3 + 1.458 × 106j0 4 3581.458 × 106 − 2.842 × 10−3j
0 4 07.290 × 105 − 4.308 × 10−3j0 4 14.336 × 10−3 − 7.290 × 105j
0 4 4−3.645 × 105 + 2.327 × 10−3j0 4 3−3.380 × 10−3 + 7.290 × 105j
0 4 356−3.645 × 105 + 3.997 × 10−3j0 4 3572.269 × 10−4 + 7.290 × 105j
0 4 359−2.354 × 10−3 − 7.290 × 105j
k1 k2 k3C1212k1 k2 k3C1312k1 k2 k3C1313
0 0 01.604 × 107 + 0j0 0 05.001 × 10−4 + 0j0 0 01.166 × 107 + 0j
0 0 4−2.187 × 106 + 1.968 × 10−3j0 2 1−2.863 × 10−3 − 1.458 × 106j0 2 0−2.916 × 106 + 1.449 × 10−3j
0 2 02.916 × 106 − 3.666 × 10−4j0 2 31.106 × 10−4 − 1.458 × 106j0 2 2−1.458 × 106 − 2.397 × 10−3j
0 2 41.458 × 106 − 3.097 × 10−3j0 2 357−1.002 × 10−3 − 1.458 × 106j0 2 358−1.458 × 106 − 2.518 × 10−3j
0 2 3561.458 × 106 − 2.027 × 10−3j0 2 359−1.417 × 10−3 − 1.458 × 106j0 4 0−2.916 × 106 + 2.298 × 10−4j
0 4 07.290 × 105 + 3.823 × 10−3j0 4 1−7.437 × 10−4 − 7.290 × 105j0 4 21.458 × 106 + 1.402 × 10−3j
0 4 4−3.645 × 105 + 1.073 × 10−3j0 4 3−6.925 × 10−4 + 7.290 × 105j0 4 3581.458 × 106 + 5.424 × 10−4j
0 4 356−3.645 × 105 + 3.916 × 10−3j0 4 357−3.324 × 10−4 + 7.290 × 105j
0 4 359−8.323 × 10−5 − 7.290 × 105j
k1 k2 k3C2212k1 k2 k3C2213k1 k2 k3C2222
0 0 0−2.227 × 10−3 + 0j0 0 01.855 × 10−3 + 0j0 0 03.645 × 107 + 0j
0 0 2−8.014 × 10−4 − 1.458 × 106j0 2 1−1.458 × 106 + 2.242 × 10−3j0 0 2−2.916 × 106 + 2.652 × 10−3j
0 0 41.646 × 10−3 + 2.187 × 106j0 2 3−1.458 × 106 − 9.985 × 10−4j0 0 42.187 × 106 − 2.554 × 10−3j
0 2 4−1.432 × 10−3 − 1.458 × 106j0 2 3571.458 × 106 − 9.966 × 10−4j0 2 02.916 × 106 + 1.752 × 10−3j
0 2 356−1.473 × 10−3 + 1.458 × 106j0 2 3591.458 × 106 + 2.525 × 10−3j0 2 4−1.458 × 106 − 8.290 × 10−5j
0 4 21.079 × 10−4 + 7.290 × 105j0 4 17.290 × 105 − 9.776 × 10−4j0 2 356−1.458 × 106 − 3.054 × 10−3j
0 4 41.413 × 10−3 + 3.645 × 105j0 4 37.290 × 105 + 4.229 × 10−3j0 4 02.187 × 106 − 7.525 × 10−4j
0 4 3562.592 × 10−4 − 3.645 × 105j0 4 357−7.290 × 105 − 1.011 × 10−4j0 4 21.458 × 106 − 2.064 × 10−3j
0 4 358−4.406 × 10−4 − 7.290 × 105j0 4 359−7.290 × 105 − 3.638 × 10−3j0 4 43.645 × 105 − 2.895 × 10−4j
0 4 3563.645 × 105 + 8.692 × 10−4j
0 4 3581.458 × 106 + 1.344 × 10−3j
k1 k2 k3C2223k1 k2 k3C2233k1 k2 k3C2312
0 0 0−3.273 × 10−4 + 0j0 0 05.832 × 106 + 0j0 0 0−5.371 × 10−4 + 0j
0 2 16.770 × 10−5 − 1.458 × 106j0 0 22.916 × 106 − 8.777 × 10−4j0 2 11.458 × 106 − 9.363 × 10−4j
0 2 31.475 × 10−3 + 1.458 × 106j0 4 0−2.916 × 106 + 4.489 × 10−3j0 2 3−1.458 × 106 − 4.469 × 10−4j
0 2 3571.473 × 10−3 + 1.458 × 106j0 4 2−1.458 × 106 + 3.561 × 10−3j0 2 3571.458 × 106 − 3.113 × 10−4j
0 2 359−1.745 × 10−3 − 1.458 × 106j0 4 358−1.458 × 106 + 2.759 × 10−3j0 2 359−1.458 × 106 + 2.436 × 10−3j
0 4 1−2.036 × 10−3 − 2.187 × 106j0 4 17.290 × 105 − 3.767 × 10−4j
0 4 34.109 × 10−4 − 7.290 × 105j0 4 37.290 × 105 + 7.465 × 10−4j
0 4 3571.317 × 10−4 − 7.290 × 105j0 4 357−7.290 × 105 − 4.772 × 10−4j
0 4 359−1.982 × 10−3 − 2.187 × 106j0 4 359−7.290 × 105 + 6.635 × 10−4j
k1 k2 k3C2313k1 k2 k3C2323k1 k2 k3C3312
0 0 0−4.698 × 10−4 + 0j0 0 01.166 × 107 + 0j0 0 05.534 × 10−4 + 0j
0 2 22.663 × 10−3 + 1.458 × 106j0 2 0−2.916 × 106 + 3.338 × 10−3j0 0 22.781 × 10−4 + 2.916 × 106j
0 2 358−1.665 × 10−3 − 1.458 × 106j0 2 21.458 × 106 + 2.699 × 10−3j0 4 2−6.452 × 10−3 − 1.458 × 106j
0 4 2−3.239 × 10−3 − 1.458 × 106j0 2 3581.458 × 106 + 3.103 × 10−3j0 4 3585.710 × 10−3 + 1.458 × 106j
0 4 3582.674 × 10−3 + 1.458 × 106j0 4 0−2.916 × 106 + 1.859 × 10−3j
0 4 2−1.458 × 106 + 6.743 × 10−4j
0 4 358−1.458 × 106 − 9.184 × 10−4j
k1 k2 k3C3313k1 k2 k3C3323k1 k2 k3C3333
0 0 03.737 × 10−4 + 0j0 0 03.272 × 10−4 + 0j0 0 03.499 × 107 + 0j
0 4 1−2.916 × 106 + 1.377 × 10−3j0 4 12.203 × 10−3 + 2.916 × 106j0 4 05.832 × 106 + 2.951 × 10−11j
0 4 3592.916 × 106 − 2.212 × 10−3j0 4 3593.092 × 10−4 + 2.916 × 106j
P ˜ k 1 , k 2 , k 3 B
P ˜ k 1 , k 2 , k 3 B
k1 k2 k3C1111k1 k2 k3C1112k1 k2 k3C1113
0 0 01.021 × 107 + 0j0 0 0−1.346 × 10−3 + 0j0 0 0−3.474 × 10−4 + 0j
0 0 2−2.916 × 106 + 4.497 × 10−4j0 0 21.024 × 10−3 + 1.458 × 106j0 2 1−1.458 × 106 + 1.752 × 10−3j
0 0 4−2.187 × 106 − 3.121 × 10−3j0 0 43.287 × 10−4 + 2.187 × 106j0 2 3−1.458 × 106 + 9.201 × 10−4j
0 2 0−2.916 × 106 + 4.497 × 10−4j0 2 4−5.437 × 10−4 − 1.458 × 106j0 2 3571.458 × 106 − 8.970 × 10−4j
0 2 41.458 × 106 + 1.246 × 10−3j0 2 356−3.103 × 10−4 + 1.458 × 106j0 2 3591.458 × 106 − 1.470 × 10−3j
0 2 3561.458 × 106 + 5.795 × 10−4j0 4 2−5.819 × 10−4 − 7.290 × 105j0 4 1−2.187 × 106 + 3.252 × 10−3j
0 4 0−2.187 × 106 − 3.121 × 10−3j0 4 41.197 × 10−3 + 3.645 × 105j0 4 37.290 × 105 − 7.335 × 10−4j
0 4 21.458 × 106 + 1.246 × 10−3j0 4 3563.065 × 10−4 − 3.645 × 105j0 4 357−7.290 × 105 − 1.229 × 10−3j
0 4 4−3.645 × 105 − 3.397 × 10−4j0 4 3581.105 × 10−3 + 7.290 × 105j0 4 3592.187 × 106 − 1.318 × 10−3j
0 4 356−3.645 × 105 + 7.221 × 10−12j
0 4 3581.458 × 106 − 5.795 × 10−4j
k1 k2 k3C1122k1 k2 k3C1123k1 k2 k3C1133
0 0 01.895 × 107 + 0j0 0 0−1.405 × 10−3 + 0j0 0 02.916 × 107 + 0j
0 0 42.187 × 106 + 2.054 × 10−3j0 2 13.890 × 10−3 + 4.374 × 106j0 0 2−2.916 × 106 + 2.518 × 10−3j
0 2 0−8.748 × 106 − 5.686 × 10−3j0 2 31.772 × 10−3 + 1.458 × 106j0 2 05.832 × 106 − 5.718 × 10−4j
0 2 4−1.458 × 106 − 1.370 × 10−4j0 2 3571.828 × 10−3 + 1.458 × 106j0 2 22.916 × 106 − 3.981 × 10−3j
0 2 356−1.458 × 106 + 2.013 × 10−3j0 2 3591.475 × 10−3 + 4.374 × 106j0 2 3582.916 × 106 + 4.465 × 10−4j
0 4 0−7.290 × 105 − 1.318 × 10−2j0 4 11.098 × 10−3 + 7.290 × 105j0 4 02.916 × 106 − 5.226 × 10−4j
0 4 43.645 × 105 + 5.098 × 10−3j0 4 33.059 × 10−3 − 7.290 × 105j0 4 2−1.458 × 106 + 9.562 × 10−3j
0 4 3563.645 × 105 + 4.815 × 10−3j0 4 357−1.042 × 10−3 − 7.290 × 105j0 4 358−1.458 × 106 + 7.009 × 10−3j
0 4 359−8.954 × 10−4 + 7.290 × 105j
k1 k2 k3C1212k1 k2 k3C1312k1 k2 k3C1313
0 0 0−4.374 × 106 + 0j0 0 08.664 × 10−4 + 0j0 0 0−5.832 × 106 + 0j
0 0 42.187 × 106 − 7.102 × 10−4j0 2 13.505 × 10−3 − 1.458 × 106j0 0 22.916 × 106 − 1.366 × 10−5j
0 2 02.916 × 106 − 3.869 × 10−3j0 2 3−2.036 × 10−3 + 1.458 × 106j0 4 02.916 × 106 − 3.465 × 10−3j
0 2 4−1.458 × 106 + 2.976 × 10−3j0 2 3571.196 × 10−3 + 1.458 × 106j0 4 2−1.458 × 106 + 2.873 × 10−3j
0 2 356−1.458 × 106 + 2.764 × 10−3j0 2 359−2.458 × 10−3 − 1.458 × 106j0 4 358−1.458 × 106 + 2.842 × 10−3j
0 4 0−7.290 × 105 + 4.309 × 10−3j0 4 1−4.336 × 10−3 + 7.290 × 105j
0 4 43.645 × 105 − 2.328 × 10−3j0 4 33.380 × 10−3 − 7.290 × 105j
0 4 3563.645 × 105 − 3.997 × 10−3j0 4 357−2.269 × 10−4 − 7.290 × 105j
0 4 3592.354 × 10−3 + 7.290 × 105j
k1 k2 k3C2212k1 k2 k3C2213k1 k2 k3C2222
0 0 02.227 × 10−3 + 0j0 0 01.553 × 10−3 + 0j0 0 01.021 × 107 + 0j
0 0 28.014 × 10−4 + 1.458 × 106j0 2 1−4.374 × 106 − 8.811 × 10−5j0 0 22.916 × 106 − 2.431 × 10−3j
0 0 4−1.646 × 10−3 − 2.187 × 106j0 2 31.458 × 106 − 2.859 × 10−3j0 0 4−2.187 × 106 + 2.792 × 10−3j
0 2 41.432 × 10−3 + 1.458 × 106j0 2 357−1.458 × 106 + 6.060 × 10−4j0 2 0−2.916 × 106 − 2.517 × 10−3j
0 2 3561.473 × 10−3 − 1.458 × 106j0 2 3594.374 × 106 + 1.237 × 10−4j0 2 41.458 × 106 + 1.090 × 10−3j
0 4 2−1.079 × 10−4 − 7.290 × 105j0 4 1−7.290 × 105 − 6.559 × 10−4j0 2 3561.458 × 106 + 1.432 × 10−3j
0 4 4−1.413 × 10−3 − 3.645 × 105j0 4 3−7.290 × 105 − 3.227 × 10−3j0 4 0−2.187 × 106 − 5.437 × 10−4j
0 4 356−2.592 × 10−4 + 3.645 × 105j0 4 3577.290 × 105 + 2.268 × 10−4j0 4 2−1.458 × 106 + 2.247 × 10−3j
0 4 3584.406 × 10−4 + 7.290 × 105j0 4 3597.290 × 105 + 3.400 × 10−3j0 4 4−3.645 × 105 − 1.823 × 10−3j
0 4 356−3.645 × 105 + 1.969 × 10−4j
0 4 358−1.458 × 106 − 3.514 × 10−3j
k1 k2 k3C2223k1 k2 k3C2233k1 k2 k3C2312
0 0 03.273 × 10−4 + 0j0 0 02.916 × 107 + 0j0 0 0−1.855 × 10−3 + 0j
0 2 1−6.770 × 10−5 + 1.458 × 106j0 0 22.916 × 106 + 2.189 × 10−3j0 2 11.458 × 106 − 2.242 × 10−3j
0 2 3−1.475 × 10−3 − 1.458 × 106j0 2 05.832 × 106 + 3.827 × 10−4j0 2 31.458 × 106 + 9.985 × 10−4j
0 2 357−1.473 × 10−3 − 1.458 × 106j0 2 2−2.916 × 106 − 8.557 × 10−4j0 2 357−1.458 × 106 + 9.966 × 10−4j
0 2 3591.745 × 10−3 + 1.458 × 106j0 2 358−2.916 × 106 + 2.756 × 10−3j0 2 359−1.458 × 106 − 2.525 × 10−3j
0 4 12.036 × 10−3 + 2.187 × 106j0 4 02.916 × 106 − 1.242 × 10−3j0 4 1−7.290 × 105 + 9.776 × 10−4j
0 4 3−4.109 × 10−4 + 7.290 × 105j0 4 21.458 × 106 − 5.018 × 10−3j0 4 3−7.290 × 105 − 4.229 × 10−3j
0 4 357−1.317 × 10−4 + 7.290 × 105j0 4 3581.458 × 106 − 7.378 × 10−3j0 4 3577.290 × 105 + 1.011 × 10−4j
0 4 3591.982 × 10−3 + 2.187 × 106j0 4 3597.290 × 105 + 3.638 × 10−3j
k1 k2 k3C2313k1 k2 k3C2323k1 k2 k3C3312
0 0 0−5.534 × 10−4 + 0j0 0 0−5.832 × 106 + 0j0 0 0−1.994 × 10−3 + 0j
0 0 2−2.781 × 10−4 − 2.916 × 106j0 0 2−2.916 × 106 + 8.777 × 10−4j0 0 22.216 × 10−3 + 2.916 × 106j
0 4 26.452 × 10−3 + 1.458 × 106j0 4 02.916 × 106 − 4.489 × 10−3j0 2 2−3.090 × 10−3 − 2.916 × 106j
0 4 358−5.710 × 10−3 − 1.458 × 106j0 4 21.458 × 106 − 3.561 × 10−3j0 2 358−1.514 × 10−3 + 2.916 × 106j
0 4 3581.458 × 106 − 2.759 × 10−3j0 4 23.003 × 10−3 + 1.458 × 106j
0 4 358−1.255 × 10−4 − 1.458 × 106j
k1 k2 k3C3313k1 k2 k3C3323k1 k2 k3C3333
0 0 0−3.737 × 10−4 + 0j0 0 0−3.272 × 10−4 + 0j0 0 01.166 × 107 + 0j
0 4 12.916 × 106 − 1.377 × 10−3j0 4 1−2.203 × 10−3 − 2.916 × 106j0 4 0−5.832 × 106 − 1.367 × 10−3j
0 4 359−2.916 × 106 + 2.212 × 10−3j0 4 359−3.092 × 10−4 − 2.916 × 106j
P ˜ k 1 , k 2 , k 3 D
P ˜ k 1 , k 2 , k 3 D
k1 k2 k3C1111k1 k2 k3C1112k1 k2 k3C1113
0 0 06.561 × 106 + 0j0 0 03.823 × 10−3 + 0j0 0 0−1.282 × 10−3 + 0j
0 0 2−4.374 × 106 + 1.160 × 10−3j0 0 2−2.355 × 10−3 + 2.187 × 106j0 2 1−2.187 × 106 + 1.105 × 10−3j
0 0 41.093 × 106 − 2.072 × 10−3j0 0 41.071 × 10−3 − 1.093 × 106j0 2 37.290 × 105 − 2.142 × 10−3j
0 2 0−4.374 × 106 + 1.160 × 10−3j0 2 26.666 × 10−4 − 1.458 × 106j0 2 357−7.290 × 105 − 2.591 × 10−3j
0 2 22.916 × 106 − 1.262 × 10−3j0 2 43.187 × 10−4 + 7.290 × 105j0 2 3592.187 × 106 + 1.453 × 10−3j
0 2 4−7.290 × 105 + 1.499 × 10−3j0 2 356−1.668 × 10−3 − 7.290 × 105j0 4 11.093 × 106 − 1.071 × 10−3j
0 2 356−7.290 × 105 − 6.593 × 10−4j0 2 3582.999 × 10−3 + 1.458 × 106j0 4 3−3.645 × 105 + 2.670 × 10−3j
0 2 3582.916 × 106 − 2.659 × 10−10j0 4 2−3.242 × 10−4 + 3.645 × 105j0 4 3573.645 × 105 + 2.440 × 10−3j
0 4 01.093 × 106 − 2.072 × 10−3j0 4 4−6.012 × 10−4 − 1.823 × 105j0 4 359−1.093 × 106 − 1.311 × 10−3j
0 4 2−7.290 × 105 + 1.499 × 10−3j0 4 3568.978 × 10−4 + 1.823 × 105j
0 4 41.823 × 105 + 9.732 × 10−5j0 4 358−1.676 × 10−3 − 3.645 × 105j
0 4 3561.823 × 105 + 9.683 × 10−12j
0 4 358−7.290 × 105 + 6.593 × 10−4j
k1 k2 k3C1122k1 k2 k3C1123k1 k2 k3C1133
0 0 02.187 × 106 + 0j0 0 0−2.897 × 10−4 + 0j0 0 02.916 × 106 + 0j
0 0 4−1.094 × 106 + 8.084 × 10−4j0 2 1−4.025 × 10−4 + 7.290 × 105j0 0 2−1.458 × 106 − 1.382 × 10−4j
0 2 0−1.458 × 106 + 1.200 × 10−3j0 2 32.409 × 10−4 − 7.290 × 105j0 4 0−1.458 × 106 − 6.165 × 10−6j
0 2 47.290 × 105 − 1.126 × 10−3j0 2 357−4.364 × 10−4 − 7.290 × 105j0 4 27.290 × 105 + 4.241 × 10−4j
0 2 3567.290 × 105 − 1.900 × 10−4j0 2 3591.795 × 10−4 + 7.290 × 105j0 4 3587.290 × 105 + 1.654 × 10−4j
0 4 03.645 × 105 − 8.533 × 10−4j0 4 18.402 × 10−4 − 3.645 × 105j
0 4 4−1.823 × 105 + 2.970 × 10−4j0 4 3−7.023 × 10−4 + 3.645 × 105j
0 4 356−1.823 × 105 + 5.990 × 10−4j0 4 357−1.036 × 10−3 + 3.645 × 105j
0 4 3594.759 × 10−4 − 3.645 × 105j
k1 k2 k3C1212k1 k2 k3C1312k1 k2 k3C1313
0 0 02.187 × 106 + 0j0 0 0−2.897 × 10−4 + 0j0 0 02.916 × 106 + 0j
0 0 4−1.094 × 106 + 8.084 × 10−4j0 2 1−4.025 × 10−4 + 7.290 × 105j0 0 2−1.458 × 106 − 1.382 × 10−4j
0 2 0−1.458 × 106 + 1.200 × 10−3j0 2 32.409 × 10−4 − 7.290 × 105j0 4 0−1.458 × 106 − 6.165 × 10−6j
0 2 47.290 × 105 − 1.126 × 10−3j0 2 357−4.364 × 10−4 − 7.290 × 105j0 4 27.290 × 105 + 4.241 × 10−4j
0 2 3567.290 × 105 − 1.900 × 10−4j0 2 3591.795 × 10−4 + 7.290 × 105j0 4 3587.290 × 105 + 1.654 × 10−4j
0 4 03.645 × 105 − 8.533 × 10−4j0 4 18.402 × 10−4 − 3.645 × 105j
0 4 4−1.823 × 105 + 2.970 × 10−4j0 4 3−7.023 × 10−4 + 3.645 × 105j
0 4 356−1.823 × 105 + 5.990 × 10−4j0 4 357−1.036 × 10−3 + 3.645 × 105j
0 4 3594.759 × 10−4 − 3.645 × 105j
k1 k2 k3C2212k1 k2 k3C2213k1 k2 k3C2222
0 0 02.792 × 10−3 + 0j0 0 03.074 × 10−5 + 0j0 0 06.561 × 106 + 0j
0 0 22.951 × 10−3 + 2.187 × 106j0 2 1−7.290 × 105 + 2.330 × 10−4j0 0 24.374 × 106 + 7.014 × 10−4j
0 0 42.232 × 10−3 + 1.093 × 106j0 2 3−7.290 × 105 + 5.117 × 10−4j0 0 41.093 × 106 + 1.030 × 10−3j
0 2 2−3.103 × 10−3 − 1.458 × 106j0 2 3577.290 × 105 + 5.010 × 10−4j0 2 0−4.374 × 106 + 9.433 × 10−4j
0 2 4−2.957 × 10−3 − 7.290 × 105j0 2 3597.290 × 105 + 5.058 × 10−4j0 2 2−2.916 × 106 + 1.423 × 10−4j
0 2 356−5.481 × 10−4 + 7.290 × 105j0 4 13.645 × 105 − 1.666 × 10−4j0 2 4−7.290 × 105 − 1.449 × 10−4j
0 2 358−1.477 × 10−3 + 1.458 × 106j0 4 33.645 × 105 + 1.777 × 10−4j0 2 356−7.290 × 105 + 8.615 × 10−4j
0 4 21.298 × 10−3 + 3.645 × 105j0 4 357−3.645 × 105 − 3.247 × 10−4j0 2 358−2.916 × 106 + 1.174 × 10−3j
0 4 41.733 × 10−3 + 1.823 × 105j0 4 359−3.645 × 105 − 1.046 × 10−3j0 4 01.093 × 106 − 2.872 × 10−3j
0 4 3562.726 × 10−4 − 1.823 × 105j0 4 27.290 × 105 − 1.485 × 10−3j
0 4 3589.600 × 10−4 − 3.645 × 105j0 4 41.823 × 105 − 3.610 × 10−4j
0 4 3561.823 × 105 + 5.363 × 10−4j
0 4 3587.290 × 105 − 2.104 × 10−3j
k1 k2 k3C2223k1 k2 k3C2233k1 k2 k3C2312
0 0 01.254 × 10−3 + 0j0 0 02.916 × 106 + 0j0 0 03.074 × 10−5 + 0j
0 2 1−1.485 × 10−3 + 2.187 × 106j0 0 21.458 × 106 − 4.806 × 10−4j0 2 1−7.290 × 105 + 2.330 × 10−4j
0 2 3−2.290 × 10−3 + 7.290 × 105j0 4 0−1.458 × 106 + 1.795 × 10−4j0 2 3−7.290 × 105 + 5.117 × 10−4j
0 2 3571.152 × 10−3 + 7.290 × 105j0 4 2−7.290 × 105 + 2.234 × 10−4j0 2 3577.290 × 105 + 5.010 × 10−4j
0 2 3593.268 × 10−4 + 2.187 × 106j0 4 358−7.290 × 105 − 4.228 × 10−4j0 2 3597.290 × 105 + 5.058 × 10−4j
0 4 16.655 × 10−4 − 1.093 × 106j0 4 13.645 × 105 − 1.666 × 10−4j
0 4 31.517 × 10−3 − 3.645 × 105j0 4 33.645 × 105 + 1.777 × 10−4j
0 4 357−1.949 × 10−3 − 3.645 × 105j0 4 357−3.645 × 105 − 3.247 × 10−4j
0 4 359−8.841 × 10−4 − 1.093 × 106j0 4 359−3.645 × 105 − 1.046 × 10−3j
k1 k2 k3C2313k1 k2 k3C2323k1 k2 k3C3312
0 0 0−1.640 × 10−3 + 0j0 0 02.916 × 106 + 0j0 0 0−1.640 × 10−3 + 0j
0 0 21.228 × 10−3 + 1.458 × 106j0 0 21.458 × 106 − 4.806 × 10−4j0 0 21.228 × 10−3 + 1.458 × 106j
0 4 2−1.516 × 10−3 − 7.290 × 105j0 4 0−1.458 × 106 + 1.795 × 10−4j0 4 2−1.516 × 10−3 − 7.290 × 105j
0 4 358−4.549 × 10−5 + 7.290 × 105j0 4 2−7.290 × 105 + 2.234 × 10−4j0 4 358−4.549 × 10−5 + 7.290 × 105j
0 4 358−7.290 × 105 − 4.228 × 10−4j
k1 k2 k3C3313k1 k2 k3C3323k1 k2 k3C3333
0 0 01.518 × 10−3 + 0j0 0 05.058 × 10−4 + 0j0 0 01.750 × 107 + 0j
0 2 1−2.916 × 106 + 1.574 × 10−3j0 2 11.420 × 10−3 + 2.916 × 106j0 2 01.166 × 107 + 1.264 × 10−1j
0 2 3592.916 × 106 − 2.529 × 10−3j0 2 3591.340 × 10−3 + 2.916 × 106j0 4 02.916 × 106 + 5.968 × 10−2j
0 4 1−1.458 × 106 − 1.493 × 10−3j0 4 1−9.979 × 10−4 + 1.458 × 106j
0 4 3591.458 × 106 + 2.746 × 10−4j0 4 359−1.115 × 10−3 + 1.458 × 106j

References

  1. Murat, F.; Tartar, L. Optimality conditions and homogenization. In Non-Linear Variational Problems; Pitman: London, UK, 1985. [Google Scholar]
  2. Cherkaev, A.V. Variational Methods for Structural Optimization; Springer: New York, NY, USA, 2000. [Google Scholar]
  3. Cherkaev, A.V.; Gibiansky, L.V. Extremal structures of multiphase heat conducting composites. Int. J. Solids Struct. 1996, 33, 2609–2623. [Google Scholar] [CrossRef]
  4. Lurie, K.A. A stable spatio-temporal G-closure and Gm-closure of a set of isotropic dielectrics with respect to one-dimensional wave propagation. Wave Motion 2004, 40, 95–110. [Google Scholar] [CrossRef]
  5. Adams, B.L.; Henrie, A.; Henrie, B.; Lyon, M.; Kalidindi, S.R.; Garmestani, H. Microstructure-sensitive design of a compliant beam. J. Mech. Phys. Solids 2001, 49, 1639–1663. [Google Scholar] [CrossRef]
  6. Kalidindi, S.R.; Houskamp, J.R.; Lyons, M.; Adams, B.L. Microstructure sensitive design of an orthotropic plate subjected to tensile load. Int. J. Plast. 2004, 20, 1561–1575. [Google Scholar] [CrossRef]
  7. Wu, X.; Proust, G.; Knezevic, M.; Kalidindi, S.R. Elastic-plastic property closures for hexagonal close-packed polycrystalline metals using first-order bounding theories. Acta Mater. 2007, 55, 2729–2737. [Google Scholar] [CrossRef]
  8. Knezevic, M.; Landry, N.W. Procedures for reducing large datasets of crystal orientations using generalized spherical harmonics. Mech. Mater. 2015, 88, 73–86. [Google Scholar] [CrossRef]
  9. Kalidindi, S.R.; Houskamp, J.R. Application of the spectral methods of microstructure design to continuous fiber reinforced composites. J. Compos. Mater. 2007, 41, 909–930. [Google Scholar] [CrossRef]
  10. Bunge, H.-J. Texture Analysis in Materials Science. Mathematical Methods; Cuvillier Verlag: Göttingen, Germany, 1993. [Google Scholar]
  11. Zuo, L.; Humbert, M.; Esling, C. Elastic properties of polycrystals in the voigt-reuss-hill approximation. J. Appl. Crystallogr. 1992, 25, 751–755. [Google Scholar] [CrossRef]
  12. Proust, G.; Kalidindi, S.R. Procedures for construction of anisotropic elastic-plastic property closures for face-centered cubic polycrystals using first-order bounding relations. J. Mech. Phys. Solids 2006, 54, 1744–1762. [Google Scholar] [CrossRef]
  13. Knezevic, M.; Kalidindi, S.R. Fast computation of first-order elastic-plastic closures for polycrystalline cubic-orthorhombic microstructures. Comput. Mater. Sci. 2007, 39, 643–648. [Google Scholar] [CrossRef]
  14. Duvvuru, H.K.; Knezevic, M.; Mishra, R.K.; Kalidindi, S.R. Application of microstructure sensitive design to fcc polycrystals. Mater. Sci. Forum 2007, 546, 675–680. [Google Scholar] [CrossRef]
  15. Knezevic, M.; Kalidindi, S.R.; Mishra, R.K. Delineation of first-order closures for plastic properties requiring explicit consideration of strain hardening and crystallographic texture evolution. Int. J. Plast. 2008, 24, 327–342. [Google Scholar] [CrossRef]
  16. Knezevic, M.; Kalidindi, S.R.; Fullwood, D. Computationally efficient database and spectral interpolation for fully plastic taylor-type crystal plasticity calculations of face-centered cubic polycrystals. Int. J. Plast. 2008, 24, 1264–1276. [Google Scholar] [CrossRef]
  17. Fuentes-Cobas, L.E.; Muñoz-Romero, A.; Montero-Cabrera, M.E.; Fuentes-Montero, L.; Fuentes-Montero, M.E. Predicting the coupling properties of axially-textured materials. Materials 2013, 6, 4967–4984. [Google Scholar] [CrossRef]
  18. Paul, B. Prediction of elastic constants of multiphase materials. Trans. Metall. Soc. AIME 1960, 218, 36–41. [Google Scholar]
  19. Hill, R. The elastic behavior of a crystalline aggregate. Proc. R. Soc. Lond. Sect. A 1952, 65, 349–354. [Google Scholar] [CrossRef]
  20. Hill, R. Elastic properties of reinforced solids: Some theoretical principles. J. Mech. Phys. Solids 1963, 11, 357–372. [Google Scholar] [CrossRef]
  21. Lyon, M.; Adams, B.L. Gradient-based non-linear microstructure design. J. Mech. Phys. Solids 2004, 52, 2569–2586. [Google Scholar] [CrossRef]
  22. Kim, I.Y.; de Weck, O.L. Adaptive weighted-sum method for bi-objective optimization: Pareto front generation. Struct. Multidiscip. Optim. 2005, 29, 149–158. [Google Scholar] [CrossRef]
  23. Fast, T.; Knezevic, M.; Kalidindi, S.R. Application of microstructure sensitive design to structural components produced from hexagonal polycrystalline metals. Comput. Mater. Sci. 2008, 43, 374–383. [Google Scholar] [CrossRef]
  24. Sintay, D.S.; Adams, B.L. Microstructure design for a rotating disk: With application to turbine engines. In Proceedings of the IDETC/CIE, 31st Design Automation Conference, Long Beach, CA, USA, 24–28 September 2005.
  25. Houskamp, J.R.; Proust, G.; Kalidindi, S.R. Integration of microstructure-sensitive design with finite element methods: Elastic-plastic case studies in fcc polycrystals. Int. J. Multiscale Comput. Eng. 2007, 5, 261–272. [Google Scholar]
  26. Saheli, G.; Garmestani, H.; Adams, B.L. Microstructure design of a two phase composite using two-point correlation functions. J. Comput. Aided Mater. Des. 2004, 11, 103–115. [Google Scholar] [CrossRef]
  27. Kalidindi, S.R.; Binci, M.; Fullwood, D.; Adams, B.L. Elastic properties closures using second-order homogenization theories: Case studies in composites of two isotropic constituents. Acta Mater. 2006, 54, 3117–3126. [Google Scholar] [CrossRef]
  28. Sigmund, O.; Torquato, S. Design of materials with extreme thermal expansion using a three-phase topology optimization method. J. Mech. Phys. Solids 1997, 45, 1037–1067. [Google Scholar] [CrossRef]
  29. Sigmund, O.; Torquato, S. Composites with extremal thermal expansion coefficients. Appl. Phys. Lett. 1996, 69, 3203–3205. [Google Scholar] [CrossRef]
  30. Torquato, S. Modeling of physical properties of composite materials. Int. J. Solids Struct. 2000, 37, 411–422. [Google Scholar] [CrossRef]
  31. Xu, B.; Arias, F.; Brittain, S.T.; Zhao, X.-M.; Grzybowski, B.; Torquato, S.; Whitesides, G.M. Making negative poisson’s ratio microstructures by soft lithography. Adv. Mater. 1999, 11, 1186–1189. [Google Scholar] [CrossRef]
  32. Mayeur, J.; Beyerlein, I.; Bronkhorst, C.; Mourad, H. The influence of grain interactions on the plastic stability of heterophase interfaces. Materials 2014, 7, 302–322. [Google Scholar] [CrossRef]
  33. Taylor, G.I. Plastic strain in metals. J. Inst. Metals 1938, 62, 307–324. [Google Scholar]
  34. Knezevic, M.; Drach, B.; Ardeljan, M.; Beyerlein, I.J. Three dimensional predictions of grain scale plasticity and grain boundaries using crystal plasticity finite element models. Comput. Methods Appl. Mech. Eng. 2014, 277, 239–259. [Google Scholar] [CrossRef]
  35. Zecevic, M.; Knezevic, M. A dislocation density based elasto-plastic self-consistent model for the prediction of cyclic deformation: Application to al6022-t4. Int. J. Plast. 2015, 72, 200–217. [Google Scholar] [CrossRef]
  36. Zecevic, M.; Knezevic, M.; Beyerlein, I.J.; Tomé, C.N. An elasto-plastic self-consistent model with hardening based on dislocation density, twinning and de-twinning: Application to strain path changes in hcp metals. Mater. Sci. Eng. A 2015, 638, 262–274. [Google Scholar] [CrossRef]
  37. Jahedi, M.; Ardeljan, M.; Beyerlein, I.J.; Paydar, M.H.; Knezevic, M. Enhancement of orientation gradients during simple shear deformation by application of simple compression. J. Appl. Phys. 2015, 117. [Google Scholar] [CrossRef]
  38. Knezevic, M.; Beyerlein, I.J.; Brown, D.W.; Sisneros, T.A.; Tomé, C.N. A polycrystal plasticity model for predicting mechanical response and texture evolution during strain-path changes: Application to beryllium. Int. J. Plast. 2013, 49, 185–198. [Google Scholar] [CrossRef]
  39. Knezevic, M.; Beyerlein, I.J.; Lovato, M.L.; Tomé, C.N.; Richards, A.W.; McCabe, R.J. A strain-rate and temperature dependent constitutive model for bcc metals incorporating non-schmid effects: Application to tantalum-tungsten alloys. Int. J. Plast. 2014, 62, 93–104. [Google Scholar] [CrossRef]
  40. Knezevic, M.; Beyerlein, I.J.; Nizolek, T.; Mara, N.A.; Pollock, T.M. Anomalous basal slip activity in zirconium under high-strain deformation. Mater. Res. Lett. 2013, 1, 133–140. [Google Scholar] [CrossRef]
  41. Knezevic, M.; Capolungo, L.; Tomé, C.N.; Lebensohn, R.A.; Alexander, D.J.; Mihaila, B.; McCabe, R.J. Anisotropic stress-strain response and microstructure evolution of textured a-uranium. Acta Mater. 2012, 60, 702–715. [Google Scholar] [CrossRef]
  42. Knezevic, M.; Carpenter, J.S.; Lovato, M.L.; McCabe, R.J. Deformation behavior of the cobalt-based superalloy haynes 25: Experimental characterization and crystal plasticity modeling. Acta Mater. 2014, 63, 162–168. [Google Scholar] [CrossRef]
  43. Knezevic, M.; McCabe, R.J.; Tomé, C.N.; Lebensohn, R.A.; Chen, S.R.; Cady, C.M.; Gray Iii, G.T.; Mihaila, B. Modeling mechanical response and texture evolution of α-uranium as a function of strain rate and temperature using polycrystal plasticity. Int. J. Plast. 2013, 43, 70–84. [Google Scholar] [CrossRef]
  44. Knezevic, M.; Zecevic, M.; Beyerlein, I.J.; Bingert, J.F.; McCabe, R.J. Strain rate and temperature effects on the selection of primary and secondary slip and twinning systems in hcp Zr. Acta Mater. 2015, 88, 55–73. [Google Scholar] [CrossRef]
  45. Lentz, M.; Klaus, M.; Beyerlein, I.J.; Zecevic, M.; Reimers, W.; Knezevic, M. In situ X-ray diffraction and crystal plasticity modeling of the deformation behavior of extruded Mg-Li-(Al) alloys: An uncommon tension-compression asymmetry. Acta Mater. 2015, 86, 254–268. [Google Scholar] [CrossRef]
  46. Lentz, M.; Klaus, M.; Wagner, M.; Fahrenson, C.; Beyerlein, I.J.; Zecevic, M.; Reimers, W.; Knezevic, M. Effect of age hardening on the deformation behavior of an Mg-Y-Nd alloy: In-situ X-ray diffraction and crystal plasticity modeling. Mater. Sci. Eng. A 2015, 628, 396–409. [Google Scholar] [CrossRef]
  47. Fromm, B.S.; Adams, B.L.; Ahmadi, S.; Knezevic, M. Grain size and orientation distributions: Application to yielding of a-titanium. Acta Mater. 2009, 57, 2339–2348. [Google Scholar] [CrossRef]
  48. Carpenter, J.S.; Nizolek, T.; McCabe, R.J.; Knezevic, M.; Zheng, S.J.; Eftink, B.P.; Scott, J.E.; Vogel, S.C.; Pollock, T.M.; Mara, N.A.; et al. Bulk texture evolution of nanolamellar Zr-Nb composites processed via accumulative roll bonding. Acta Mater. 2015, 92, 97–108. [Google Scholar] [CrossRef]
  49. Van Houtte, P. Fast calculation of average taylor factors and mandel spins for all possible strain modes. Int. J. Plast. 2001, 17, 807–818. [Google Scholar] [CrossRef]
  50. Shaffer, J.B.; Knezevic, M.; Kalidindi, S.R. Building texture evolution networks for deformation processing of polycrystalline fcc metals using spectral approaches: Applications to process design for targeted performance. Int. J. Plast. 2010, 26, 1183–1194. [Google Scholar] [CrossRef]
  51. Kalidindi, S.R.; Duvvuru, H.K.; Knezevic, M. Spectral calibration of crystal plasticity models. Acta Mater. 2006, 54, 1795–1804. [Google Scholar] [CrossRef]
  52. Kalidindi, S.R.; Knezevic, M.; Niezgoda, S.; Shaffer, J. Representation of the orientation distribution function and computation of first-order elastic properties closures using discrete fourier transforms. Acta Mater. 2009, 57, 3916–3923. [Google Scholar] [CrossRef]
  53. Knezevic, M.; Al-Harbi, H.F.; Kalidindi, S.R. Crystal plasticity simulations using discrete fourier transforms. Acta Mater. 2009, 57, 1777–1784. [Google Scholar] [CrossRef]
  54. Al-Harbi, H.F.; Knezevic, M.; Kalidindi, S.R. Spectral approaches for the fast computation of yield surfaces and first-order plastic property closures for polycrystalline materials with cubic-triclinic textures. CMC Comput. Mater. Contin. 2010, 15, 153–172. [Google Scholar] [CrossRef]
  55. Zecevic, M.; McCabe, R.J.; Knezevic, M. A new implementation of the spectral crystal plasticity framework in implicit finite elements. Mech. Mater. 2015, 84, 114–126. [Google Scholar] [CrossRef]
  56. Mihaila, B.; Knezevic, M.; Cardenas, A. Three orders of magnitude improved efficiency with high-performance spectral crystal plasticity on GPU platforms. Int. J. Numer. Methods Eng. 2014, 97, 785–798. [Google Scholar] [CrossRef]
  57. Knezevic, M.; Savage, D.J. A high-performance computational framework for fast crystal plasticity simulations. Comput. Mater. Sci. 2014, 83, 101–106. [Google Scholar] [CrossRef]
  58. Savage, D.J.; Knezevic, M. Computer implementations of iterative and non-iterative crystal plasticity solvers on high performance graphics hardware. Comput. Mech. 2015, in press. [Google Scholar] [CrossRef]
  59. Bunge, H. Texture analysis in Materials Science; Butterworths: London, UK, 1982. [Google Scholar]
  60. Carslaw, H.S. Introduction to the Theory of Fourier’s Series and Integrals, 2nd ed.; MacMillan and Co. Limited: London, UK, 1921. [Google Scholar]
  61. Press, W.H.; Teukolsky, S.A.; Vetterling, W.T.; Flannery, B.P. Numerical Recipes: The Art of Scientific Computing; Cambridge University Press: Oxford, UK, 1992. [Google Scholar]
  62. Briggs, W.L.; Henson, V.E. The DFT: An Owner’s Manual for the Discrete Fourier Transform; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 1995. [Google Scholar]
  63. Brigham, E.O. The Fast Fourier Transform and Its Applications; Prentice Hall: Englewood Cliffs, NJ, USA, 1988. [Google Scholar]
  64. Cooley, J.W.; Tukey, J.W. Algorithm for the machine computation of complex fourier series. Math. Comput. 1965, 19, 297–301. [Google Scholar] [CrossRef]
  65. Duhamel, P.; Vetterli, M. Fast fourier transforms: A tutorial review and a state of the art. Signal Process. 1990, 19, 259–299. [Google Scholar]
  66. Adams, B.L.; Gao, X.; Kalidindi, S.R. Finite approximations to the second-order properties closure in single phase polycrystals. Acta Mater. 2005, 53, 3563–3577. [Google Scholar] [CrossRef]
  67. Meyers, M.A.; Chawla, K.K. Mechanical Behavior of Materials; Prentice Hall: Upper Saddle River, NJ, USA, 1998. [Google Scholar]
  68. Zecevic, M.; McCabe, R.J.; Knezevic, M. Spectral database solutions to elasto-viscoplasticity within finite elements: Application to a cobalt-based fcc superalloy. Int. J. Plast. 2015, 70, 151–165. [Google Scholar] [CrossRef]
  69. Knezevic, M.; McCabe, R.J.; Lebensohn, R.A.; Tomé, C.N.; Liu, C.; Lovato, M.L.; Mihaila, B. Integration of self-consistent polycrystal plasticity with dislocation density based hardening laws within an implicit finite element framework: Application to low-symmetry metals. J. Mech. Phys. Solids 2013, 61, 2034–2046. [Google Scholar] [CrossRef]
  70. Knezevic, M.; Lebensohn, R.A.; Cazacu, O.; Revil-Baudard, B.; Proust, G.; Vogel, S.C.; Nixon, M.E. Modeling bending of α-titanium with embedded polycrystal plasticity in implicit finite elements. Mater. Sci. Eng. A 2013, 564, 116–126. [Google Scholar] [CrossRef]
  71. Knezevic, M.; Jahedi, M.; Korkolis, Y.P.; Beyerlein, I.J. Material-based design of the extrusion of bimetallic tubes. Comput. Mater. Sci. 2014, 95, 63–73. [Google Scholar] [CrossRef]
  72. Knezevic, M.; McCabe, R.J.; Lebensohn, R.A.; Tomé, C.N.; Mihaila, B. Finite Element Implementation of a Self-Consistent Polycrystal Plasticity Model: Application to α-Uranium, Proceedings; Materials Science and Technology Division, Los Alamos National Laboratory: Los Alamos, NM, USA, 2012; Volume 2, pp. 789–796. [Google Scholar]
  73. Jahedi, M.; Paydar, M.H.; Knezevic, M. Enhanced microstructural homogeneity in metal-matrix composites developed under high-pressure-double-torsion. Mater. Charact. 2015, 104, 92–100. [Google Scholar] [CrossRef]
  74. Jahedi, M.; Paydar, M.H.; Zheng, S.; Beyerlein, I.J.; Knezevic, M. Texture evolution and enhanced grain refinement under high-pressure-double-torsion. Mater. Sci. Eng. A 2014, 611, 29–36. [Google Scholar] [CrossRef]
  75. Jahedi, M.; Knezevic, M.; Paydar, M. High-pressure double torsion as a severe plastic deformation process: Experimental procedure and finite element modeling. J. Mater. Eng. Perform. 2015, 24, 1471–1482. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Landry, N.W.; Knezevic, M. Delineation of First-Order Elastic Property Closures for Hexagonal Metals Using Fast Fourier Transforms. Materials 2015, 8, 6326-6345. https://0-doi-org.brum.beds.ac.uk/10.3390/ma8095303

AMA Style

Landry NW, Knezevic M. Delineation of First-Order Elastic Property Closures for Hexagonal Metals Using Fast Fourier Transforms. Materials. 2015; 8(9):6326-6345. https://0-doi-org.brum.beds.ac.uk/10.3390/ma8095303

Chicago/Turabian Style

Landry, Nicholas W., and Marko Knezevic. 2015. "Delineation of First-Order Elastic Property Closures for Hexagonal Metals Using Fast Fourier Transforms" Materials 8, no. 9: 6326-6345. https://0-doi-org.brum.beds.ac.uk/10.3390/ma8095303

Article Metrics

Back to TopTop