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.     

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

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

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

IntroductionThelimitanalysisofstructuresisoneofthemostpracticalandusefulbranchesinplasticity .Ithasimportantapplicationbackgroundforproblemssuchasthedeterminationofloadcarryingcapacityandplasticformingofmetal.Thepurposeofthelimitanalysisofstructuresistoprovidereliabletheoreticalbasesforengineeringdesignandsafetyassessment.Asasimplifiedmethodforelastoplasticproblems,limitanalysisneednotrequirethehistoryofloadandcancomputethelimitloadsdirectlyinsteadofelastoplasticincrementalcomputationwhichisus…     

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

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

Stress function method in shakedown analysis of axisymmetric shells

A self-equilibrated stress obtained from the stress functions of thin shells is used for the static shakedown theorem as a residual stress. In combination with the finite element method, a linear programming formulation of the shakedown analysis of axisymmetric shells is derived. The physical meaning of the stress function method is clear and its computing amount is small. Some examples of the plates and shells show that the method is reasonable and efficient.The project was supported by the National Natural Science Foundation of China     

Micro/macromechanical plastic limit analyses of composite materials and structures

The load-bearing capacities of ductile composite materials and structures are studied by means of a combined micro/macromechanics approach. Firstly, on the microscopic scale, the aim is to get the macroscopic strength domains by means of the homogenization theory of micromechanics. A representative volume element (RVE) is selected to reflect the microstructures of the composite materials. By introducing the homogenization theory into the kinematic limit theorem of plastic limit analysis, an optimization format to directly calculate the limit loads of the RVE is obtained. And the macroscopic yield criterion can be determined according to the relation between macroscopic and microscopic fields. Secondly, on the macroscopic scale, by introducing the Hill's yield criterion into the kinematic limit theorem, the limit loads of orthotropic structures such as unidirectional fiber-reinforced composite structures are worked out. The finite element modeling of the kinematic limit analysis is deduced into a nonlinear mathematical programming with equality-constraint conditions that can be solved by means of a direct iterative algorithm. Finally, some examples are illustrated to show the application of the present approach. Project supported by the National Natural Science Foundation of China (No. 19902007), the National Foundation for Excellent Doctoral Dissertation of China (No. 200025), the Fund of the Ministry of Education of China for Returned Oversea Scholars and the Basic Research Foundation of Tsinghua University.     

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

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

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

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

An analysis of the modeling of stress and strain fields in a real microstructure using a hybrid method

The distribution of stress and strain fields in a micro-structural area of a particle reinforced composite is studied by a combination of experimental and numerical method (hybrid method). With the experimental values of displacements in a micro-region as the boundary loading condition, strain and stress fields inside the micro-region are calculated by the finite element method under two different kinds of modeling, namely, as plane stress and plane strain condition. The differences between the two kinds of modeling conditions as applied to micro-structural areas are discussed. Project Supported by the National Natural Science Foundation of China (19972046) and National Overseas Study Foundation.     

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

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

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

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

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.     

SIMPLE SHAKEDOWN OF STRUCTURES UNDER VARIABLE MULTI-LOADINGS

The shakedown analysis of structures under variable multi-loadings is considered, and the corresponding simple shakedown condition is presented in this paper. Distribution of fixed stresses field is given, and the self-equilibrium of fixed stresses field is analyzed. Elastic shakedown and plastic shakedown conditions are presented based on the fixed stresses field. The theorem is convenient to evaluate the shakedown limit of structures under cyclical variable multiloadings through solving positive scalar fields and fixed stresses field factors at a series of dangerous positions of the structure, and tedious computations are avoided. Finally the theorem is applied to a thick-walled cylindrical tube under variable pressure and temperature, and the rolling contact problem. The results are in good agreement with some computational results.     

Dynamic shakedown theory of elastoplastic work-hardening structures allowing for second-order geometric effects

This paper is based on piecewise linear yield surface and discretization of structure. By allowing for inertial force, damping force and second-order geometric effects, the two generalized dynamic shakedown theorems are given for shakedown analysis of structure.The project is supported by National Natural Science Foundation of China and Jilin Provincial Applying-Basic Research Projects.     

Application of the quadrilateral area coordinate method： a new element for laminated composite plate bending problems

Recently, some new quadrilateral finite elements were successfully developed by the Quadrilateral Area Coordinate (QAC) method. Compared with those traditional models using isoparametric coordinates, these new models are less sensitive to mesh distortion. In this paper, a new displacement-based, 4-node 20-DOF (5-DOF per node) quadrilateral bending element based on the first-order shear deformation theory for analysis of arbitrary laminated composite plates is presented. Its bending part is based on the element AC-MQ4, a recent-developed high-performance Mindlin-Reissner plate element formulated by QAC method and the generalized conforming condition method; and its in-plane displacement fields are interpolated by bilinear shape functions in isoparametric coordinates. Furthermore, the hybrid post-processing procedure, which was firstly proposed by the authors, is employed again to improve the stress solutions, especially for the transverse shear stresses. The resulting element, denoted as AC-MQ4-LC, exhibits excellent performance in all linear static and dynamic numerical examples. It demonstrates again that the QAC method, the generalized conforming condition method, and the hybrid post-processing procedure are efficient tools for developing simple, effective and reliable finite element models. The project is supported by the National Natural Science Foundation of China (10502028), the Special Foundation for the Authors of the Nationwide (China) Excellent Doctoral Dissertation (200242), and the Science Research Foundation of China Agricultural University (2004016).     

A study on stress intensity factors and singular stress fields of three-dimensional interface crack

By using the finite-part integral concepts and limit technique, the hypersingular integrodifferential equations of three-dimensional (3D) planar interface crack were obtained; then the dominant-part analysis of 2D hypersingular integral was further used to investigate the stress fields near the crack front theoretically, and the accurate formulae were obtained for the singular stress fields and the complex stress intensity factors. After that, a numerical method is proposed to solve the hypersingular integrodifferential equations of 3D planar interface crack, and the problem of elliptical planar crack is then considered to show the application of the method. The numerical results obtained are satisfactory. Project supported by the Foundation of Solid Mechanics Open Research Laboratory of State Education Commission at Tongji University and the National Natural Science Foundation.     

Boundary element analysis for elastic and elastoplastic problems of 2D orthotropic media with stress concentration

Both the orthotropy and the stress concentration are common issues in modern structural engineering. This paper introduces the boundary element method (BEM) into the elastic and elastoplastic analyses for 2D orthotropic media with stress concentration. The discretized boundary element formulations are established, and the stress formulae as well as the fundamental solutions are derived in matrix notations. The numerical procedures are proposed to analyze both elastic and elastoplastic problems of 2D orthotropic media with stress concentration. To obtain more precise stress values with fewer elements, the quadratic isoparametric element formulation is adopted in the boundary discretization and numerical procedures. Numerical examples show that there are significant stress concentrations and different elastoplastic behaviors in some orthotropic media, and some of the computational results are compared with other solutions. Good agreements are also observed, which demonstrates the efficiency and reliability of the present BEM in the stress concentration analysis for orthotropic media. The project supported by the Basic Research Foundation of Tsinghua University, the National Foundation for Excellent Doctoral Thesis (200025) and the National Natural Science Foundation of China (19902007). The English text was polished by Keren Wang.     

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.     

Analysis of interlaminar stress and nonlinear dynamic response for composite laminated plates with interfacial damage

By considering the effect of interfacial damage and using the variation principle, three-dimensional nonlinear dynamic governing equations of the laminated plates with interfacial damage are derived based on the general sixdegrees-of-freedom plate theory towards the accurate stress analysis. The solutions of interlaminar stress and nonlinear dynamic response for a simply supported laminated plate with interfacial damage are obtained by using the finite difference method, and the results are validated by comparison with the solution of nonlinear finite element method. In numerical calculations, the effects of interfacial damage on the stress in the interface and the nonlinear dynamic response of laminated plates are discussed.     

Triangular grid method for stress-wave propagation in 2-D orthotropic materials

A triangular grid method is presented to calculate propagation problems of elastic stress waves in 2-D orthotropic materials. This method is based on the dynamic equilibrium equations of the computational cells formed among the auxiliary triangular grids. The solution is obtained by calculating alternately the nodal displacements and the central point stresses of the spatial grids. The numerical results are compared with the corresponding solutions of the finite element method. Comparisons show that the triangular grid method yields a higher calculational speed than the finite element method. The stress concentrations are investigated from wave-field analyses when the stress wave propagates within an orthotropic plate with a hole. Finally, the presented numerical method is used to study the features of wave propagation and diffraction in a square orthotropic plate with a hole when an impact load is applied to the top of the plate.This work was supported by National Natural Science Foundation of China (Nos. 10025212 and 10232040) and Natural Science Foundation of Liaoning province (No. 20021070).