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

作为一种基于自然邻近插值的新型无网格法,自然单元法克服了大多数无网格法难以施加本质边界条件的困难．将自然单元法与减缩基技术相结合,建立了一种轴对称结构极限下限分析的数值格式和求解算法．通过不断修正自平衡应力场,轴对称结构极限下限分析可转化为一系列的非线性数学规划子问题,并由复合形法求解．在每个非线性规划子问题中,自平衡应力场表示为一组带有待定系数的自平衡应力场基矢量的线性组合,并且这些自平衡应力场基矢量可由弹塑性增量分析的平衡迭代结果得到．算例结果表明,本文所提的轴对称结构极限下限分析方法行之有效．     

meshless local Petrov-Galerkin method for static shakedown analysis of elasto-plastic structures

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

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

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

边坡稳定分析的虚功率法

基于有限元计算所得应力场的改进极限平衡法,对渗流与有效应力耦合作用强烈的边坡或地基的稳定性分析具有优势.本文提出了边坡稳定分析的虚功率法,即基于极限分析的上限定理,利用机动许可的组合刚体滑动机构和有限元应力场,用滑动机构的速度间断面上的抗滑力功率与滑动力功率的比值计算安全系数.通过分步优化方法,获得边坡给定滑动机制的稳定安全系数.对2个典型折线滑动机构的边坡案例进行了分析,比较了采用静力平衡应力场和静力许可应力场对安全系数的影响.指出基于土体线弹性本构模型所得的有效应力场计算的稳定安全系数,也是边坡稳定安全性的一个不错的度量.算例中本文计算所得的边坡稳定安全系数,与文献中推荐的答案很接近,其滑动机制与有限元强度折减法分析所得的滑动机制基本一致,安全系数也接近,表明本文提出的方法是合理的边坡稳定分析新方法,为边坡和地基的稳定性分析提供了新的选择.     

边坡稳定分析的虚功率法

基于有限元计算所得应力场的改进极限平衡法, 对渗流与有效应力耦合作用强烈的边坡或地基的稳定性分析具有优势. 本文提出了边坡稳定分析的虚功率法, 即基于极限分析的上限定理, 利用机动许可的组合刚体滑动机构和有限元应力场, 用滑动机构的速度间断面上的抗滑力功率与滑动力功率的比值计算安全系数. 通过分步优化方法, 获得边坡给定滑动机制的稳定安全系数. 对2个典型折线滑动机构的边坡案例进行了分析, 比较了采用静力平衡应力场和静力许可应力场对安全系数的影响. 指出基于土体线弹性本构模型所得的有效应力场计算的稳定安全系数, 也是边坡稳定安全性的一个不错的度量. 算例中本文计算所得的边坡稳定安全系数, 与文献中推荐的答案很接近, 其滑动机制与有限元强度折减法分析所得的滑动机制基本一致, 安全系数也接近, 表明本文提出的方法是合理的边坡稳定分析新方法, 为边坡和地基的稳定性分析提供了新的选择.      

变曲率连续滚滑接触安定性的数值方法研究

针对接触表面变曲率的特点,引入局部坐标系,构造出局部坐标下残余应力应变场的分布状态,建立了变曲率连续啮合过程中安定状态残余应力的计算方法。该数值方法将弹塑性问题分解为弹性问题和特征应变决定的残余问题,并采用增量映射方法求解特征应变决定的残余问题,可直接得到接触安定状态下的接触残余应力,并随之进行安定极限的判定。采用该数值方法计算了不同曲率处接触点的安定极限,给出了安定极限与摩擦因数之间的关系,并与有关数值结果相比较,验证了该算法的有效性。     

用MVM准则求解轴对称弹塑性问题

本文采用MVM 屈服准则,用相关联流动法则建立材料的本构关系.对于实际工程中常见的轴对称问题(平面应力、平面应变),进行弹塑性分析,给出求解问题的一组微分方程.采用Prager 假设,给出应力场和位移场.在分析中可以看出:对于平面应变问题,当v≠0.5时,求解应力场的问题是非静定的;当v=0.5或在平面应力问题中,求解应力的问题是静定的,方程组易于求解.通过数值计算考察SD 效应对结构的影响.结果表明,在压缩过程中,SD 效应增强了结构抵抗塑性变形的能力.     

结构安定分析的Galerkin边界元方法

基于Melan静力安定定理,利用Galerkin边界元方法建立了多组交变载荷作用下结构安定分析的下限计算格式.在给定载荷域的载荷角点所对应载荷作用下,采用Galerkin边界元法计算相应的虚拟弹性应力场,并且利用结构在Galerkin边界元弹塑性增量计算中同一增量步中不同迭代步之间的应力差作为自平衡应力场的基矢量,通过这些基矢量的线性组合构造了自平衡应力场,大大降低了所形成的数学规划问题的未知变量数.并通过复合形法对非线性规划问题直接进行求解,得到了结构在交变载荷作用下的下限安定乘子.计算结果表明,所采用的方法具有较高的精度和计算效率.     

土质边坡失稳的突变性分析

对土质边坡圆弧滑动稳定性问题，采用应力软化模型，由系统能量导出极限平衡方程，运用突变理论方法对土坡失稳进行了分析。分析表明，边坡参数变化为某一临界值时，土坡的滑动位移突然增大，发生突变失稳现象，由刚体极限平衡稳定理论无法体现应力软化土质边坡的这种失稳的突变性。     

对“分区混合有限元法计算应力强度因子”的讨论

用分区混合有限元法计算应力强度因子,是分区混合能量原理实际应用的一个成功例子。该方法是把裂纹尖端附近作为Ⅰ区,采用一个应力奇异单元,应力场取裂纹尖端附近渐近解的第一项,以应力强度因子作未知量;把其余部分作为Ⅱ区,采用位移型常规单元,以结点位移作为未知量。     

Kinematic limit analysis of frictional materials using nonlinear programming

In this paper, a nonlinear numerical technique is developed to calculate the plastic limit loads and failure modes of frictional materials by means of mathematical programming, limit analysis and the conventional displacement-based finite element method. The analysis is based on a general yield function which can take the form of the Mohr–Coulomb or Drucker–Prager criterion. By using an associated flow rule, a general nonlinear yield criterion can be directly introduced into the kinematic theorem of limit analysis without linearization. The plastic dissipation power can then be expressed in terms of kinematically admissible velocity fields and a nonlinear optimization formulation is obtained. The nonlinear formulation only has one constraint and requires considerably less computational effort than a linear programming formulation. The calculation is based entirely on kinematically admissible velocities without calculation of the stress field. The finite element formulation of kinematic limit analysis is developed and solved as a nonlinear mathematical programming problem subject to a single equality constraint. The objective function corresponds to the plastic dissipation power which is then minimized to give an upper bound to the true limit load. An effective, direct iterative algorithm for kinematic limit analysis is proposed in this paper to solve the resulting nonlinear mathematical programming problem. The effectiveness and efficiency of the proposed method have been illustrated through a number of numerical examples.     

Extended bounding theorems of limit analysis

This paper studies the bounding problems of the complete solu-tion of limit analysis for a rigid-perfectly plastic medium,allowing for the discontinuity of plastic flow.A generalizedvariational principle involving conditions of the rigid-plas-tic interface and the discontinuous surface of a velocityfield has been advanced for the mixed-boundary value problem.Based on this principle,a set of variational formulae of li-mit analysis is established.The safety factors obtained bythese formulae lie between the upper and lower bounds obtainedby the classical bounding theorems with the same kinematicallyand statically admissible field.Moreover,extended bounding theorems have been derivedand proved,which hold a broader stress and velocity field thanthe statically and kinematically admissible field.The corol-laries of these theorems indicate the relationship between thevariational solution and the complete solution of limit analy-sis.Applications of these theorems show that a close approxi-mation can be obtained     

A bipotential approach for plastic limit loads of strip footings with non-associated materials

This paper is concerned with a bipotential approach for estimating the plastic collapse loads of a half-space made with a non-associated Mohr–Coulomb material and indented by a rigid punch. In geotechnics, this problem is called the bearing capacity of shallow strip footing for which the analytical solution is derived by Prandtl (1920) [46] and Hill (1950) [35] in the context of associated plasticity. However, when the plastic model is not associated, no analytical methods have yet been developed. Here we explore this issue in a rigorous mathematical framework coupling the bipotential concept and limit analysis. First, the method proposed makes use of the method of characteristics to build a statically and plastically admissible stress field that enables a lower estimate of the plastic limit loads. Next, the extended kinematic theorem of limit analysis to non-standard plasticity is applied to derive an upper quasi-bound of the collapse loads. For this aim, the internal rate of plastic dissipation is obtained thanks to the bipotential functional depending on both a trial stress field and a Prandtl-like collapse mechanism. The analytic estimates are compared to the formulae and numerical results provided in literature.     

A nonlinear programming approach to kinematic shakedown analysis of frictional materials

This paper develops a novel nonlinear numerical method to perform shakedown analysis of structures subjected to variable loads by means of nonlinear programming techniques and the displacement-based finite element method. The analysis is based on a general yield function which can take the form of most soil yield criteria (e.g. the Mohr–Coulomb or Drucker–Prager criterion). Using an associated flow rule, a general yield criterion can be directly introduced into the kinematic theorem of shakedown analysis without linearization. The plastic dissipation power can then be expressed in terms of the kinematically admissible velocity and a nonlinear formulation is obtained. By means of nonlinear mathematical programming techniques and the finite element method, a numerical model for kinematic shakedown analysis is developed as a nonlinear mathematical programming problem subject to only a small number of equality constraints. The objective function corresponds to the plastic dissipation power which is to be minimized and an upper bound to the shakedown load can be calculated. An effective, direct iterative algorithm is then proposed to solve the resulting nonlinear programming problem. The calculation is based on the kinematically admissible velocity with one-step calculation of the elastic stress field. Only a small number of equality constraints are introduced and the computational effort is very modest. The effectiveness and efficiency of the proposed numerical method have been validated by several numerical examples.     

LOWER BOUND LIMIT ANALYSIS OF THREE-DIMENSIONAL ELASTOPLASTIC STRUCTURES BY BOUNDARY ELEMENT METHOD

IntroductionThelimitanalysisofstructuresisoneofthemostpracticalandusefulbranchesinplasticity .Ithasimportantapplicationbackgroundforproblemssuchasthedeterminationofloadcarryingcapacityandplasticformingofmetal.Thepurposeofthelimitanalysisofstructuresistoprovidereliabletheoreticalbasesforengineeringdesignandsafetyassessment.Asasimplifiedmethodforelastoplasticproblems,limitanalysisneednotrequirethehistoryofloadandcancomputethelimitloadsdirectlyinsteadofelastoplasticincrementalcomputationwhichisus…     

An extension of Neumann’s method for shape control of force-induced elastic vibrations by eigenstrains

This paper is concerned with the force-induced vibrations of linear elastic solids and structures. We seek a transient distribution of actuating stresses produced by additional eigenstrain, such that the vibrations produced by a given set of imposed forces are exactly compensated. This problem, known as dynamic shape control problem in structural engineering, or as dynamic displacement compensation problem in automatic control, is inverse to the usual direct problem of determining displacements due to imposed forces and actuation stresses. In the present paper, we extend a method, which was introduced by F.E. Neumann for demonstrating the uniqueness of direct elastodynamic problems. We use this extended Neumann method in order to show that the distribution of the actuating stresses for shape control must be equal to any statically admissible stress distribution that is in temporal equilibrium with the imposed forces. We furthermore discuss the role of stresses corresponding to this class of solutions in some detail, emphasizing the non-unique nature of a statically admissible stress. As an analytical justification of our formulations, we show that our method reveals some static results by J.M.C. Duhamel and by W. Voigt and D.E. Carlson. Particularly, our method can be interpreted as a dynamic extension of the Duhamel body-force analogy. We moreover present numerical results for a dynamically loaded, irregularly shaped domain in a state of plane strain. These finite element computations give excellent evidence for the validity of the presented method of shape control for both, the case of a step-input and the case of a harmonic excitation.     

极限下限分析的正交基无单元Galerkin法

基于极限分析的下限定理，建立了用正交基无单元Galerkin法进行理想弹塑性结构极 限分析的整套求解算法．下限分析所需的虚拟弹性应力场可由正交基无单元Galerkin法直接 得到，所需的自平衡应力场由一组带有待定系数的自平衡应力场基矢量的线性组合进行模 拟．这些自平衡应力场基矢量可由弹塑性增量分析中的平衡迭代得到．通过对自平衡应力场 子空间的不断修正，整个问题的求解将化为一系列非线性数学规划子问题，并通过复合形法 进行求解．算例表明该方法有效地克服了维数障碍问题，使计算效率得到了充分的提高，是 切实可行的．     

Analytical formulation of stresses in curved composite beams

Summary This paper outlines an analytical method for computing normal and shear stresses generated in a curved laminated beam under bending loads. Each cross section is assumed to be symmetrical and loads are applied in the plane of symmetry. We build a statically admissible stress field in order to plot normal and shear stress distributions. Received 5 March 1997; accepted for publication 18 September 1997     

Limit analysis of composite materials based on an ellipsoid yield criterion

Employing repeating unit cell (RUC) to represent the microstructure of periodic composite materials, this paper develops a numerical technique to calculate the plastic limit loads and failure modes of composites by means of homogenization technique and limit analysis in conjunction with the displacement-based finite element method. With the aid of homogenization theory, the classical kinematic limit theorem is generalized to incorporate the microstructure of composites. Using an associated flow rule, the plastic dissipation power for an ellipsoid yield criterion is expressed in terms of the kinematically admissible velocity. Based on nonlinear mathematical programming techniques, the finite element modelling of kinematic limit analysis is then developed as a nonlinear mathematical programming problem subject to only a small number of equality constraints. The objective function corresponds to the plastic dissipation power which is to be minimized and an upper bound to the limit load of a composite is then obtained. The nonlinear formulation has a very small number of constraints and requires much less computational effort than a linear formulation. An effective, direct iterative algorithm is proposed to solve the resulting nonlinear programming problem. The effectiveness and efficiency of the proposed method have been validated by several numerical examples. The proposed method can provide theoretical foundation and serve as a powerful numerical tool for the engineering design of composite materials.