Elasticity, stability, and ideal strength of beta-SiC in plane-wave-based ab initio calculations

On the basis of the pseudopotential plane-wave method and the local-density-functional theory, this paper studies energetics, stress-strain relation, stability, and ideal strength of beta-SiC under various loading modes, where uniform uniaxial extension and tension and biaxial proportional extension are considered along directions [001] and [111]. The lattice constant, elastic constants, and moduli of equilibrium state are calculated and the results agree well with the experimental data. As the four SI-C bonds along directions [111], [(1) over bar 11], [11(1) over bar] and [111] are not the same under the loading along [111], internal relaxation and the corresponding internal displacements must be considered. We find that, at the beginning of loading, the effect of internal displacement through the shuffle and glide plane diminishes the difference among the four Si-C bonds lengths, but will increase the difference at the subsequent loading, which will result in a crack nucleated on the {111} shuffle plane and a subsequently cleavage fracture. Thus the corresponding theoretical strength is 50.8 GPa, which agrees well with the recent experiment value, 53.4 GPa. However, with the loading along [001], internal relaxation is not important for tetragonal symmetry. Elastic constants during the uniaxial tension along [001] are calculated. Based on the stability analysis with stiffness coefficients, we find that the spinodal and Born instabilities are triggered almost at the same strain, which agrees with the previous molecular-dynamics simulation. During biaxial proportional extension, stress and strength vary proportionally with the biaxial loading ratio at the same longitudinal strain.

A new type of boundary integral-equation for plane problems of elasticity including rotational forces

Based on the idea proposed by Hu [Scientia Sinica Series A XXX, 385-390 (1987)], a new type of boundary integral equation for plane problems of elasticity including rotational forces is derived and its boundary element formulation is presented. Numerical results for a rotating hollow disk are given to demonstrate the accuracy of the new type of boundary integral equation.

Plane wave approximation in linear elasticity

We consider the approximation of solutions of the time-harmonic linear elastic wave equation by linear combinations of plane waves. We prove algebraic orders of convergence both with respect to the dimension of the approximating space and to the diameter of the domain. The error is measured in Sobolev norms and the constants in the estimates explicitly depend on the problem wavenumber. The obtained estimates can be used in the h- and p-convergence analysis of wave-based finite element schemes.

Optimum distribution of material in sandwich plates loaded in their plane. Including plates with varying modulus of elasticity.

Mode of access: Internet.

Local plastic mechanisms in thin steel plates under in-plane compression

Thin-walled steel plates subjected to in-plane compression develop two types of local plastic mechanism, namely the roof-shaped mechanism and the so-called flip-disc mechanism, but the intriguing question of why two mechanisms should develop was not answered until recently. It was considered that the location of first yield point shifted from the centre of the plate to the midpoint of the longitudinal edge depending on the b/t ratio, imperfection level, and yield stress of steel, which then decided the type of mechanism. This paper has verified this hypothesis using analysis and laboratory experiments. An elastic analysis using Galerkin's method to solve Marguerre's equations was first used to determine the first yield point, based on which the local plastic mechanism/imperfection tolerance tables have been developed which give the type of mechanism as a function of b/t ratio, imperfection level and yield stress of steel. Laboratory experiments of thin-walled columns verified the imperfection tolerance tables and thus indirectly the hypothesis. Elastic and rigid-plastic curves were them used to predict the effect on the ultimate load due to the change of mechanism. A finite element analysis of selected cases also confirmed the results from simple analyses and experiments.

Beneficial defects : exploiting the intrinsic polishing-induced wafer roughness for the catalyst-free growth of Ge in-plane nanowires

We outline a metal-free fabrication route of in-plane Ge nanowires on Ge(001) substrates. By positively exploiting the polishing-induced defects of standard-quality commercial Ge(001) wafers, micrometer-length wires are grown by physical vapor deposition in ultra-high-vacuum environment. The shape of the wires can be tailored by the epitaxial strain induced by subsequent Si deposition, determining a progressive transformation of the wires in SiGe faceted quantum dots. This shape transition is described by finite element simulations of continuous elasticity and gives hints on the equilibrium shape of nanocrystals in the presence of tensile epitaxial strain.

A micropolar peridynamic theory in linear elasticity

A state-based micropolar peridynamic theory for linear elastic solids is proposed. The main motivation is to introduce additional micro-rotational degrees of freedom to each material point and thus naturally bring in the physically relevant material length scale parameters into peridynamics. Non-ordinary type modeling via constitutive correspondence is adopted here to define the micropolar peridynamic material. Along with a general three dimensional model, homogenized one dimensional Timoshenko type beam models for both the proposed micropolar and the standard non-polar peridynamic variants are derived. The efficacy of the proposed models in analyzing continua with length scale effects is established via numerical simulations of a few beam and plane-stress problems. (C) 2015 Elsevier Ltd. All rights reserved.

In-plane elastic responses of circular cell honeycombs and bundled circular tubes in a diamond array structure

Assemblages of circular tubes and circular honeycombs in close packed arrangement are presently both competing and complementing regular honeycomb structures (HCS). The intrinsic isotropy of bundled tubes/rings in hexagonal arrays restricts their use to applications with isotopic need. With the aim of extending the utility of tubes/rings assemblages to anisotropic needs, this paper explores the prospects of bundled tubes and circular honeycombs in a general diamond array structure (DAS) to cater these needs. To this end, effective transverse Young's moduli and Poisson's ratio for thick/thin DAS are obtained theoretically. Analysis frameworks including thin ring theory (TRT), curved beam theory (CBT) and elasticity formulations are tested and corroborated by FEA employing contact elements. Results indicate that TRT and CBT are reasonable for thin tubes and honeycombs. Nevertheless, TRT yields compact formulae to study the anisotropy ratio, moduli spectrum and sensitivity of the assemblage as a function of thicknesses and array structure. These formulae supplement designers as a guide to tailor the structures. On the other hand, elasticity formulation can estimate over a larger range including very thick tubes/rings. In addition, this formulation offers to estimate refined transverse strengths of assemblages. (C) 2015 Elsevier Ltd. All rights reserved.

In-plane effective shear modulus of generalized circular honeycomb structures and bundled tubes in a diamond array structure

Use of circular hexagonal honeycomb structures and tube assemblies in energy absorption systems has attracted a large number of literature on their characterization under crushing and impact loads. Notwithstanding these, effective shear moduli (G*) required for complete transverse elastic characterization and in analyses of hierarchical structures have received scant attention. In an attempt to fill this void, the present study undertakes to evaluate G* of a generalized circular honeycomb structures and tube assemblies in a diamond array structure (DAS) with no restriction on their thickness. These structures present a potential to realize a spectrum of moduli with minimal modifications, a point of relevance for manufactures and designers. To evaluate G* in this paper, models based on technical theories - thin ring theory and curved beam theory - and rigorous theory of elasticity are investigated and corroborated with FEA employing contact elements. Technical theories which give a good match for thin HCS offer compact expressions for moduli which can be harvested to study sensitivity of moduli on topology. On the other hand, elasticity model offers a very good match over a large range of thickness along with exact analysis of stresses by employing computationally efficient expressions. (C) 2015 Elsevier Ltd. All rights reserved.

Ab initio investigation of the elasticity and stability of aluminium

On the basis of the pseudopotential plane-wave (PP-PW) method in combination with the local density functional theory (LDFT), complete stress-strain curves for the uniaxial loading and uniaxial deformation along the [001] and [111] directions, and the biaxial proportional extension along [010] and [001] for aluminium are obtained. During the uniaxial loading, certain general behaviours of the energy versus the stretch and the load versus the stretch are confirmed; in each case, there exist three special unstressed structures: f.c.c., b.c.c., and f.c.t. for [001]; f.c.c., s.c., and b.c.c. for [111]. Using stability criteria, we find that all of these states are unstable, and always occur together with shear instability, except the natural f.c.c. structure. A Pain transformation from the stable f.c.c. structure to the stable b.c.c. configuration cannot be obtained by uniaxial compression along any equivalent [001] and [111] direction. The tensile strengths are similar for the two directions. For the higher energy barrier of the [111] direction, the compressive strength is greater than that for the [001] direction. With increase in the ratio of the biaxial proportional extension, the stress and tensile strength increase; however, the critical strain does not change significantly. Our results add to the existing ab initio database for use in fitting and testing interatomic potentials.

Some problems in nonlinear elasticity

Two separate problems are discussed: axisymmetric equilibrium configurations of a circular membrane under pressure and subject to thrust along its edge, and the buckling of a circular cylindrical shell.

An ordinary differential equation governing the circular membrane is imbedded in a family of n-dimensional nonlinear equations. Phase plane methods are used to examine the number of solutions corresponding to a parameter which generalizes the thrust, as well as other parameters determining the shape of the nonlinearity and the undeformed shape of the membrane. It is found that in any number of dimensions there exists a value of the generalized thrust for which a countable infinity of solutions exist if some of the remaining parameters are made sufficiently large. Criteria describing the number of solutions in other cases are also given.

Donnell-type equations are used to model a circular cylindrical shell. The static problem of bifurcation of buckled modes from Poisson expansion is analyzed using an iteration scheme and pertubation methods. Analysis shows that although buckling loads are usually simple eigenvalues, they may have arbitrarily large but finite multiplicity when the ratio of the shell's length and circumference is rational. A numerical study of the critical buckling load for simple eigenvalues indicates that the number of waves along the axis of the deformed shell is roughly proportional to the length of the shell, suggesting the possibility of a "characteristic length." Further numerical work indicates that initial post-buckling curves are typically steep, although the load may increase or decrease. It is shown that either a sheet of solutions or two distinct branches bifurcate from a double eigenvalue. Furthermore, a shell may be subject to a uniform torque, even though one is not prescribed at the ends of the shell, through the interaction of two modes with the same number of circumferential waves. Finally, multiple time scale techniques are used to study the dynamic buckling of a rectangular plate as well as a circular cylindrical shell; transition to a new steady state amplitude determined by the nonlinearity is shown. The importance of damping in determining equilibrium configurations independent of initial conditions is illustrated.

Elasticity, band-gap bowing, and polarization of AlxGa1-xN alloys

Elastic constants, the bulk modulus, Young's modulus, band-gap bowing coefficients, spontaneous and piezoelectric polarizations, and piezoelectric coefficients of hexagonal AlxGa1-xN ternary alloys are calculated using first-principles methods. The fully relaxed structures and the structures subjected to homogeneous biaxial and uniaxial tension are investigated. We show that the biaxial tension in the plane perpendicular to the c axis and the uniaxial tension along the c axis all reduce the bulk modulus, whereas they reduce and enhance Young's modulus, respectively. We find that the biaxial and uniaxial tension can enhance the bowing coefficients. We also find that the biaxial tension can enhance the total polarization, while the uniaxial tension will suppress the total polarization. (C) 2008 American Institute of Physics.

The effect of anisotropic elasticity on the yielding characteristics of overconsolidated natural clay

Quantitative application of elastoplastic theory to the yielding behaviour of natural soils has always been uncertain. Part of the reason is that the theory was developed for reconstituted materials with isotropic structure, in contrast to natural soils that are usually anisotropic. The approach considered in this study assumes that pre-yielding behaviour is governed by the theory of linear anisotropic elasticity and that yield loci in the mean effective stress ( p') – deviator stress (q) plane are aligned approximately along the coefficient of earth pressure (K0) line. The assumption of a rotated yield locus associated with anisotropic elastic behaviour within the state boundary surface indicates that the elastic wall within the state boundary surface is inclined. The form of the state boundary surface has been determined mathematically in terms of anisotropic elastic and Cam-Clay soil parameters. Stress path tests were conducted on samples of Belfast Upper Boulder Clay removed from a depth of 28 m below ground surface. Good agreement was found between predicted and measured yield loci. The study also examined the influence of subsequent isotropic compression on the yielding characteristics of the natural clay. The indications are that the anisotropy developed during deposition disappears when the sample is loaded to a stress level at least twice the stress generated during the original deposition process. The methods developed in the paper have also been applied to test results reported previously on Winnipeg clay, and good agreement was obtained.

A smooth evolutionary structural optimization procedure applied to plane stress problem

Topological optimization problems based on stress criteria are solved using two techniques in this paper. The first technique is the conventional Evolutionary Structural Optimization (ESO), which is known as hard kill, because the material is discretely removed; that is, the elements under low stress that are being inefficiently utilized have their constitutive matrix has suddenly reduced. The second technique, proposed in a previous paper, is a variant of the ESO procedure and is called Smooth ESO (SESO), which is based on the philosophy that if an element is not really necessary for the structure, its contribution to the structural stiffness will gradually diminish until it no longer influences the structure; its removal is thus performed smoothly. This procedure is known as "soft-kill"; that is, not all of the elements removed from the structure using the ESO criterion are discarded. Thus, the elements returned to the structure must provide a good conditioning system that will be resolved in the next iteration, and they are considered important to the optimization process. To evaluate elasticity problems numerically, finite element analysis is applied, but instead of using conventional quadrilateral finite elements, a plane-stress triangular finite element was implemented with high-order modes for solving complex geometric problems. A number of typical examples demonstrate that the proposed approach is effective for solving problems of bi-dimensional elasticity. (C) 2014 Elsevier Ltd. All rights reserved.

Some surface profiles of a strip after plane-strain indentation by rigid bodies with serrated surfaces

This paper is concerned with the surface profiles of a strip after rigid bodies with serrated (saw-teeth) surfaces indent the strip and are subsequently removed. Plane-strain conditions are assumed. This has application in roughness transfer of final metal forming process. The effects of the semi-angle of the teeth, the depth of indentation and the friction on the contact surface on the profile are considered.