循环荷载作用下饱和黏性土的弹塑性双面模型

基于作者新近提出的一种混合塑性硬化准则和临界状态土力学理论建立了一个三维形式表述的适于饱和黏性土的弹塑性双面模型,此硬化准则是著名的M asing方法在三维应力空间中的推广。该模型的二维形式已采用饱和黏性土样室内动态三轴不排水剪切试验的结果进行了验证并已在另文发表。本文所建的三维弹塑性动本构模型与一个半经验的孔隙水压力增长模式相结合,通过有限元程序ABAQUS*对一个循环荷载作用下的条形基础下部地基土层内的位移、应力分布和超孔隙水压力随时间的变化进行了分析,计算结果表明模型预测的地基土的力学响应是符合实际的。     

基于弹塑性有限单元法的齿轮轴数值模拟

基于弹塑性有限单元法的一般原理,建立了齿轮轴热处理的热力耦合分析模型,对12Cr齿轮轴热处理工艺中温度、应力以及相变的分布情况进行了模拟,研究了加热规范的成形性,并与实际结果比较证实本模型的有效性,为有限单元法在热处理工艺模拟中的应用奠定了一定基础。     

Elastoplastic Finite Element Analysis and Strength Evaluation of Adhesive Butt Joints of Similar and Dissimilar Hollow Shafts Subjected to External Bending Moments

Stress distributions and deformation of adhesive butt joints are analyzed by an elastoplastic finite element method when the joints of similar and dissimilar shafts are subjected to external bending moments. The effects of the ratio of Young's modulus for the adherends to that for an adhesive and the effects of the adhesive thickness on the interface stress distribution are investigated. Joint strength is predicted by using the elastoplastic interface stress distributions. It is found that the singular stress at the edge of the interfaces increases with an increase of the ratio of Young's modulus. Measurement of strains in joints and experiments on joint strength were conducted. The numerical results are in fairly good agreement with the experimental results. It is observed that the joint strength for dissimilar shafts are smaller than those for similar shafts. A fracture of dissimilar adhesive up-bonded shafts occurred from the interface of the adherend with smaller Young's modulus. It is seen that joint strength increases as the adhesive thickness increases.     

Nonclassical constitutive model involving void evolution of casting magnesium alloy

The void evolution equation and the elastoplastic constitutive model of casting magnesium alloy were investigated. The void evolution equation consists of the void growth and the void nucleation equations. The void growth equation was obtained based on the continuous supposition of the material matrix, and the void nucleation equation was derived by assuming that the void nucleation follows a normal distribution. A softening function related to the void evolution was given. After the softening function was embedded to a nonclassical elastoplastic constitutive equation, a constitutive model involving void evolution was obtained. The numerical algorithm and the finite element procedure related to the constitutive model were developed and applied to the analysis of the distributions of the stress and the porosity of the notched cylindrical specimens of casting magnesium alloy ZL305. The computed results show satisfactory agreement with the experimental data.     

Efficient numerical methods for the large‐scale,parallel solution of elastoplastic contact problems

Quasi‐static elastoplastic contact problems are ubiquitous in many industrial processes and other contexts, and their numerical simulation is consequently of great interest in accurately describing and optimizing production processes. The key component in these simulations is the solution of a single load step of a time iteration. From a mathematical perspective, the problems to be solved in each time step are characterized by the difficulties of variational inequalities for both the plastic behavior and the contact problem. Computationally, they also often lead to very large problems. In this paper, we present and evaluate a complete set of methods that are (1) designed to work well together and (2) allow for the efficient solution of such problems. In particular, we use adaptive finite element meshes with linear and quadratic elements, a Newton linearization of the plasticity, active set methods for the contact problem, and multigrid‐preconditioned linear solvers. Through a sequence of numerical experiments, we show the performance of these methods. This includes highly accurate solutions of a three‐dimensional benchmark problem and scaling our methods in parallel to 1024 cores and more than a billion unknowns. Copyright © 2015 John Wiley & Sons, Ltd.     

基于Cosserat连续体理论的疲劳寿命概率可靠性分析

基于Cosserat连续体理论建立考虑应变梯度效应的粉末高温合金弹塑性模型.该模型考虑微尺度长度效应,并能消除缺陷局部化问题求解时的网格依赖性.基于参变量变分原理提出一种高效的求解方法,该方法将原非线性问题转化为互补问题求解,可大大提高求解效率和收敛性.针对一种假想涡轮盘模型,将Cosserat尺度参数和Basquin斜率作为随机变量,分别采用蒙特卡洛法与二阶矩法进行疲劳寿命概率可靠度计算.结果表明二阶矩法在保证精度的情况下具有较高的效率.     

连续大跨输煤栈桥在不同地震强度下的时程分析

大跨度输煤栈桥跨度大、受荷复杂,为全面了解此类栈桥结构的抗震性能和寻求合理的设计方法,利用ABAQUS对某连续大跨输煤栈桥进行了不同地震强度下的弹性及弹塑性时程分析。通过对比不同地震动强度下结构的整体反应及构件损伤特征,得出了该栈桥在给定地震波下出现的破坏模式、塑性发展特点及抗震薄弱部位,对其抗震性能做出了结论性评价。     

A four‐node corotational quadrilateral elastoplastic shell element using vectorial rotational variables

A four‐node corotational quadrilateral elastoplastic shell element is presented. The local coordinate system of the element is defined by the two bisectors of the diagonal vectors generated from the four corner nodes and their cross product. This local coordinate system rotates rigidly with the element but does not deform with the element. As a result, the element rigid‐body rotations are excluded in calculating the local nodal variables from the global nodal variables. The two smallest components of each nodal orientation vector are defined as rotational variables, leading to the desired additive property for all nodal variables in a nonlinear incremental solution procedure. Different from other existing corotational finite‐element formulations, the resulting element tangent stiffness matrix is symmetric owing to the commutativity of the local nodal variables in calculating the second derivative of strains with respect to these variables. For elastoplastic analyses, the Maxwell–Huber–Hencky–von Mises yield criterion is employed, together with the backward‐Euler return‐mapping method, for the evaluation of the elastoplastic stress state; the consistent tangent modulus matrix is derived. To eliminate locking problems, we use the assumed strain method. Several elastic patch tests and elastoplastic plate/shell problems undergoing large deformation are solved to demonstrate the computational efficiency and accuracy of the proposed formulation. Copyright © 2013 John Wiley & Sons, Ltd.     

高效实用的板带变形求解方法

针对轧件弹塑性变形与轧辊弹性变形计算分离的弊端,用轧制力分布系数将这2个模型结合起来,并解决了此2个模型之间的耦合问题.理论计算结果在1700mm热连轧机组样品进行验证,理论和实践两者符合较好.