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

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

动载荷下弹塑性随动强化结构的安定问题

本文讨论了弹塑性随动强化结构的动力安定问题,其加载面在应力空间中是任意形状的凸曲面。文中给出了相应的安定定理,指出如果能找到任何一个与时间无关的和虚设残余应力而使结构安定,则该结构将一定是动力安定的。此外,本文还给出了塑性功和塑性应变的上界估计式,并以在交变载荷作用下的Mises随动强化材料的简支圆板为例,对定理的应用进行了说明。     

循环接触下安定状态问题的研究

基于线性随动强化理论,运用算子分离技术,研究将弹塑性问题转换为弹性问题和残余问题的分析方法,且针对循环载荷接触安定状态,建立了计算机分析程序,该研究能够分析计算弹塑性接触载荷在安定状态下的应力、残余累积应变及残余应力,分析计算了不同载荷的安定状态,并探讨其残余应力场的分析方法。     

应变强化模型的安定准则研究

基于Melan经典的安定理论和von Mises屈服准则,建立了塑性应变强化条件下结构安定的数学模型,根据与时间无关的应力场的特性,对结构中与时间无关的应力场进行了合理的数学变换,将其与载荷变化系数联系起来,推导出与其对应的结构安定极限范围的表达式,给出塑性应变强化模型安定性存在的简化条件.该结论有利于简化应变强化条件下结构的安定分析.     

考虑夹杂相互作用的复合陶瓷夹杂界面的断裂分析

复合材料中夹杂含量较高时,夹杂间的相互作用能显著改变材料细观应力应变场分布,基体和夹杂中的平均应力应变水平也会发生较大变化,导致复合材料强度等力学性能发生显著变化. 为修正单一夹杂模型运用在实际材料中的误差,基于相互作用直推估计法,建立一种考虑含夹杂相互作用的夹杂界面裂纹开裂模型. 首先根据相互作用直推估计法,得到残余应力和外载应力共同作用下夹杂中的平均应力,再计算无限大基体中相同的夹杂达到相同应力场时的等效加载应力,将此加载应力作为含界面裂纹夹杂的等效应力边界条件,在此边界条件下求得界面裂纹尖端的应力强度因子,进而得到界面裂纹开裂的极限加载条件,并分析了夹杂弹性性能、含量、热残余应力、夹杂尺寸等因素对界面裂纹开裂条件的影响. 结果表明,方法能够有效修正单夹杂模型运用在实际材料中的误差,较大的残余应力对界面裂纹开裂有重要的影响,夹杂刚度的影响并非单调且比较复杂;在残余应力较小时,降低柔性夹杂刚度或者增大刚性夹杂刚度都有利于提高材料强度;扩大夹杂尺寸将导致裂纹开裂极限应力显著降低,从而降低材料强度.      

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

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

棘轮安定曲线在不锈钢压力容器应变强化技术中的应用

应变强化是不锈钢压力容器结构实现轻型化的重要途径,而应变强化内压的确定则是应变强化技术的核心。为了能够更准确有效地达到结构应变强化的目的,对循环加载的应变强化方式进行了研究。通过304不锈钢材料的室温单轴棘轮试验,建立了应力比R0条件下的棘轮安定曲线。根据Mises等效原理,利用该曲线通过一次弹塑性有限元分析直接获得结构在循环载荷作用下的强化内压和产生的塑性应变,与试验结果吻合较好,说明运用循环加载的应变强化方式可以有效地达到应变强化的目的。在达到相同的应变强化程度要求下,该方法降低了强化内压,因此可以减小过载加压的风险。     

非线性随动强化条件下的安定定理

基于经典的安定理论与随动强化模型的一般性质,将结构在强化过程中的背应力计入Von Mises屈服准则,建立了随动强化条件下结构的静力安定定理;将背应力与对应的塑性应变率的点积在一个载荷循环内的积分计入塑性耗散功,建立了随动强化条件下结构的机动安定定理,扩展了经典安定理论的应用范围。针对两种定理的存在格式进行了理论证明,并以推论形式给出了结构在随动强化条件下静力安定和机动安定另外两种存在格式。结果表明,随动强化材料的安定状态和安定极限不受强化过程的影响,只取决于材料的初始屈服应力和最终屈服应力。     

关于安定定理的推理

1.静力安定定理的推理对于理想弹塑性物体,若以σ_(ij)、ε_(ij)表示该物体处于弹塑性状态时的实际应力和实际应变,以σ_(ij)~e、ε_(ij)~e表示在加载过程中相应于完全弹性体中的应力和应变,以σ_(ij)~r、ε_(ij)~r表示其实际的残余应力和残余应变,则有     

弹塑性板壳结构非线性有限元分析

本文根据塑性流动理论的基本公式，由隐式积分导出了与路径无关的变量更新算法和一致切线模量。采用单元广义应力应变直接离散塑性流动定律，构造了杂交应力单元一致切线刚度矩阵的显式表达式，编制了结构有限元程序ＳＡＦＥ，数值算例表明：本文的计算方法和计算程序是正确可靠的，可用于弹塑性板壳结构的非线性分析，计算结果屈曲临界载荷和极限承载能力。     

Shakedown analysis of engineering structures with limited kinematical hardening

We present a numerical method for the computation of shakedown loads of engineering structures with limited kinematical hardening under thermo-mechanical loading. The method is based on Melan’s statical shakedown theorem, which results in a nonlinear convex optimization problem. This is solved by an interior-point algorithm recently developed by the authors, specially designed for lower bound shakedown analysis of large-scale problems. Limited kinematical hardening is taken into account by use of a two-surface model, such that both alternating plasticity and incremental collapse can be captured. For the yield surface as well as for the bounding surface the von Mises criterion is used. The proposed method is validated by two examples, where numerical results are compared to those of literature where available.     

Theoretical and computational aspects in the shakedown analysis of finite elastoplasticity

A fully nonlinear shakedown analysis is considered for structures undergoing large elastic-plastic strains. The underlying kinematics of finite elastoplasticity are based on the multiplicative decomposition of the deformation gradient into elastic and plastic parts. It is Shown that the notion of a fictitious, self-equilibrated residual stress field of Melan's linear shakesdown theorem has to be replaced by the notion of real, self-equilibrated residual state. Path-dependent and path-independent shakedown theorems are presented that can be realized in an incremental step-by-step procedure using Finite Element codes. The numerical implementation is considered for highly nonlinear truss structures undergoing large cyclic deformations with ideal-plastic, isotropic and kinematic hardening material behavior. Path-dependency of the residual states in the case of non-adaptation and path-independency in the case of shakedown are shown, and the shakedown domain is determined taking into account also the stability boundaries of the structure.     

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.     

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.     

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

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

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.