基于正交基无单元Galerkin法和非线性规划的安定分析方法

基于安定分析的下限定理,用正交基无单元Galerkin法建立了交交载荷作用下理想弹塑性结构安定分析的下限计算格式.在给定载荷域的载荷角点所对应的载荷作用下,采用正交基无单元Galerkin法计算相应的虚拟弹性应力场.并且利用结构在正交基无单元Galerkin法弹塑性增量分析中平衡迭代结果计算得到自平衡应力场基矢量,然后由这些基矢量的线性组合模拟自平街应力场.安定分析问题最终被归结为一系列未知变量较少的非线性数学规划子问题,通过复合形法求解.算例表明该方法的计算结果是令人满意的,并且对初始复合形顶点和用于构造自平衡应力场基矢量的载荷增量是非常不敏感的.     

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

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

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

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

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

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

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

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

非协调元在复合材料安定下限分析中的应用

采用应力函数法,结合均匀化理论和应变法,在细观层次上研究了复合材料的极限和安定分析,获取复合材料代表性体积元在载荷加载历史未知下的容许承载域。利用8节点非协调等参元离散结构,获取弹性应力场和自平衡残余应力场,建立复合材料在细观层次上安定下限的优化格式。在满足计算精度的同时,大大降低了优化规模。以周期性纤维增强金属基复合材料为例,验证了该单元在安定下限分析中的有效性和可靠性。     

薄板结构的安定性分析

本文关于薄板的分析采用样条函数插值及结构中的残余应力场用温度参数模拟，解决了薄板结构在交变载荷作用下的安定问题。由于温度参数的引用，使得研究的问题的自由度大幅度减少，文中还对薄板的Ｍｉｓｅｓ型屈服函数进行了线性化处理，它是安定定理的求解线性化的基础，算例表明本文所提出的方法及程序的可行和有效。     

固支圆板安定性的实验研究

本文利用机动法求解固支圆板在循环集中力作用下的安定问题,得到了安定载荷的理论解;利用激光全息的实验方法研究固支圆板的塑性累积现象,测得了安定载荷的上下限,并对实验结果进行了分析研究.     

弹塑性疲劳微弯裂纹尖端塑性区问题

主要研究疲劳载荷作用下弹塑性裂纹弯曲延伸扩展问题.通过分析论证,比较精确地研究了疲劳载荷作用下弯曲延伸裂纹尖端塑性区域边界上交变应力的分布状况.综合考虑了疲劳作用应力,塑性区域交变应力,利用二阶摄动方法,研究计算了弯曲延伸裂纹尖端塑性区域的范围,并预测了疲劳载荷作用下弹塑性裂纹扩展路径.     

扭力轴三维裂纹扩展寿命仿真研究

对疲劳载荷谱作用下三维表面裂纹,采用双重边界元理论求解裂纹前沿的应力应变场,运用Forman理论、最小应变能密度法和Elber模型,计算裂纹前沿各点的扩展长度、扩展方向和应力强度因子等特征量.根据增量步下裂纹几何形状的改变,对裂纹面进行网格重划分和迭代计算,模拟三维裂纹的扩展,预测裂纹扩展寿命.扭力轴表面裂纹扩展的仿真结果表明该方法合理可行.     

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.     

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.     

Shakedown reduced kinematic formulation,separated collapse modes,and numerical implementation

Our shakedown reduced kinematic formulation is developed to solve some typical plane stress problems, using finite element method. Whenever the comparisons are available, our results agree with the available ones in the literature. The advantage of our approach is its simplicity, computational effectiveness, and the separation of collapse modes for possible different treatments. Second-order cone programming developed for kinematic plastic limit analysis is effectively implemented to study the incremental plasticity collapse mode. The approach is ready to be used to solve general shakedown problems, including those for elastic–plastic kinematic hardening materials and under dynamic loading.     

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.     

On shakedown of three-dimensional elastoplastic strain-hardening structures

The present article considers the shakedown problem of structures made of either kinematic or mixed strain-hardening materials. Some basic and useful shakedown properties of elastoplastic strain-hardening structures are proved mathematically. It is impossible for a kinematic strain-hardening structure to be involved in incremental plastic collapse, and so its only possible failure mode is that of alternating plasticity. A time-independent self-equilibrium stress field has no influence on the shakedown of a kinematic strain-hardening structure although it contributes to the magnitude of plastic deformation. The sufficient shakedown conditions for either kinematic or mixed strain-hardening structures are deduced, from which the lower bound of shakedown load domain can be obtained via a mathematical programming problem. It should be pointed out that, to guarantee the safety of an elastoplastic strain-hardening structure, the damage analysis is also necessary to determine the maximum load factor the structure can bear. The shakedown analysis of strain-hardening structures can be simplified by the conclusions obtained in this article, as is illustrated by two simple examples.     

Lower bound limit analysis of three-dimensional elastoplastic structures by boundary element method

IntroductionThelimitanalysisofstructuresisoneofthemostpracticalandusefulbranchesinplasticity .Ithasimportantapplicationbackgroundforproblemssuchasthedeterminationofloadcarryingcapacityandplasticformingofmetal.Thepurposeofthelimitanalysisofstructuresistoprovidereliabletheoreticalbasesforengineeringdesignandsafetyassessment.Asasimplifiedmethodforelastoplasticproblems,limitanalysisneednotrequirethehistoryofloadandcancomputethelimitloadsdirectlyinsteadofelastoplasticincrementalcomputationwhichisus…