基于梯度塑性理论的动力软化问题分析

对于动力应变软化问题,采用梯度塑性模型进行分析,该模型能够有效地克服软化材料在有限元分析中的网格依赖性问题。对于动力非线性软化方程的求解则利用于基于参变量变分原理的参数二次规划方法。对于结构动力方程时域上的求解则采用传统的Newmark方法。本文的算法与传统方法相比在求解基于梯度塑性模型的非线性动力软化问题时保证了已有参数二次规划算法的优良特性,有实现简单与稳定性好等优点。给出的数值算例证实了本文的理论工作与所研制程序的正确性。     

随机结构动力可靠度分析的概率密度演化方法

基于随机结构动力反应分析的概率密度演化方法，提出了一类新的随机结构动力可靠度分析方法。在随机结构动力反应概率密度演化方程的基础上，对于首次超越问题，根据所给的首次超越破坏准则施加相应的吸收壁边界条件，求解具有吸收壁边界条件的概率密度演化方程并在安全域内积分．给出结构的动力可靠度。结合精细时程积分方法和具有TVD性质的差分格式，讨论了计算结构动力可靠度的数值方法。以八层框架结构为例进行了动力可靠度分析并与随机模拟分析结果进行了比较。     

动力弹塑性刚体界面元法

本文探讨了用刚体界面元法求解二维动力弹塑性问题的理论．首先，从连续介质的质点运动微分方程出发，利用其等效积分方程，通过引入刚体模式的位移场函数建立了刚体界面元的基本动力方程．为了以渐进解法求解弹塑性问题，将界面元简化的应力矢量引入塑性流动理论，从而给出了界面元的弹塑性本构关系．     

二维正交各向异性结构弹塑性问题的边界元分析

给出了二维正交各向异性结构弹塑性问题的边界元分析方法, 包括相应边界积分方程、内点应力公式、边界元求解格式以及弹塑性应力计算方法。在弹塑性分析中, 引入了Hill-Tsai 屈服准则, 采用初应力法和切向预测径向返回法确定实际应力状态。通过具体算例分析了二维正交各向异性结构的弹塑性应力和塑性区分布情况, 部分数值结果与已有结果进行了比较, 两者基本吻合。结果表明, 本文中给出的边界元法可以有效地用于求解二维正交各向异性结构的弹塑性问题。      

受演变随机激励结构响应的扩展精细积分方法

摘要：对于受演变随机激励的线性多自由度体系,给出了计算其非平稳响应的扩展精细积分方法。首先采用虚拟激励法,将随机荷载转化成确定性荷载,然后采用Duhamel积分的精细计算方法,构造出统一形式的精确、高效递推格式。本文方法避免了矩阵的求逆运算,不依赖于系统矩阵或其动力矩阵的性态,提高了数值稳定性和应用范围。本文方法具有与混合型时程精细积分方法同样高的数值精度,而效率上要高于增维精细积分方法。算例验证了本文算法的优越性。     

非线性动力方程的一种改进精细积分单步方法EI北大核心CSCD

微分求积法和单步块方法都是单步多级数值方法,但是直接应用于求解非线性动力方程时的计算量比较巨大,为此提出了一种基于单步块方法的改进精细积分单步方法。结合精细积分法,该方法采用s级的单步块方法的第s个方程对Duhamel积分项进行数值积分。具体采用四阶Runge-Kutta法获得待求变量的预估值,并采用新四点积分公式计算Duhamel积分项。相对于现有的单步方法,该改进算法在数值精度和稳定性上更优。通过非线性动力方程的典型算例验证了该算法的优势。     

基础隔震结构地震响应时程分析新方法

针对隔震层恢复力模型为折线型的层间剪切动力方程，研究了基础隔震结构在地震作用下的弹塑性时程反应。通过多种方法比较，说明了各种方法的优缺点。采用把结构特性和外部激励分开进行精细积分的增维分块精细积分法，避免了每一步迭代过程中的指数矩阵精细积分求解，在保持计算精度的同时，提高了计算效率；考虑地震作用特点后，增维分块精细积分法的计算效率有所提高。因此，增维分块积分法在计算基础隔震结构地震响应上是非常有效的。     

非线性动力方程的精细时空有限元方法

根据Hamilton变作用定律构造了时空有限元矩阵;并根据传递矩阵原理,利用时间的一维性将时空有限元矩阵变换为时间方向的传递矩阵,将初值问题转化为一般矩阵相乘问题以方便求解。为了保证计算的稳定性,参考了精细积分的思想提出精细时空有限元方法,并给出线性问题在时间级数荷载作用下的计算式。数值分析结果证明该方法在线性问题分析上非常准确并可以推广到非线性动力方程的求解;只需将非线性解看作初始解和增量解的叠加,通过精细时空有限元线性求解方法计算增量解,逐步修正后即可得到非线性解。结果表明该方法是一个有效的求解非线性动力方程的方法。     

复合随机振动系统的动力可靠度分析

建议了一类新的复合随机振动系统动力可靠度分析方法。基于复合随机振动系统反应分析的密度演化方法,根据首次超越破坏准则,对密度演化方程施加相应的边界条件,进而求解密度演化方程,在安全域内积分给出结构的动力可靠度。结合精细时程积分方法与具有TVD性质的差分格式,研究了基于密度演化方法求解结构动力可靠度问题的数值方法。以受到随机地震作用、具有随机参数的八层层间剪切型结构为例,进行了结构动力可靠度分析并与随机模拟结果进行了比较。研究表明,建议的方法具有较高的精度和效率。     

具有输入约束的移动机器人路径跟踪预测控制

针对移动机器人设定速度与实际速度不匹配的问题,提出了一种基于预测控制的移动机器人路径跟踪方法。首先,基于运动学关系建立了移动机器人路径跟踪误差模型,给出了期望路径的参数化方法及更新方程。其次,通过定义与路径参数相关的预测性能指标,并结合状态空间方程给出了预测模型向量描述,得到了具有不等式约束的二次规划优化问题。进而,采用有效集二次规划方法求解具有不等式约束的优化问题获得最优控制量。最后,通过仿真实验分析所提算法的有效性,并设计移动机器人路径跟踪控制系统实验平台验证所提算法的有效性。     

压气机过盈配合的弹塑性有摩擦接触的研究

增压器压气机叶轮、轴套和轴采用过盈配合技术联成一体,这是三维多体弹塑性有摩擦接触问题,是两种非线性相互耦合的边值待定问题。采用有限元参数二次规划法,并结合多重子结构技术,充分利用两种方法各自的长处分析求解柴油机涡轮增压器叶轮与轴套、轴套与轴的三维弹塑性有摩擦接触问题,针对不同的过盈量、摩擦因数、转速和轴套壁厚进行了大量计算,获得了叶轮、轴套与轴之间接触应力的相应分布规律。轴套与轴的装配过盈量是影响轮轴接触应力的重要因素。在选择叶轮、轴套和轴三者装配尺寸时,尤其采用压力组装法时应严格控制轴套与轴的过盈量。研究结果表明,此方法对压气机弹塑性接触特性是有效的,能够反映接触区域的接触法向应力、变形以及摩擦力的大小和分布情况。     

Elasto‐plasticity revisited: numerical analysis via reproducing kernel particle method and parametric quadratic programming

Aiming to simplify the solution process of elasto‐plastic problems, this paper proposes a reproducing kernel particle algorithm based on principles of parametric quadratic programming for elasto‐plasticity. The parametric quadratic programming theory is useful and effective for the assessment of certain features of structural elasto‐plastic behaviour and can also be exploited for numerical iteration. Examples are presented to illustrate the essential aspects of the behaviour of the model proposed and the flexibility of the coupled parametric quadratic programming formulations with the reproducing kernel particle method. Copyright © 2002 John Wiley & Sons, Ltd.     

Parametric quadratic programming method for elastic contact fracture analysis

A solution procedure for elastic contact fracture mechanics has been proposed in this paper. The procedure is based on the quadratic programming and finite element method (FEM). In this paper, parametric quadratic programming method for two-dimensional contact mechanics analysis is applied to the crack problems involving the crack surfaces in frictional contact. Based on a linear complementary contact condition, the parametric variational principle and FEM, a linear complementary method is extended to analyze contact fracture mechanics. The near-tip fields are properly modeled in the analysis using special crack tip elements with quarter-point nodes. Stress intensity factor solutions are presented for some frictional contact fracture problems and are compared with known results where available.     

基于参变量变分原理的三维Cosserat体模型弹塑性分析与应变局部化模拟

基于参变量变分原理,该文发展了三维Cosserat连续体模型弹塑性有限元分析的二次规划算法。由于Cosserat连续体模型的本构方程中存在材料内尺度参数,该模型可以解决经典连续介质理论在分析应变软化问题时病态的有限元网格依赖性问题。数值结果表明所发展的三维Cosserat连续介质弹塑性有限元模型可以有效的模拟应变局部化现象并且该算法具有很好的数值稳定性,同时获得的数值结果具有良好的非网格依赖性。     

A fast first-order optimization approach to elastoplastic analysis of skeletal structures

It is classical that, when the small deformation is assumed, the incremental analysis problem of an elastoplastic structure with a piecewise-linear yield condition and a linear strain hardening model can be formulated as a convex quadratic programming problem. Alternatively, this paper presents a different formulation, an unconstrained nonsmooth convex optimization problem, and proposes to solve it with an accelerated gradient-like method. Specifically, we adopt an accelerated proximal gradient method, that has been developed for a regularized least squares problem. Numerical experiments show that the presented algorithm is effective for large-scale elastoplastic analysis. Also, a simple warm-start strategy can speed up the algorithm when the path-dependent incremental analysis is carried out.     

Global Constitutive Behavior of Periodic Assemblies of Inelastic Bodies in Contact

This work attempts to capture the macroscopic behavior of inelastic bodies in contact by means of a numerical homogenized constitutive relation. The analysis is restricted to small strains, plane problems, and monotonic proportional loadings. An important feature of this work is the quasi-static frictional contact analysis of the microstructure composed of deformable inelastic bodies by means of a parametric quadratic programming principle and its corresponding algorithm in numerical analysis. Two numerical examples are given to demonstrate the efficiency of the algorithm presented in this article.     

Hybrid-Trefftz stress elements for elastoplasticity

The stress model of the hybrid finite element formulation is applied to the analysis of quasi-static, gradient-dependent elastoplastic structural problems. The finite element approximation consists in the direct estimate of the stress and plastic multiplier fields in the domain of the element and of the displacements and plastic multiplier gradients on its boundary. The finite element equations are derived directly from the relevant fundamental structural conditions, namely equilibrium, compatibility, elasticity and gradient-dependent plasticity. The finite element solving system for the finite step incremental analysis is encoded as a recursive sequence of symmetric parametric linear complementarity problems (SPLCP). The sequence of SPLCP is solved using a direct extension of the restricted basis linear programming algorithm. The implementation of the formulation and of the algorithm is illustrated with numerical applications. © 1998 John Wiley & Sons, Ltd.     

A finite element model for contact analysis of multiple Cosserat bodies

The objective of this paper is to develop a finite element model for multi-body contact analysis of Cosserat materials. Based on the parametric virtual work principle, a quadratic programming method is developed for finite element analysis of contact problems. The contact problem with friction between two Cosserat bodies is treated in the same way as in plastic analysis. The penalty factors, that are normally introduced into the algorithm for contact analysis, have a direct influence on accuracy of solution. There is no available rule for choosing a reasonable value of these factors for simulation of contact problems of Cosserat materials, and they are therefore cancelled through a special technique so that the numerical results can be of high accuracy. Compared with the conventional work on Cosserat elasticity, the newly developed model is on the contact analysis of the Cosserat materials and is seldom found in the existing literatures. Four examples are computed to illustrate the validity and importance of the model developed.     

A second order cone complementarity approach for the numerical solution of elastoplasticity problems

In this paper we present a new approach for solving elastoplastic problems as second order cone complementarity problems (SOCCPs). Specially, two classes of elastoplastic problems, i.e. the J ₂ plasticity problems with combined linear kinematic and isotropic hardening laws and the Drucker-Prager plasticity problems with associative or non-associative flow rules, are taken as the examples to illustrate the main idea of our new approach. In the new approach, firstly, the classical elastoplastic constitutive equations are equivalently reformulated as second order cone complementarity conditions. Secondly, by employing the finite element method and treating the nodal displacements and the plasticity multiplier vectors of Gaussian integration points as the unknown variables, we obtain a standard SOCCP formulation for the elastoplasticity analysis, which enables the using of general SOCCP solvers developed in the field of mathematical programming be directly available in the field of computational plasticity. Finally, a semi-smooth Newton algorithm is suggested to solve the obtained SOCCPs. Numerical results of several classical plasticity benchmark problems confirm the effectiveness and robustness of the SOCCP approach.