On the elastoplastic cyclic analysis of plane beam structures using a flexibility-based finite element approach

This paper presents the extension of a flexibility-based large increment method (LIM) for the case of cyclic loading. In the last few years, LIM has been successfully tested for solving a range of non-linear structural problems involving elastoplastic material models under monotonic loading. In these analyses, the force-based LIM algorithm provided robust solutions and significant computational savings compared to the displacement-based finite element approach by using fewer elements and integration points. Although in cyclic analysis a step-by-step solution procedure has to be adopted to account for the plastic history, LIM will still have many advantages over the traditional finite element method. Before going into the basic idea of this extension, a brief discussion regarding LIM governing equations is presented followed by the proposed solution procedure. Next, the formulation is specified for the treatment of the elastic perfectly plastic beam element. The local stage for the beam behavior is discussed in detail and the required improvement for the LIM methodology is described. Illustrative truss and beam examples are presented for different non-linear material models. The results are compared with those obtained from a standard displacement method and again highlight the potential benefits of the proposed flexibility-based approach.     

Elastoplastic Analysis with Limited Plastic Deformations and Displacements*

ABSTRACT

This paper presents solution methods for elastoplastic and shakedown analysis of linearly elastic, perfectly plastic bodies for which the conventional classical formulations of these problems are completed by constraints on overall plastic deformation and elastoplastic displacement. The methods are described in terms of nonlinear mathematical programming and provide solutions when the plastic reserves of the body are not fully exhausted, and the plastic performance and the plastic deformations are controlled. Application of the method is illustrated by an example.     

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

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

基于数字散斑相关实验测量的材料构型力的计算方法

材料构型力学主要研究材料中的缺陷(夹杂、空穴、位错、裂纹、塑性区等)的构型(形状、尺寸和位置)改变时,所引起的系统自由能的变化。本研究将基于数字散斑相关技术,实验测量材料试件的位移场分布,随后通过材料构型力的定义式,计算求得弹塑性材料中缺陷构型力的分布。其方法概括如下:位移场通过数字图像相关技术测得;应变及位移梯度场利用三次样条拟合获得;线弹性材料应力通过简单线弹性本构方程获取,而塑性材料的表面应力场通过Ramberg-Osgood本构方程计算求得;弹塑性应变能密度分布则由应力-应变曲线数值积分获得。该方法对普通弹性材料或者弹塑性材料均适用,可以用于各种不同的缺陷及缺陷群的材料构型力测量。     

非完全相关荷载下钢筋混凝土框架结构体系可靠度分析

目前框架结构体系可靠度分析大多是在荷载完全相关的假定下完成的,对于非完全相关荷载作用下的体系可靠度研究则鲜有涉及,对于非理想弹塑性或弹脆性的结构尤其如此.本文在随机系统分析的概率密度演化理论的基础上,结合结构非线性全过程分析的位移控制算法,推导了基于归一化位移参数的结构非线性发展概率密度演化方程和承载力裕度的概率密度演化方程.采用纤维梁柱单元进行非线性分析,研究了钢筋混凝土框架结构在非完全相关荷载下的体系可靠度,并与Monte Carlo法进行了对比分析,证明了文中建议方法的可行性.基于文中方法,分析了荷载相关性对可靠度的影响,计算结果表明,荷载相关性对体系可靠度有比较明显的影响.     

基于等效静态载荷法的非线性动态拓扑优化

等效静态载荷法是目前解决动态拓扑优化问题的一种有效方法,但是由于未考虑塑性变形中应力应变关系的时效性,会导致等效载荷偏大和载荷方向偏差,因此无法有效处理材料的非线性问题.据此,提出一种基于模量比率因子的等效静态载荷计算方法.首先,根据结构应力、应变和杨氏模量之间的函数关系,计算各时刻的计算杨氏模量;然后,得到各时刻的模量比率因子,根据位移等效原则计算修正等效静态载荷,进而进行多工况静态拓扑优化,并更新设计变量再次进行动力学分析,重复直至收敛.数值算例表明,修正后的等效静态载荷法能够解决弹塑性材料的非线性问题,在提高等效精度的同时拓展了在非线性材料方面的应用边界.工程算例表明,在相同约束下,采用基于模量比率因子的等效静态载荷法在轻量化的同时能得到更小的塑性应变.     

无额外自由度广义有限元非线性分析

将无额外自由度的广义有限元法由线弹性分析扩展到弹塑性大变形分析.局部强化函数的构建依赖于已有节点,不引入额外自由度,避免了线性相关性问题.在更新拉格朗日框架下,通过控制方程弱形式的线性化推导得到了节点内力的率形式,并分为材料和几何两部分.考虑超弹性和亚弹-塑性两种材料模型,采用Newton-Raphson迭代求解,给出了相关的一致切线刚度阵.三个典型算例的数值结果表明,本文发展的非线性无额外自由度广义有限元方法不仅能够准确求解超弹性和弹塑性大变形问题,同时相比于传统的线性有限元方法具有更高的精度.本文工作进一步拓宽了无额外自由度广义有限元方法的应用领域.     

用摄动法分析准静态载荷下弹塑性结构的安定问题

本文借助于摄动法用分段线性加载面讨论了弹塑性结构在准静态载荷下的安定问题,给出了一个普遍适用的不等式,由此导出了广义的米兰(Melan,1938)安定准则以及位移与塑性变形的界,并举例进行了说明。     

Identification of parameters of models of nonlinear deformation of isotropic and composite materials on the basis of calculations and experiments aimed at analyzing the dynamic behavior of cylindrical metal–plastic shells

A method for identification of material parameters of the constitutive relations of elastoplastic and viscoelastic deformation of isotropic and composite materials is developed. The method is based on minimizing the functional of the residue of results of numerical and experimental analysis of unsteady deformation of structural elements made of examined materials. The method is tested, and prospects of its application for determining material parameters of viscoelastic and elastoplastic models of nonlinear deformation of cylindrical metal–plastic shells under explosive loading are demonstrated.     

弹塑性力学问题的虚单元法

具有有限差分法特征的虚单元法，可视为是有限元法向任意多边形单元的扩展。在材料细观力学性能表征、非均质材料力学分析等非线性问题方面，传统的弹塑性有限元法具有网格数目多、效率低下等不足之处，而虚单元法使网格划分更加灵活，为材料的弹塑性力学分析等非线性问题提供了新的思路。基于增量法弹塑性力学原理和双线性投影算子，建立了弹塑性力学问题的虚单元法求解技术，提出了弹塑性力学问题虚单元法的应力更新方案，研究了弹性力学问题虚单元法的精度和收敛性，讨论了虚单元法求解弹塑性力学问题的网格依赖性。同时，开展了任意多边形和凹多边形单元的数值试验研究，结果表明，虚单元法无须分割多边形，仅需节点自由度便可求得单元刚度矩阵和应力等效荷载，程序实现简单，计算精度高，改善了传统有限元的网格依赖性和塑性区的网格奇异性。     

空间杆系与叠层板组合结构的非线性稳定分析

本文用非线性有限元法研究一种空间杆系与复合材料叠层板组合结构的非线性稳定问题。将杆系离散成空间梁柱单元，将复合材料叠层板离散成复合材料叠层板单元，并考虑了几何非线性和板的成层正交各向异性性质。为了获得较高的计算效率并顺利通过稳定临界点，求解非线性平衡方程采用了变增量步长和柱面弧长法迭代策略。在考证了程序的正确性之后，具体分析了几个工程模型，得出了于工程设计有益的结论。     

有初始几何缺陷的壳体弹塑性大变形分析

本文在等温小变形弹塑性内时本构方程偏量形式的基础上，导出了适用于大位移、大转动、小应变分析的弹塑性内时本构方程，进一步推导出了带有初始几何缺陷的非线性弹塑性问题的有限元方程，可用于分析缺陷对结构非线性弹塑性反应的影响，也可用于带缺陷的非线性问题求解及稳定性分析．     

Variational multiscale analysis of elastoplastic deformation using meshfree approximation

A variational multiscale method has been presented for efficient analysis of elastoplastic deformation problems. Severe deformation occurs in plastic region and leads to high gradient displacement. Therefore, solution needs to be refined to properly capture local deformation in plastic region. In this work, scale decomposition based on variational formulation is presented. A coarse scale and a fine scale are introduced to represent global and local behavior, respectively. The displacement is decomposed into a coarse and a fine scale. Subsequently the problem is also decomposed into a coarse and a fine scale from the variational formulation. Each scale variable is approximated using meshfree method. Adaptivity can easily and nicely be implemented in meshfree method. As a method of increasing resolution, extrinsic enrichment of partition of unity is used. Each scale problem is solved iteratively and conversed results are obtained consequently. Iteration procedure is indispensable for the elastoplastic deformation analysis. Therefore iterative solution procedure of each scale problem is naturally adequate. The proposed method is applied to the Prandtl’s punch test and shear band problem. The results are compared with those of other methods and the validity of the proposed method is demonstrated.     

考虑支承的大跨度索膜结构整体分析与设计

基于修正的拉格朗日法,建立非线性有限元平衡方程,介绍了空间梁单元非线性刚度矩阵。针对此类结构在形态分析及设计中的难点问题,提出利用定力单元、定形单元和微调单元建立数值模型。采用AutoCAD二次开发工具ObjectARX编制有限元软件,能够进行考虑支承的大跨度索膜结构整体分析与设计。以芜湖体育场为工程背景,进行了形态分析和长期荷载组合下的反应。计算结果表明,本文所提方法和有限元软件显著地提高了设计工作的效率,并能准确地获得满足设计要求的结构初始形态及受力分析数据,便于实际工程的设计应用。     

考虑局部切向约束接触问题的直接刚度法

针对工程中大量存在的切向滑移受到约束的接触问题,提出了基于Lagrange乘子的点-面及点-点接触直接刚度法,该方法生成的接触协调条件可以直接组装到结构刚度阵中,从而可直接用于考虑存在结构几何非线性及材料非线性的接触问题分析中。采用这一算法进行了考虑拱坝横缝张合效应的地震响应分析,结果表明考虑切向滑移约束对接触中的张开度有明显影响。     

Multi-scale dual-phase cyclic plasticity with quantitative transitions and size effects of layered structures

A method is developed for cyclic elastoplastic analysis across micro/meso/macro scales which is effective for the quantitative transition of physical variables and for evaluating the size effects of microstructures. By using the improved self-consistent scheme proposed by Fan^[1] and carrying out a delicate mesoscopic analysis based on a shear-lag model, the size effects including the thickness of hard and soft layers relative to the inclusion dimension are obtained on the overall elastoplastic responses of materials up to 50 cycles. The dominant characteristics of the analysis are that the characteristic dimensions of a microstructure such as the thickness of the layers and the inclusion dimension can be explicitly incorporated into the formulation. Results of numerical analysis using only 4 plastic constants show that the thicker the layer relative to the inclusion size, the softer the material in producing more plastic strain values for a given applied stress amplitude. This is in agreement with the well-known experimental rule that the yield strength of layered structures is inversely proportional to the square root of the spacing between layers. It is found that ratcheting depends very much on the size of the layered-structure and that the thinner the relative thickness of the layer the less the ratcheting displacement. This finding may be used to explain why phenomenological models on ratcheting are not quite successful so far, indicating the significance of across scale analysis in understanding issues which have existed for a long time.     

基于近场动力学的薄板弯曲大变形及断裂分析

薄板结构仅在较小的荷载下就能产生较大的位移、旋转,甚至引发结构产生裂纹并扩展,进而发生结构整体断裂,因此,建立薄板结构在大变形过程中的裂纹扩展及断裂仿真模型,具有重要的工程实际意义.文章建立了用于薄板结构几何大变形和断裂分析的近场动力学(PD)和连续介质力学(CCM)耦合模型.首先基于冯·卡门假设,采用更新的拉格朗日法得到薄板在几何大变形增量步下的虚应变能密度增量公式,并利用虚功原理和均质化假设求出几何大变形微梁键的本构模型参数;接着分别建立几何大变形薄板PD模型与CCM模型的虚应变能密度增量,并建立了薄板几何大变形PD-CCM耦合模型;最后模拟了薄板结构在横向变形作用下的渐进断裂过程,得到与实验结果高度一致的仿真结果,验证了所提出的几何非线性PD-CCM耦合模型的精度.结果表明:本文所提出的薄板PD-CCM耦合模型具有简单高效,无需考虑材料参数限制和边界效应的特点,可以很好地用于预测薄板结构在几何大变形过程中的局部损伤和结构断裂,有利于薄板结构的断裂安全评价和理论发展.     

On the plastic buckling of curved carbon nanotubes

This research, for the first time, predicts theoretically static stability response of a curved carbon nanotube(CCNT) under an elastoplastic behavior with several boundary conditions. The CCNT is exposed to axial compressive loads. The equilibrium equations are extracted regarding the Euler–Bernoulli displacement field by means of the principle of minimizing total potential energy.The elastoplastic stress-strain is concerned with Ramberg–Osgood law on the basis of deformation and flow theories of plasticity. To seize the nano-mechanical behavior of the CCNT, the nonlocal strain gradient elasticity theory is taken into account. The obtained differential equations are solved using the Rayleigh–Ritz method based on a new admissible shape function which is able to analyze stability problems. To authorize the solution, some comparisons are illustrated which show a very good agreement with the published works. Conclusively, the best findings confirm that a plastic analysis is crucial in predicting the mechanical strength of CCNTs.     

功能梯度材料结构分析的半解析数值方法研究

一种典型的半解析数值方法——线法被引入功能梯度材料的结构分析。首先推导了功能梯度材料位移形式的平衡方程和边界条件,然后阐述了线法功能梯度材料结构分析的基本步骤和数值原理。该方法的基本思想是通过有限差分将问题的控制方程半离散为定义在沿梯度方向离散节线上的常微分方程组,然后应用B样条函数Gauss配点法求解该常微分方程组得到问题的解答。为演示线法在功能梯度材料结构分析中的应用,给出了线性梯度和指数梯度功能梯度材料板分别受恒定位移、均匀拉伸载荷和弯曲载荷作用的数值算例。与相应问题解析解和其他数值方法的比较表明,线法的计算结果具有很高的精度,而且不需要任何特殊的考虑就能够有效模拟材料内部物性参数的连续变化,也无需事先选取满足特定条件的待定场函数,是一种非常适合功能梯度材料结构形式和材料特点的半解析数值方法。     

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

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