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 ...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.展开更多
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 au...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.展开更多
文摘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.
基金Supported by the National Social Science Foundation of China (No. 09BTJ012)Scientific Research Fund ofHunan Provincial Education Department (No. 09c390)+1 种基金supported in part by a HKUSeed Funding Program for Basic Research (Project No. 2009-1115-9042)a grant from Hong Kong ResearchGrant Council-General Research Fund (Project No. HKU779210M)
文摘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.
