含凹坑缺陷圆柱壳的数值极限分析

使用文［１］提出的三维结构塑性极限分析的一般计算方法，我们对含凹坑缺陷的圆柱壳进行了数值极限分析．对凹坑和筒体各种组合的几何参数，本文给出了筒壳极限压力的上限．计算结果与现有的理论、实验和数值解进行了比较．本文调查和评估了各种形状和尺寸的凹坑对筒壳极限承载能力的影响规律，研究了对应于不同凹坑尺寸的筒壳两种典型的破坏模式．根据以上数值结果，本文采用几何参数Ｇ来反映凹坑各参数对筒壳极限压力的综合影响，并给出了估计带凹坑筒体极限压力的拟合公式．本文结果对含凹坑缺陷压力容器的安全评估具有重要参考价值     

多孔材料塑性极限载荷及其破坏模式分析

运用塑性力学中的机动极限分析理论，研究韧性基体多孔材料的塑性极限承载能力和破坏模式。以多孔材料的细观结构为研究对象，将细观力学中的均匀化理论引入到塑性极限分析中，并结合有限元技术，建立细观结构极限载荷的一般计算格式，并提出相应的求解算法。数值算例表明：细观孔洞对材料的宏观强度影响明显；在单向拉伸作用下，孔洞呈现膨胀扩大规律；多孔材料破坏源于基体塑性区的贯通。     

含缺陷结构的塑性极限分析

结合极限分析中的数学规划理论和有限元技术,提出了三维含结构极限分析的数学规划方法,并采用罚函数法引入塑性不可压条件,对于考虑多组独立变化载荷作用的情况,提出了加载路径射线辐射求解方案,并基于这种射线射状的加载路径,推导了多组载荷联合作用下结构塑性极限上限分析的数学规划格式,编制了相应的有限元程序,文中的数值了该方法的正确性与有效性。     

平面刚架弹塑性大位移分析的多刚体离散元法

本文基于多刚体－弹簧系统模型，给出了求解平面刚架结构弹塑性、大位移极限承载力分析的多刚体离散元法。文中首先推导了多刚体离散元法在总体坐标下的切线刚度阵，建立多刚体离散元法的增量平衡方程；而后推导了多刚体离散元的弹塑性弹簧系数矩阵，建立了多刚体离散元内力屈服面塑性铰法的增量求解格式，成功地进行了平面钢框架的弹塑性、大位移极限承载力分析。计算结果与其他数值方法或实验结果吻合良好，显示了多刚体离散元方法进行结构极限承载力分析这一复杂问题的优越性     

土质边坡极限荷载上限计算方法研究

以边坡极限分析的上限解析式为标准,在无重力情况下利用不同屈服准则的数值方法计算了相同工况下边坡的极限荷载。针对一个实际边坡工况使用解析和FLAC方法计算了边坡的极限荷载上限,并且基于强度折减法使用FLAC确定了各种荷载下边坡的安全系数。通过与解析方法计算结果的比较发现,基于Mohr-Coulomb准则的FLAC软件计算出的边坡极限荷载更加接近解析解。不同边坡顶角或不同内摩擦角的计算结果显示FLAC计算出的边坡极限荷均稍大于解析解,并且内摩擦角越大差值越大,但两者的差别最大不超过5%。     

极限分析的无搜索数学规划算法

本文研究理想刚塑性介质极限载荷因子的计算方法。根据极限分权理论的上限定理,建立了计算极限载荷因子的一般数学规划有限元格式。针对这种格式的特点,提出了一个求解极限载荷因子的无搜索迭代算法。这个算法中采用逐步识别刚性、塑性分区,不断修正目标函数的方案,克服了目标函数非光滑所导致的困难。本文提出的算法建立于位移模式有限元基础上,有较广的适用范围,且具有计算效率高,稳定性好,格式简单易于程序实现等优点。     

轴对称结构极限下限分析的自然单元法

作为一种基于自然邻近插值的新型无网格法,自然单元法克服了大多数无网格法难以施加本质边界条件的困难.将自然单元法与减缩基技术相结合,建立了一种轴对称结构极限下限分析的数值格式和求解算法.通过不断修正自平衡应力场,轴对称结构极限下限分析可转化为一系列的非线性数学规划子问题,并由复合形法求解.在每个非线性规划子问题中,自平衡...     

半潜式钻井平台管道钢构支架极限强度研究

钢构支架是半潜式钻井平台管道支吊架的主要类型之一,钢构支架的极限强度是管道系统正常工作的重要保障。研究结构极限强度的方法有理论分析、有限元数值仿真和实验分析。在分析极限理论的基础上由静力法计算四种钢构支架试件的极限载荷;利用MSC软件对试件极限载荷进行数值仿真分析,并对有限元模型化技术进行讨论;对试件进行实验极限载荷测试,对比分析了三种方法测得的极限载荷。结果表明,三种分析方法计算得到结构的极限载荷基本一致,对于结构形式较为简单的结构通过理论分析可以得到简化解析解析解;数值仿真分析中采用合理的有限元模型化技术(结构有限元模型、边界、约束等)可得到精度较高的计算结果。     

弹性连杆机构广义刚度可靠性分析的数值模拟法

首先对响应面法进行了改进并应用于弹性连杆机构刚度可靠性分析，以迭代的格式和选择最优插值点的响应面法确定弹性连杆机构刚度可靠性分析的极限状态函数，编制了相应的有限元程序。然后在考虑安全、失效状态模糊性的基础上建立了弹性连杆机构的广义可靠性分析模型；提出了以重要抽样法与描述性抽样法相结合来求解弹性连杆机构广义失效概率的方法，此方法可以大量减少抽样时间，并且可以大大提高抽样效率，从而加快结果的收敛速度。     

极限分析的无搜索数学规划算法

本文研究理想刚塑性介质极限载荷因子的计算方法。根据极限分权理论的上限定理,建立了计算极限载荷因子的一般数学规划有限元格式。针对这种格式的特点,提出了一个求解极限载荷因子的无搜索迭代算法。这个算法中采用逐步识别刚性、塑性分区,不断修正目标函数的方案,克服了目标函数非光滑所导致的困难。本文提出的算法建立于位移模式有限元基础上,有较广的适用范围,且具有计算效率高,稳定性好,格式简单易于程序实现等优点。     

A SOLUTION PROCEDURE FOR LOWER BOUND LIMIT AND SHAKEDOWN ANALYSIS BY SGBEM

The symmetric Galerkin boundary element method (SGBEM) instead of the finite element method is used to perform lower bound limit and shakedown analysis of structures. The self-equilibrium stress fields are constructed by a linear combination of several basic self-equilibrium stress fields with parameters to be determined. These basic self-equilibrium stress fields are expressed as elastic responses of the body to imposed permanent strains and obtained through elastic-plastic incremental analysis. The complex method is used to solve nonlinear programming and determine the maximal load amplifier. The limit analysis is treated as a special case of shakedown analysis in which only the proportional loading is considered. The numerical results show that SGBEM is efficient and accurate for solving limit and shakedown analysis problems. Project supported by the National Natural Science Foundation of China (No. 19902007), the National Foundation for Excellent Doctorial Dissertation of China (No. 200025) and the Basic Research Foundation of Tsinghua University.     

Upper Bound Shakedown Analysis of Plates Utilizing the C1 Natural Element Method

This paper proposes a numerical solution method for upper bound shakedown analysis of perfectly elasto-plastic thin plates by employing the C1 natural element m...     

含缺陷轴对称体的安定与极限分析

利用静力安全定理得到了计算轴对称体安定与极限载荷的统一格式，采用温度参数法构造安定分析所需的残余应力场，为了克服对工程实际问题进行安定分析时解题规模与计算精度的矛盾，针对轴对称体的特点，采用两种线性化方案对屈服面进行线性化处理，即直接内接法和在降维应力偏量空间中对屈服面的线性化处理，使安定分析转化为一线性规划问题，在简化过程中合理选择线性化方案以便使应力校核点接近精确的屈服面；为了减小计算量，在求     

Shakedown analysis for a class of strengthening materials within the framework of gradient plasticity

The classical shakedown theory is extended to a class of perfectly plastic materials with strengthening effects (Hall–Petch effects). To this aim, a strain gradient plasticity model previously advanced by Polizzotto (2010) is used, whereby a featuring strengthening law provides the strengthening stress, i.e. the increase of the yield strength produced by plastic deformation, as a degree-zero homogeneous second-order differential form in the accumulated plastic strain with associated higher order boundary conditions. The extended static (Melan) and kinematic (Koiter) shakedown theorems are proved together with the related lower bound and upper bound theorems. The shakedown limit load problem is addressed and discussed in the present context, and its solution uniqueness shown out. A simple micro-scale structural system is considered as an illustrative example. The shakedown limit load is shown to increase with decreasing the structural size, which is a manifestation of the classical Hall–Petch effects in a context of cyclic loading.     

弹塑性结构安定下限分析的无网格局部

将基于Voronoi结构的无网格局部Petrov-Galerkin法与减缩基技术相结合,建立了一种安定下限分析的新方法．为了克服移动最小二乘近似难以准确施加本质边界条件的缺点,采用了自然邻近插值构造试函数.通过引入基准载荷域上载荷角点的概念,消除了安定下限分析中由时间参数所引起的求解困难．利用减缩基技术,将安定分析问题化为一系列未知变量较少的非线性规划子问题．在每个非线性规划子问题中,自平衡应力场由一组带有待定系数的自平衡应力场基矢量的线性组合进行模拟,而这些自平衡应力场基矢量可应用弹塑性增量分析中的平衡迭代结果得到．算例结果证明了提出的分析方法的有效性．      

Shakedown analyses for rolling and sliding contact problems

There is a range of problems where repeated rolling and sliding contact occurs over a half space of an elastic–perfectly plastic material. For such problems shakedown and limit analysis provide significant advantages over other forms of analysis when a global understanding of deformation behaviour is required. In this paper, a recently developed numerical upper bound method, the Linear Matching Method (LMM), for shakedown analyses is applied to the solution of a problem previously considered by Ponter et al. [Ponter, A.R.S., Hearle, A.D., Johnson, K.L., 1985. J. Mech. Phys. Solids 33 (4), 339–362] for a moving Hertzian contact, with sliding friction. This semi-analytic solution is an upper bound based on certain specific kinematic assumptions. We show that the Ponter, Hearle and Johnson solution is a reasonable approximate solution for a circular contact area but is less accurate for an elliptic contact area. For an elliptic contact area LLM solutions converge to the line contact solution. The effect of the non-coincidence of the direction of travel and slide is also investigated.     

Lower bound shakedown analysis by the symmetric Galerkin boundary element method

In this paper, the static shakedown theorem is reformulated making use of the symmetric Galerkin boundary element method (SGBEM) rather than of finite element method. Based on the classical Melan’s theorem, a numerical solution procedure is presented for shakedown analysis of structures made of elastic-perfectly plastic material. The self-equilibrium stress field is constructed by linear combination of several basis self-equilibrium stress fields with parameters to be determined. These basis self-equilibrium stress fields are expressed as elastic responses of the body to imposed permanent strains obtained through elastic–plastic incremental analysis. The lower bound of shakedown load is obtained via a non-linear mathematical programming problem solved by the Complex method. Numerical examples show that it is feasible and efficient to solve the problems of shakedown analysis by using the SGBEM.     

Limit and shakedown analysis under hydrogen embrittlement condition

A modified shakedown theorem and its solving technique are presented to involve hydrogen embrittlement of steel into limit and shakedown analysis. Firstly, the shakedown theorem for hydrogen embrittled material is derived from a limited kinematic hardening shakedown theorem and hydrogen enhanced localized plasticity mechanism of hydrogen embrittlement. In the presented theorem, hydrogen’s effect is taken into account by the synergistic action of both strength reduction and stress redistribution. Secondly, a novel solving technique is developed based on the basis reduction method, in which the complicated constraints in the resulting nonlinear mathematical programming are released. At last, three numerical examples are carried out to verify the performance of the proposed method and to reveal hydrogen’s effect on the limit and shakedown load of structure. The numerical results are discussed and compared with those from literatures, which proves the accuracy and high efficiency of the introduced solving technique. It is concluded that the proposed theorem can predict the limit and shakedown load of hydrogen embrittled structure reasonably.     

Bounds for shakedown of cohesive-frictional materials under moving surface loads

In this paper, shakedown of a cohesive-frictional half space subjected to moving surface loads is investigated using Melan’s static shakedown theorem. The material in the half space is modelled as a Mohr–Coulomb medium. The sliding and rolling contact between a roller and the half space is assumed to be plane strain and can be approximated by a trapezoidal as well as a Hertzian load distribution. A closed form solution to the elastic stress field for the trapezoidal contact is derived, and is then used for the shakedown analysis. It is demonstrated that, by relaxing either the equilibrium or the yield constraints (or both) on the residual stress field, the shakedown analysis leads to various bounds for the elastic shakedown limit. The differences among the various shakedown load factors are quantitatively compared, and the influence of both Hertzian and trapezoidal contacts for the half space under moving surface loads is studied. The various bounds and shakedown limits obtained in the paper serve as useful benchmarks for future numerical shakedown analysis, and also provide a valuable reference for the safe design of pavements.