Yield criterion and elasto-plastic damage constitutive model for frozen sandy soil

A series of triaxial compression tests was carried out on a frozen sandy soil under confining pressures of 0–18 MPa at −6 °C. The experimental results indicate that, the strength of frozen sandy soil increases versus the increase in the confining pressures when σ₃ ? 3 MPa, but decreases when σ₃ > 3 MPa. This phenomenon is called the strengthening and weakening effects of confining pressures. A yield function, considering both effects, is proposed using the experimental method according to Drucker’s postulate, and the mathematical expression of the hardening parameter, which can describe the softening and hardening phenomenon, is provided. An elasto-plastic constitutive model for frozen sandy soil is developed. Based on the continuum damage theory, the cross anisotropic damage variables are deduced and their change regularities are investigated. Then the elasto-plastic damage constitutive model is proposed by introducing damage variables into elasto-plastic constitutive model. The validity of the model is verified by comparing its modeling results with experimental results obtained from triaxial tests. It is found that, this model can predict the deformation regularity of frozen soil exactly. It can simulate the stress–strain process under high confining pressures when pressure melting phenomena appear especially well.     

Stress space and strain space formulation of the elasto-plastic constitutive relations for singular yield surface

This paper gives the stress space and the strain space formulations of the elastoplastic constitutive relations at a singular point on a yield surface, discusses the parallelism of the two space formulations and points out that the strain space formulation has a wider range of applicability.     

Evolution of subsequent yield surfaces and elastic constants with finite plastic deformation. Part-I: A very low work hardening aluminum alloy (Al6061-T6511)

In the present study, the initial and subsequent yield surfaces in Al 6061-T6511, based on 10 με deviation from linearity definition of yield, are presented. The subsequent yield surfaces are determined during tension, free end torsion, and combined tension–torsion proportional loading paths after reaching different levels of strains. The yield surfaces are also obtained after linear, bi-linear and non-linear unloading paths after finite plastic deformation. The initial yield surface is very close to the von-Mises yield surface and the subsequent yield surfaces undergo translation and distortion. In the case of this low work hardening material, the size of the yield surfaces is smaller and negative cross-effect is observed with finite plastic deformation. The subsequent yield have a usual “nose” in the loading direction and flattened shape in the reverse loading direction; the observed nose is more dominant in the case of tension and combined tension–torsion loading than in torsional loading. The size of the yield surfaces after unloading is smaller than the initial yield surface but larger than subsequent yield surfaces obtained during prior loading, show much smaller cross-effect, and the shape of these yield surfaces depends strongly on the loading and unloading paths. Elastic constants (Young’s and shear moduli) are also measured within each subsequent yield surfaces. Evolution of these constants with finite deformation is also presented. The decrease of the two moduli is found to be much smaller than reported earlier in tension by Cleveland and Ghosh [Cleveland, R.M., Ghosh, A.K., 2002. Inelastic effects on springback in metals. Int. J. Plast. 18, 769–785]. Part-II and III [(Khan et al., 2009a) and (Khan et al., 2009b)] of the papers will include experimental results on annealed 1100 Al (a very high work hardening material) and on both Al alloys (Al6061-T6511 and annealed 1100 Al) in tension- tension stress space, respectively. The results for both cases are quite different than the ones that are presented in this paper.     

Evolution of subsequent yield surfaces and elastic constants with finite plastic deformation. Part II: A very high work hardening aluminum alloy (annealed 1100 Al)

Results are presented on the evolution of subsequent yield surfaces with finite deformation in a very high work hardening annealed 1100 aluminum alloy. In Part I [Khan, A.S., Kazmi, R., Stoughton, T., Pandey, A., 2009a. Evolution of subsequent yield surfaces and elastic constants with finite plastic deformation. Part 1: a very low work hardening aluminum alloy (Al-6061–T6511) 25, 1611–1625.] of this paper, similar results are presented for a very low work hardening aluminum alloy. Those results were very different from the present ones, and all the results were for proportional loading paths. The subsequent yield surfaces are determined in tension, free end torsion and combined tension–torsion proportional and non-proportional loading paths, using 10 με deviation from linearity definition of yield. Yield surfaces are also determined after linear, bi-linear, and non-linear unloading paths after finite deformation under tension, free end torsion, and combined tension–torsion loading. The initial yield surface is closer to the von-Mises surface and the subsequent yield surfaces show distortion, expansion, positive cross-effect, and “nose” in the loading direction. Additionally, the subsequent yield surfaces after non-proportional loading paths show shrinkage and compounded distortion. The yield surfaces after unloading depict strong anisotropy, positive cross-effect and exhibits different proportion of distortion in each loading conditions. The Young’s and shear modulus decrease with plastic deformation and this decrease is much less than those reported in the published literature.     

一种岩石损伤本构模型在地下强爆炸中的应用

为了描述地下爆炸波的传播,建立了一种考虑塑性硬化、剪胀和损伤软化效应的岩石本构模型,并应用于地下强爆炸自由场的数值计算,获得的速度和位移波形以及峰值速度衰减曲线等与国外地下试验测量数据和同类计算相比十分接近,从而验证了本构模型的有效性。     

基于柔度法与截面轴力弯矩屈服面本构的梁柱单元

基于截面层次的梁柱单元能兼顾计算精度与效率,但现有截面上的本构模型不能准确反映轴力存在与变化对弯矩的影响,且柔度法在构建梁柱单元中未发挥优势,为此提出了基于截面轴力与弯矩屈服面本构和柔度法的梁柱单元模型。首先,根据截面组合思想构建截面上轴力与弯矩的屈服面;然后,依托柔度法和塑性理论进行单元状态和截面状态确定;最后,利用提出的单元进行平面钢框架的动力弹塑性分析。结果表明,所选用的基于截面轴力弯矩屈服面的本构能很好考虑轴力存在及变化对弯矩的影响;相对于刚度法,使用一个由柔度法构建的单元就足够描述一个梁-柱构件,其计算精度接近基于柔度法的纤维单元模型。