Quantitative law of diffusion induced fracture

Through dimension analysis, an almost analytical model for the maximum diffusion induced stress（DIS）and critical temperature（or concentration） difference at which cracks begin to initiate in the diffusion process is developed. It interestingly predicts that the spacing of diffusioninduced cracks is constant, independent of the thickness of specimen and the temperature difference. These conclusions are validated by our thermal shock experiments on alumina plates. Furthermore, the proposed model can interpret observed hierarchical crack patterns for high temperature jump cases, and a three-stage relation between the residual strength and the temperature difference. The prediction for crack spacing can guide the biomimetic thermal-shockfailure proof design, in which the hard platelets smaller than the predicted diffusion induced by constant crack-spacing are embedded in a soft matrix, and, therefore, no fracture will happen. This may guide the design of the thermal protection system and the lithium ion battery. Finally we present the maximum normalized DISes for various geometry and boundary conditions by single-variable curves for the stressindependent diffusion process and two-variable contour plots for the stress-dependent diffusion process, which can provideengineers and materialists a simple and easy way to quickly evaluate the reliability of related materials and devices.

Efficient Algorithms for Generating Truncated Multivariate Normal Distributions

Sampling from a truncated multivariate normal distribution （TMVND） constitutes the core computational module in fitting many statistical and econometric models. We propose two efficient methods, an iterative data augmentation （DA） algorithm and a non-iterative inverse Bayes formulae （IBF） sampler, to simulate TMVND and generalize them to multivariate normal distributions with linear inequality constraints. By creating a Bayesian incomplete-data structure, the posterior step of the DA Mgorithm directly generates random vector draws as opposed to single element draws, resulting obvious computational advantage and easy coding with common statistical software packages such as S-PLUS, MATLAB and GAUSS. Furthermore, the DA provides a ready structure for implementing a fast EM algorithm to identify the mode of TMVND, which has many potential applications in statistical inference of constrained parameter problems. In addition, utilizing this mode as an intermediate result, the IBF sampling provides a novel alternative to Gibbs sampling and elimi- nares problems with convergence and possible slow convergence due to the high correlation between components of a TMVND. The DA algorithm is applied to a linear regression model with constrained parameters and is illustrated with a published data set. Numerical comparisons show that the proposed DA algorithm and IBF sampler are more efficient than the Gibbs sampler and the accept-reject algorithm.

ELASTIC CALCULATION OF STRESSES IN RINGS USING AIRY STRESS FUNCTION

Stress calculation formulae for a ring have been obtained by using Airy stress function of the plane strain field with the decomposition of the solutions for normal stresses of Airy biharmonic equation into two parts when it is loaded under two opposite inside forces along a diameter. One part should fulfill a constraint condition about normal stress distribution along the circumference at an energy valley to do the minimum work. Other part is a stress residue constant. In order to verify these formulae and the computed results, the computed contour lines of equi-maximal shear stresses were plotted and quite compared with that of photo-elasticity test results. This constraint condition about normal stress distribution along circumference is confirmed by using Greens’ theorem. An additional compression exists along the circumference of the loaded ring, explaining the divorcement and displacement of singularity points at inner and outer boundaries.