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.     

基于均匀化理论韧性复合材料塑性极限分析

运用细观力学中的均匀化方法分析了韧性复合材料的塑性极限承载能力.从反映复合材料细观结构的代表性胞元入手,将均匀化理论运用到塑性极限分析中,计算由理想刚塑性、Mises组分材料构成的复合材料的极限承载能力.运用机动极限方法和有限元技术,最终将上述问题归结为求解一组带等式约束的非线性数学规划问题,并采用一种无搜索直接迭代算法求解.为复合材料的强度分析提供了一个有效手段.     

自适应分析在确定裂纹尖端塑性区中的应用

在分析裂纹扩展及材料强度时对塑性区的估算是很重要的,本文提出采用自适应有限元分析来确定裂纹尖端塑性区的方法,计算结果表明这种方法通过网格自动加密能够有效地跟踪出弹塑性的交界面。     

确定复合材料宏观屈服准则的细观力学方法

运用细观力学中的均匀化方法，分析了含周期性微结构复合材料的宏观屈服准则，并对Hill-Tsai准则进行了修正。从基于复合材料细观结构的代表性胞元入手，运用塑性极限理论中的机动分析以及有限元方法，计算了细观结构的极限载荷域。通过宏细观尺度对应关系，得到复合材料的宏观屈服准则。     

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.     

多孔材料塑性极限载荷及其破坏模式分析

运用塑性力学中的机动极限分析理论，研究韧性基体多孔材料的塑性极限承载能力和破坏模式。以多孔材料的细观结构为研究对象，将细观力学中的均匀化理论引入到塑性极限分析中，并结合有限元技术，建立细观结构极限载荷的一般计算格式，并提出相应的求解算法。数值算例表明：细观孔洞对材料的宏观强度影响明显；在单向拉伸作用下，孔洞呈现膨胀扩大规律；多孔材料破坏源于基体塑性区的贯通。     

A general algorithm for plastic flow simulation by finite element limit analysis

Limit analysis has been rendered versatile in many structural and metal forming problems. In metal forming analysis, the slip-line method and the upper bound method have filled the role of limit analysis. As a breakthrough of the previous work, a computational approach to limit solutions is considered as the most challenging area.In the present work, a general algorithm for limit solutions of plastic flow is developed with the use of finite element limit analysis. The algorithm deals with a generalized Hölder inequality, a duality theorem, and combined smoothing and successive approximation in addition to a general procedure for finite element analysis. The algorithm is robust such that from any initial trial solution, the first iteration falls into a convex set which contains the exact solution (s) of the problem. The idea of the algorithm for limit solutions is extended from rigid⧹perfectly plastic materials to work-hardening materials by the nature of the limit formulation, which is also robust with numerically stable convergence and highly efficient computing time.     

基于连续体损伤理论的硼钢高温成形极限确定法

为了应对节能减排及汽车安全性的双重要求,先进高强钢的开发和应用受到愈来愈高的重视. 超高强度钢板热冲压成形技术是减轻车身重量、提高汽车抗冲击和防撞性能的重要途径,目前已经成为世界汽车制造行业的热门技术. 成形极限是硼钢热冲压成形过程中一个非常重要的表征参数,但由于高温半球形凸模胀形试验(Nakazima) 的复杂性造成实验测试面临诸多困难. 针对硼钢热冲压成形过程特点,提出将韧性损伤模型与有限元数值模拟技术结合起来预测硼钢高温成形极限的方法. 在连续体损伤力学基础上,建立了考虑有效应力和等效塑性应变对微孔损伤影响的勒迈特(Lemaitre) 韧性损伤演化方程,在损伤耗散势函数中引入了等效塑性应变因子,从而寻找到能够准确反映硼钢在奥氏体状态下的成形特性的勒迈特耗散势函数. 选取基于遗传算法的基于非支配排序的多目标优化算法(NSGA-II) 与有限元分析软件“FORGETM” 相结合的方法优化韧性损伤模型中材料的损伤因子,并利用优化后的损伤因子准确预测到板料成形过程的拉伸失稳和断裂现象. 将建立的韧性损伤模型引入到硼钢高温半球形凸模胀形有限元仿真模型中,通过预测的成形极限与实验测定结果的比较,验证了所提韧性损伤模型预测硼钢高温成形极限的可靠性.      

板材多点成形过程的有限元分析

多点成形过程采用静力隐式格式进行数值模拟是比较合适的。本文建立了用于多点成形过程分析的静力隐式弹塑性大变形有限元方法 ,给出了对稳定迭代收敛过程效果较好的板壳有限单元模型、处理多点不连续接触边界的接触单元方法以及增量变形过程中应力及塑性应变计算的多步回映计算方法。基于这些方法编制了计算软件 ,应用该软件进行了矩形板的液压胀形过程及球形模具拉伸成形过程的有限元分析 ,数值计算结果与典型的实验结果及计算结果吻合很好。最后给出了球形、圆柱形目标形状的实际多点成形过程的数值模拟结果。     

A mathematical programming algorithm for limit analysis

This paper deals with the limit analyses of perfect rigid-plastic continua. Based on the kinematic theorem of the limit analysis theory, a mathematical programming finite element formula for determining the upper bound load multiplier has been established, and an iteration algorithm proposed accordingly. In this algorithm the plastic and rigid zones are distinguished for every iteration step, and the goal function is modified gradually. The difficulties caused by the nonsmoothness of the goal function are overcome. Some examples solved by this algorithm are presented. The project supported by National Natural Science Foundation of China.     

A plasticity model for cellular materials with open-celled structure

Based on a rigid-plastic material model that obeys the von Mises yield criterion, the plastic behavior of foams with an open-celled structure is studied in this paper using a single unit cell. An approximate continuum plasticity model is developed within the framework of the upper bound theorem of plasticity to describe the yield behavior of foams. The microscopic velocity fields are derived for the unit cell, which satisfy the incompressibility and the kinematic boundary conditions, and expressed in macroscopic rate of deformation. From the microscopic velocity fields, a macroscopic yield function is developed for foams under multi-axial stresses and includes the effects of the hydrostatic stress due to the void presence and growth. The dependency of the derived yield surfaces of foams on their relative densities is studied. The plastic behavior of foams is also studied numerically using the finite element method. The newly developed plasticity model is compared with the finite element analysis results and other available foam models and then correlated with the finite element results.     

STRENGTH PROPERTIES OF JOINTED ROCK MASSES BASED ON THE HOMOGENIZATION METHOD

This paper aims to determine the strength properties of jointed rock masses by means of the homogenization method.To reflect the microstructure of jointed rock masses,a representative element volume (R...     

基于自然单元法的极限上限分析

自然单元法是一种基于离散点集的Voronoi图和Delaunay三角化几何信息,以自然邻近插值为试函数的新型数值方法.相对于一般无网格法中常采用的移动最小二乘近似而言,自然邻近插值不涉及到复杂的矩阵求逆运算,更不需要任何人为的参数,可以提高计算效率.采用该方法构造的形函数满足Delta函数的性质,可以像有限元一样准确地施加边界条件,可以方便处理场函数及其导数的不连续性的问题.论文将自然单元法应用到极限上限分析中,编制了相应的计算程序,通过极限分析的几个经典算例进行了验证,同时采用类似于分片应力磨平的方式,编制相应的磨平程序,由计算点上的塑性耗散功外推得到了节点上的塑性耗散功的值,从而画出了极限状态下结构的塑性耗散功的分布云图.计算结果表明采用自然单元法求解极限上限分析具有稳定性好,精度高,收敛快等优点.     

Mechanical formula for the plastic limit pressure of stent during expansion

The mechanics of cardiovascular stents during the process of expansion are very important for stent function and safety. In general, finite element method （FEM） or experi- ments are major methods used to ascertain mechanical prop- erties of the stent. In this paper, we develop a theoretical model of the tubular stent, derive formulas for the axial forces and moments on the stent end, and propose formu- las for the plastic limit pressure vs. the stent＇s radius during expansion. Examples covering different geometrical param- eters and material parameters are provided, and the plastic limit pressures calculated by FEM and the present method are compared, proving that the present formulas are acceptable and meaningful for the design and innovation of the stent.     

半潜式钻井平台管道钢构支架极限强度研究

钢构支架是半潜式钻井平台管道支吊架的主要类型之一,钢构支架的极限强度是管道系统正常工作的重要保障。研究结构极限强度的方法有理论分析、有限元数值仿真和实验分析。在分析极限理论的基础上由静力法计算四种钢构支架试件的极限载荷;利用MSC软件对试件极限载荷进行数值仿真分析,并对有限元模型化技术进行讨论;对试件进行实验极限载荷测试,对比分析了三种方法测得的极限载荷。结果表明,三种分析方法计算得到结构的极限载荷基本一致,对于结构形式较为简单的结构通过理论分析可以得到简化解析解析解;数值仿真分析中采用合理的有限元模型化技术(结构有限元模型、边界、约束等)可得到精度较高的计算结果。     

Numerical study of the plastic behaviour in tension of welds in high strength steels

The influence of the mismatch between material properties and constraint on the plastic deformation behaviour of the heat affected zone of welds in high strength steels is investigated in this study, using finite element simulations. An elastoplastic implicit three-dimensional finite element code (EPIM3D) was used in the analysis. The paper presents the mechanical model of the code and the methodology used for the numerical simulation of the tensile test of welded joints. Numerical results of the tensile test of welded samples with different hypothetical widths for the Heat Affected Zone and various material mismatch levels are shown. The analysis concerns the overall strength and ductility of the joint and in relation to the plastic behaviour of the heat affected zone. The influence of the yield stress, tensile strength and constraint on the stress and plastic strain distribution in the soft heat affected zone is also discussed.     

等径通道挤压过程三维有限元模拟

利用三维有限元模型对等径通道挤压过程中变形分布的不均匀性进行了模拟,通过对不同摩擦系数下截面等效塑性应变分布比较发现:沿三个方向截面上的变形分布都不是均匀的,这说明二维有限元模型不能真实模拟等径通道挤压过程中试样中的变形分布。此外,摩擦对等效塑性应变分布及挤压力的影响较大,截面变形不均匀参数和最大等效塑性应变出现的位置随摩擦系数的增大而变化。压力-位移曲线的变化可以结合试样的变形过程来解释。     

Non-linear and failure behaviour of spotwelds: a “global” finite element and experiments in pure and mixed modes I/II

The paper deals with joint element model used in crashworthiness simulations. The first part of the paper is dedicated to the formulation of a new “global” finite element for spotweld modelling. The mechanical behaviour of the joint is elastic–plastic type and damage is taken into account to model the failure of the welded area. The second part of the paper concerns a new experimental procedure for joint strength analysis in pure and mixed modes I/II and for joint model characterisation. Experiment is based on Arcan principle and results are compared to open literature. In the last part of the paper, the parameters of the new joint model are identified using experiments and used for several shapes of spot-welded specimens. The model predicts reasonably the elastic–plastic part of the response but is unable to predict the post-peak response observed especially in the case of pure shear.     

Experimental and FEM study of windshield subjected to high speed bird impact

Both experimental and finite element model (FEM) results are presented for the dynamic strength behavior of windshield subjected to bird impact. The experimental data taken from a series of high speed photographs are compared with the numerical results predicted by using FEM in which the windshield was modeled entirely with solid elements and the bird body was approximately simulated by an elastic-plastic material with failure element behavior. Effective plastic strain and element pressure were adopted as the failure criteria and once the pressure or the effective plastic strain of an element reached the critical value, the element would lose the tensile resistance capability completely. The deflection and stress distribution in the windshield were obtained. It is shown that the result from the finite element analysis agrees with those from the full-scale bird impact test. The project supported by the National Natural Science Foundation of China (10272011)