Faster computation of the Karhunen-Loeve expansion using its domain independence property

The goal of this work is to reduce the cost of computing the coefficients in the Karhunen-Loeve (KL) expansion. The KL expansion serves as a useful and efficient tool for discretizing second-order stochastic processes with known covariance function. Its applications in engineering mechanics include discretizing random field models for elastic moduli, fluid properties, and structural response. The main computational cost of finding the coefficients of this expansion arises from numerically solving an integral eigenvalue problem with the covariance function as the integration kernel. Mathematically this is a homogeneous Fredholm equation of second type. One widely used method for solving this integral eigenvalue problem is to use finite element (FE) bases for discretizing the eigenfunctions, followed by a Galerkin projection. This method is computationally expensive. In the current work it is first shown that the shape of the physical domain in a random field does not affect the realizations of the field estimated using KL expansion, although the individual KL terms are affected. Based on this domain independence property, a numerical integration based scheme accompanied by a modification of the domain, is proposed. In addition to presenting mathematical arguments to establish the domain independence, numerical studies are also conducted to demonstrate and test the proposed method. Numerically it is demonstrated that compared to the Galerkin method the computational speed gain in the proposed method is of three to four orders of magnitude for a two dimensional example, and of one to two orders of magnitude for a three dimensional example, while retaining the same level of accuracy. It is also shown that for separable covariance kernels a further cost reduction of three to four orders of magnitude can be achieved. Both normal and lognormal fields are considered in the numerical studies. (c) 2014 Elsevier B.V. All rights reserved.

Stress-strain fields and the effectiveness shear properties for three-phase composites with imperfect interface

Based on the 'average stress in the matrix' concept of Mori and Tanaka (:Mori, T., Tanaka, K., 1973. Average stress in matrix and average elastic energy of materials with misfitting inclusion. Acta Metall. 21, 571-580) a micromechanical model is presented for the prediction of the elastic fields in coated inclusion composites with imperfect interfaces. The solutions of the effective elastic moduli for this kind of composite are also obtained. In two kinds of composites with coated particulates and fibers, respectively, the interface imperfections are takes to the assumption that the interface displacement discontinues are linearly related to interface tractions like a spring layer of vanishing thickness. The resulting effective shear modulus for each material and the stress fields in the composite are presented under a transverse shear loading situation.

Micro/Nanomechanical and Tribological Properties of thin Diamond-Like Carbon Coatings

The diamond-like carbon (DLC) films with different thicknesses on 9Crl8 bearing steels were prepared using vacuum magnetic-filtering arc plasma deposition. Vickers indentation. nanoin-dentation and nanoscratch tests were used to characterize the DLC films with a wide range of applied loads. Mechanical and tribological behaviors of these submicron films were investigated and interpreted. The hardnesses of 9Crl8 and DLC, determined by nanoindentation, are approximately 8GPa and 60GPa respectively; their elastic moduli are approximately 25OGPa and 600GPa respectively. The friction coefficients of 9Crl8, DLC. organic coating, determined by nanoscratch, are approximately 0. 35, 0. 20 and 0. 13 respectively. It is demonstrated that nanoindentation and nanoscratch tests can provide more information about the near-surface elastic-plastic deformation, friction and wear properties. The correlation of mechanical properties and scratch resistance of DLC films on 9Crl8 steels can provide an assessment for the load-carrying capacity and wear resistance

Effect of inclusions and holes on the stiffness and strength of honeycombs

A finite element study has been performed on the effects of holes and rigid inclusions on the elastic modulus and yield strength of regular honeycombs under biaxial loading. The focus is on honeycombs that have already been weakened by a small degree of geometrical imperfection, such as a random distribution of fractured cell walls, as these imperfect honeycombs resemble commercially available metallic foams. Hashin-Shtrikman lower and upper bounds and self-consistent estimates of elastic moduli are derived to provide reference solutions to the finite element calculations. It is found that the strength of an imperfect honeycomb is relatively insensitive to the presence of holes and inclusions, consistent with recent experimental observations on commercial aluminum alloy foams.

Embedded cell model of three-dimensional two-phase composites

An embedded cell model is presented to obtain the effective elastic moduli and the elastic-plastic stress-strain relations of three-dimensional two-phase particulate composites. Each cell consists of an ellipsoidal inclusion surrounded by a finite ellipsoidal matrix that embedded in an infinite matrix. When both matrix and particle are elastic, the effective elastic moduli are derived which is an exact analytic formula without any simplified approximation that can be expressed in an explicit form. Further, the elastic-plastic stress-strain relations are obtained for spherical cells and oblate spheroid cells, in which the matrix is elastic and the particle is elastic-plastic. In addition, the macroscopic elastic-plastic constitutive relation of particle reinforced composites (PRC) is investigated by a systematic approach [1] in which the matrix is elastic-plastic and the particle is elastic.

Novel Mechanical Behavior Of Zno Nanorods

A novel stress-strain relation with two stages of linear elastic deformation is observed in [0 0 0 1]-oriented ZnO nanorods under uniaxial tensile loading. This phenomenon results from a phase transformation from wurtzite (WZ, P6(3)mc space group) to a body-centered tetragonal structure with four-atom rings (denoted as BCT-4) belonging to the P4(2)/mnm space group. The analysis here focuses on the effects of nanorod size and temperature on the phase transformation and the associated mechanical behavior. It is found that as size is increased from 19.5 to 45.5 angstrom, the critical stress for nucleation of the transformation decreases by 25% from 21.90 to 16.50 GPa and the elastic moduli of the WZ- and BCT-4-structured nanorods decrease by 24% (from 299.49 to 227.51 GPa) and 38% (from 269.29 to 166.86 GPa), respectively. A significant temperature effect is also observed, with the critical stress for transformation initiation decreasing 87.8% from 17.89 to 2.19 GPa as temperature increases from 300 to 1500 K. (c) 2007 Elsevier B.V. All rights reserved.

A New Model for Sediment Transport

The problem of predicting sediment transportation by water waves is treated analytically with the rate of wave energy dissipation or wave damping. With resorting to the theory of shallow water waves and the basis of Yamamoto’s Coulomb-damped poroelastic model, the Boussinesq-type equation has been derived over a variation depth bed. For convenience Cnoidal wave is just discussed, The Cnoidal wave with complex wave length and wave velocity, which are as a function of wave frequency, water depth, permeability, Poisson’s ratio and complex elastic moduli of bed soil, is applied to analyse the rate of sediment transportation. Considering the sediment transportation depended on the shear stress near-bed or the horizontal velocity, the conclusion of Yamamoto’s experiment in clay bed has been extended to general situation. It could be figured out that the model should provide a method to avoid the undistinguishable factors during sediment transport processes and relate mass transport with the sediment peculiarities.

Numerical investigation of crack growth in concrete subjected to compression by the generalized beam lattice model

The beam lattice-type models, such as the Euler-Bernoulli (or Timoshenko) beam lattice and the generalized beam (GB) lattice, have been proved very effective in simulating failure processes in concrete and rock due to its simplicity and easy implementation. However, these existing lattice models only take into account tensile failures, so it may be not applicable to simulation of failure behaviors under compressive states. The main aim in this paper is to incorporate Mohr-Coulomb failure criterion, which is widely used in many kinds of materials, into the GB lattice procedure. The improved GB lattice procedure has the capability of modeling both element failures and contact/separation of cracked elements. The numerical examples show its effectiveness in simulating compressive failures. Furthermore, the influences of lateral confinement, friction angle, stiffness of loading platen, inclusion of aggregates on failure processes are respectively analyzed in detail.

Short-range-order effects on intrinsic plasticity of metallic glasses

Through a systematical analysis of the elastic moduli for 137 metallic glasses (MGs) and 56 polycrystalline metals, we use a simple model developed by Knuyt et al. [J. Phys. F: Met. Phys. 16 (1986) p.1989; Phil. Mag. B 64 (1991) p.299] based on a Gaussian distribution for the first-neighbor distance to reveal the short-range-order (SRO) structural conditions for plasticity of MGs. It is found that the SRO structure with dense atomic packing, large packing dispersion and a significant anharmonicity of atomic interaction within an MG is favorable for its global plasticity. Although these conditions seem paradoxical, their perfect matching is believed to be a key for designing large plastic bulk MGs not only in compression but also in tension.

Mechanics of capillary forming of aligned carbon nanotube assemblies.

Elastocapillary self-assembly is emerging as a versatile technique to manufacture three-dimensional (3D) microstructures and complex surface textures from arrangements of micro- and nanoscale filaments. Understanding the mechanics of capillary self-assembly is essential to engineering of properties such as shape-directed actuation, anisotropic wetting and adhesion, and mechanical energy transfer and dissipation. We study elastocapillary self-assembly (herein called "capillary forming") of carbon nanotube (CNT) microstructures, combining in situ optical imaging, micromechanical testing, and finite element modeling. By imaging, we identify sequential stages of liquid infiltration, evaporation, and solid shrinkage, whose kinetics relate to the size and shape of the CNT microstructure. We couple these observations with measurements of the orthotropic elastic moduli of CNT forests to understand how the dynamic of shrinkage of the vapor-liquid interface is coupled to the compression of the forest. We compare the kinetics of shrinkage to the rate of evporation from liquid droplets having the same size and geometry. Moreover, we show that the amount of shrinkage during evaporation is governed by the ability of the CNTs to slip against one another, which can be manipulated by the deposition of thin conformal coatings on the CNTs by atomic layer deposition (ALD). This insight is confirmed by finite element modeling of pairs of CNTs as corrugated beams in contact and highlights the coupled role of elasticity and friction in shrinkage and stability of nanoporous solids. Overall, this study shows that nanoscale porosity can be tailored via the filament density and adhesion at contact points, which is important to the development of lightweight multifunctional materials.

Structural, mechanical, thermodynamic, and electronic properties of thorium hydrides from first-principles

We perform first-principles calculations of the structural, electronic, mechanical, and thermodynamic properties of thorium hydrides (ThH2 and Th4H15) based on the density functional theory with generalized gradient approximation. The equilibrium geometries, the total and partial densities of states, charge density, elastic constants, elastic moduli, Poisson's ratio, and phonon dispersion curves for these materials are systematically investigated and analyzed in comparison with experiments and previous calculations. These results show that our calculated equilibrium structural parameters are well consistent with experiments. The Th-H bonds in all thorium hydrides exhibit weak covalent character, but the ionic properties for ThH2 and Th4H15 are different due to their different hydrogen concentration. It is found that while in ThH2 about 1.5 electrons transfer from each Th atom to H, in Th4H15 the charge transfer from each Th atom is around 2.1 electrons. Our calculated phonon spectrum for the stable body-centered tetragonal phase of ThH2 accords well with experiments. In addition we show that ThH2 in the fluorite phase is mechanically and dynamically unstable.

岩石中地震波速度和衰减的应力变化响应研究

Stress change is one of key factors in seismic nucleating and triggering; therefore for understanding and forecasting earthquakes, it is necessary to research on stress status and its changes in rocks. Propagating in underground structures, wave velocity and attenuation contain information on stress changes of the Earth’s interior. For a better understanding of relationship between seismic data and stress changes, modeling and ultrasonic test supply significant references. In this article, acoustoelastic theory is introduced to explain nonlinear elastic characteristics of rocks. Based on the acoustoelastic theory, a solid-fluid coupled model is given to calculate velocity under different stress for porous and liquid fulfilled rocks. Except for the stress-velocity relationship, effects of pore pressure induced stress changes on ultrasonic coda attenuation are also studied. Intrinsic attenuation quality factors are calculated for a comparison purpose. Finally, the relationship between elastic constants and stress changes is thoroughly investigated, a mixture model from two phases of Hooke media is introduced to explain the differences between dynamic and static moduli, a relation among wave length, wave velocities and elastic moduli considering dimension of microstructure, dimension and state of surface between phases is presented. The most important aspect of this work is exploring and establishing relationships between the seismic properties of rocks and changes of their stress conditions, which will have its application in earthquake forecast and seismic hazard.

边坡动力响应分析及应用研究

As a marginal subject, dynamic responses of slopes is not only an important problem of engineering geology (Geotechnical problem), but also of other subjects such as seismology, geophysics, seismic engineering and engineering seismic and so on. Owning to the gulf between different subjects, it is arduous to study dynamic responses of slopes and the study is far from ripeness. Studying on the dynamic responses of slopes is very important in theories as well as practices. Supported by hundreds of bibliographies, this paper systemically details the development process of this subject, introduces main means to analyze this subject, and then gives brief remarks to each means respectively. Engineering geology qualitative analysis is the base of slopes dynamic responses study. Because of complexity of geological conditions, engineering geology qualitative analysis is very important in slopes stability study, especially to rock slopes with complex engineering geology conditions. Based on research fruits of forerunners, this paper summarizes factors influencing slopes dynamic stability into five aspects as geology background, stratums, rock mass structure, and topography as well as hydrogeology condition. Based on rock mass structure controlling theory, engineering geology model of the slope is grouped into two typical classes, one is model with obvious controlling discontinuities, which includes horizontal bedded slope, bedding slope, anti-dip slope, slide as well as slope with base rock and weathered crust; the other is model without obvious controlling discontinuities, which includes homogeneous soil slope, joint rock mass slope. Study on slope failure mechanism under dynamic force, the paper concludes that there are two effects will appear in slope during strong earthquake, one is earthquake inertia force, the other is ultra pore pressure buildup. The two effects lead to failure of the slope. To different types of slope failure, the intensity of two effects acting on the slope is different too. To plastic flow failure, pore pressure buildup is dominant; to falling rock failure and toppling failure, earthquake inertia force is dominant in general. This paper briefly introduces the principle of Lagrangian element method. Through a lot of numerical simulations with FLAC3D, the paper comprehensively studies dynamic responses of slopes, and finds that: if the slope is low, displacement, velocity and acceleration are linear enlarging with elevation increasing in vertical direction; if the slope is high enough, displacement, velocity and acceleration are not linear with elevation any more, on the other hand, they fluctuate with certain rhythm. At the same time, the rhythm appears in the horizontal direction in the certain area near surface of the slope. The distribution form of isoline of displacement, velocity and acceleration in the section of the slope is remarkably affected by the slope angle. In the certain area near the slope surface, isoline of displacement, velocity and acceleration is parallel to the surface of the slope, in the mean time, the strike direction of the extreraum area is parallel to the surface of the slope too. Beyond this area, the isoline direction and the strike direction of the extremum area turn to horizontal with invariable distance. But the rhythm appearing or not has nothing to with the slope angle. The paper defines the high slope effect and the low slope effect of slopes dynamic responses, discusses the threshold height H^t of the dynamic high slope effect, and finds that AW is proportional to square root of the dynamic elastic moduli El P , at the same time, it is proportional to period Tof the dynamic input. Thus, the discriminant of H^t is achieved. The discriminant can tell us that to a slope, if its height is larger than one fifth of the wavelength, its response regular will be the dynamic high slope effect; on the other hand, its response regular will be the dynamic low slope effect. Based on these, the discriminant of different slopes taking on same response under the same dynamic input is put forward in this paper. At the same time, the paper studies distribution law of the rhythm extremum point of displacement, velocity and acceleration, and finds that there exists relationship of N = int among the slope height H, the number of the rhythm extremum VHlhro) point N and ffthre- Furthermore, the paper points out that if N^l, the response of the slope will be dynamic high slope effect; \fN

Stress correlations in granular materials: An entropic formulation

We study the response of dry granular materials to external stress using experiment, simulation, and theory. We derive a Ginzburg-Landau functional that enforces mechanical stability and positivity of contact forces. In this framework, the elastic moduli depend only on the applied stress. A combination of this feature and the positivity constraint leads to stress correlations whose shape and magnitude are extremely sensitive to the nature of the applied stress. The predictions from the theory describe the stress correlations for both simulations and experiments semiquantitatively. © 2009 The American Physical Society.

Effect of the crown design and interface lute parameters on the stress-state of a machined crown-tooth system: A finite element analysis

The effect of preparation design and the physical properties of the interface lute on the restored machined ceramic crown-tooth complex are poorly understood. The aim of this work was to determine, by means of three-dimensional finite element analysis (3D FEA) the effect of the tooth preparation design and the elastic modulus of the cement on the stress state of the cemented machined ceramic crown-tooth complex. The three-dimensional structure of human premolar teeth, restored with adhesively cemented machined ceramic crowns, was digitized with a micro-CT scanner. An accurate, high resolution, digital replica model of a restored tooth was created. Two preparation designs, with different occlusal morphologies, were modeled with cements of 3 different elastic moduli. Interactive medical image processing software (mimics and professional CAD modeling software) was used to create sophisticated digital models that included the supporting structures; periodontal ligament and alveolar bone. The generated models were imported into an FEA software program (hypermesh version 10.0, Altair Engineering Inc.) with all degrees of freedom constrained at the outer surface of the supporting cortical bone of the crown-tooth complex. Five different elastic moduli values were given to the adhesive cement interface 1.8 GPa, 4 GPa, 8 GPa, 18.3 GPa and 40 GPa; the four lower values are representative of currently used cementing lutes and 40 GPa is set as an extreme high value. The stress distribution under simulated applied loads was determined. The preparation design demonstrated an effect on the stress state of the restored tooth system. The cement elastic modulus affected the stress state in the cement and dentin structures but not in the crown, the pulp, the periodontal ligament or the cancellous and cortical bone. The results of this study suggest that both the choice of the preparation design and the cement elastic modulus can affect the stress state within the restored crown-tooth complex.