Upper-bound limit analysis based on the natural element method

The natural element method (NEM) is a newly-developed numerical method based on Voronoi diagram and Delaunay triangulation of scattered points, which adopts natural neighbour interpolation to construct trial functions in the framework of Galerkin method. Owing to its distinctive advantages, the NEM is used widely in many problems of computational mechanics. Utilizing the NEM, this paper deals with numerical limit analysis of structures made up of perfectly rigid-plastic material. According to kinematic theorem of plastic limit analysis, a mathematical programming natural element formulation is established for determining the upper bound multiplier of plane problems, and a direct iteration algorithm is proposed accordingly to solve it. In this algorithm, the plastic incompressibility condition is handled by two different treatments, and the nonlinearity and nonsmoothness of the goal function are overcome by distinguishing the rigid zones from the plastic zones at each iteration. The procedure implementation of iterative process is quite simple and effective because each iteration is equivalent to solving an associated elastic problem. The obtained limit load multiplier is proved to monotonically converge to the upper bound of true solution. Several benchmark examples are investigated to validate the significant performance of the NEM in the application field of limit analysis.     

弯管结构塑性极限分析的数值方法及应用

从塑性极限分析数值计算的角度,分析了多组载荷联合作用下弯管结构的塑性极限承载能力。为了克服塑性极限上限分析中目标函数非线性非光滑所导致的数值困难,提出了一种弯管结构塑性极限上限分析的无搜索优化迭代算法;采用一种改进的弯管单元并利用加载路径的径向射求解方案处理多组载荷系统。通过对典型弯管结构进行塑性极限分析得出了一些有价值的结论。     

基于噪声和残缺投影图像的重建算法研究

滤波反投影法和迭代法是SR-CT重建算法中的两大类算法,松弛迭代法是迭代法中的一种;本文主要研究滤波反投影方法和松弛迭代法.为了研究两种方法的影响因素,讨论了滤波反投影法中滤波函数的选取,分析了松弛迭代法中重要参数,如松弛因子、迭代次数对重建效果的影响,给出了重建效果和这些参数的影响关系;对比了两种方法在处理噪声和残缺...     

Structure buckling hybrid reliability analysis of a supercavitating projectile using a model with truncated probability and multi-ellipsoid convex set uncertainties

Structure buckling problems for supercavitating projectiles are often observed in high underwater velocity operating conditions. As a result, it is necessary to perform a structure buckling reliability analysis. It is common that \hboxprobabilistic and nonprobabilistic uncertainty information exist simultaneously. Also, it is reasonable that probability distributions of most random variables in engineering are treated as truncated probability distributions. In this paper, a new hybrid model is proposed, which deals with structure performance function with both truncated probability distribution variables and multi-ellipsoid convex set variables. The model discussed here is adapted for two separate cases in the standard super-sphere space, i.e., limit state surface interferences with a unit super-sphere region or not. The hybrid reliability index is calculated using a modified limit step length iteration algorithm to ensure convergence. Good convergence and validity of the iteration algorithm are verified using numerical examples with highly nonlinear structure performance function. The hybrid model is applied to the structure buckling hybrid reliability analysis of a supercavitating projectile. Results show that structure buckling hybrid reliability index increases with the increase in the ratio of base diameter to cavitator diameter, and decreases with the increase of initial launch velocity. Also, uncertainty degree of cavitator drag coefficient and initial launch velocity should be controlled in demonstration for high structure buckling reliability of supercavitating projectiles.     

结构可靠度FORM方法的混沌动力学分析

引入混沌动力学理论讨论了FORM收敛失败的非线性动力学根源. 给出了几个典 型函数在参数区间上的可靠指标分岔图，展示了极限状态函数经过FORM迭代成为 非线性映射后计算结果的周期振荡、分岔和混沌等复杂动力学现象，计算了非线性映射的 Lyapunov指数. 结果表明，极限状态函数设计点的曲率大小与FORM的收敛性没有简单的联 系，判别FORM迭代计算收敛性的指标是非线性映射的Lyapunov指数.     

Computational limit analysis of rigid-plastic bodies in plane strain

This paper deals with the computational methods for limit analysis of plane strain problems. The finite element mathematical programming formula (FE-MPF) for determining the upper bound load multiplier established by the authors earlier is adopted and modified for plane strain problems. The penalty method is used to impose the incompressibility constraint. The FE-MPF is solved by a direct iteration procedure without the need of a searching process. This algorithm is not sensitive to the volumetric locking effect. And it can be easily extended to the limit analysis of three dimensional problems. The results of numerical examples are satisfactory and show the stable convergency of the present algorithm.     

多变量单元结构共振解耦动力优化设计方法

针对截面多变量单元结构的动力优化问题,建立了带频率约束的结构动力优化设计模型,对隐式非线性频率约束函数进行Taylor近似展开,给出了截面多变量单元频率梯度函数显式表达式,基于Kuhn-Tucker条件构建迭代算法,其中拉格朗日乘子通过建立联合方程组求解,形成了含多变量单元共振解耦优化设计方法。以矩形截面单元结构为算例的结果表明,所给单元多变量算法具有良好的准确性;截面变量对频率的贡献存在主次之分,区分指标可采用梯度值;次要变量的修正因子迭代中可采用定值,且其下限应尽可能降低,利于节约成本。本文工作对多变量复杂截面结构动力优化设计可提供理论指导,提高结构动力优化方法的适用性。     

接触分析的光滑模型及迭代算法

利用变分不等式和基于信息熵的凝聚函数把有摩擦接触问题模型化为一个标准的凸二次规划问题,极大地简化了这一复杂的问题,同时引入摩擦方向约束并构造了以无摩擦解为初值的迭代算法,在较摩擦系数时计算也能保证收敛,算例表明算法高效可靠。     

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.     

Non-linear analysis of structures using two-point method

Non-linear algebraic equations must be solved by an iterative method, the non-linear equations being linearized by evaluating the non-linear terms with the known solution from the preceding iteration. The Newton-Raphson method, which is based on the Taylor series expansion and uses the tangent stiffness matrix, has been extensively used to solve non-linear problems. In this paper, a new Newton-Raphson algorithm is developed for analyses involving non-linear behavior. Our method, here named as a two-point method, is constructed as a predictor-corrector one, most frequently taking Newton's method in the first iteration. It should be noted that our concern in this research ignores the problem of passing limit points. The presented method incorporates the known information at each stage of the loading process to determine the subsequent unknown variables. Compared with the classic Newton-Raphson algorithm, it offers a strategy that can be deployed to reduce both the number of the iterations and the computing time involved in non-linear analysis of structures.     

一种求算结构可靠度指标的新方法

为解决结构功能函数在结构设计点附近非线性程度较高时,一次可靠度计算方法(如JC法)不收敛的问题,提出一种新的可靠度计算方法.该算法根据结构可靠指标的几何意义,先将迭代点靠近极限状态面,再通过在极限状态面上搜索下一迭代点,逐渐减小极限状态面上的点到标准正态空间中坐标原点的距离,从而达到收敛的目的.与其他方法相比,该方法不...     

捷联惯导姿态更新的双子样二次迭代优化算法研究

推导了捷联惯导姿态更新算法中双子样迭代算法，提出以锥误差最小为优化标准的双子样迭代优化算法，并对双子样，三子样，双子样迭代算法在锥运动条件下的情况做了仿真分析。     

离散变量结构形状优化设计的综合算法

提出一种离散变量结构形状优化设计的方法,该方法在优化过程中将两类变量统一考虑,根据两类变量的相对差商值确定搜索的变量;并且对形状变量施加了运动极限的限制,从而保证了收敛速度快,收敛平稳;将该算法应用于几个经典的结构优化算例：运算结果显示了该方法计算效率高,优化效果也比较满意．     

Geometrical nonlinear analysis of structures using residual variables

The capability of nonlinear analysis methods in tracing the equilibrium path depends on how to return to static status. In this paper, some residual factors existed in the iteration steps are employed. By setting the residual load parameter to zero, minimizing its factor and the residual displacement parameter, three novel constraint equations are obtained. The new formulas are proved by numerical examples. All obtained results illustrate the minimum residual load scheme's robustness in comparison with the cylindrical arc-length algorithm and other previous strategies. Moreover, the capacities of new procedures in tracing the load and displacement limit points are assessed.     

应用响应面方法结合空间映射方法的可靠性评价

应用响应面结合空间映射方法,在第1次迭代拟合极限状态函数,其它迭代 应用映射技术在第1迭代响应面基础上映射调整得到新的极限状态函数,并进行可靠性分析. 这样就改变了序列响应面方法评价可靠性时需要反复对模型进行试验设计、分析并拟合极限状 态函数的执行过程,从而大大降低了模型分析的计算量.     

基于单纯形寻优的响应面可靠性分析方法

结构可靠度计算常采用经典的响应面法拟合隐式功能函数或高维功能函数,而对于强非线性功能函数的实际工程问题,尽管其能够计算出结构可靠度的结果,但此时多项式响应面的拟合精度不够,很容易造成不收敛的现象。为了解决上述问题,将响应面法与单纯形寻优的思路进行结合来探求一种有效的计算方法。本文利用单纯形算法对每次迭代的验算点进行优化;再以优化后的设计验算点为中心进行取样,利用响应面法循环迭代计算;最后,沿着真实响应面逐渐逼近最终的验算点。该方法能够解决高维非线性的隐式极限状态方程可靠度计算收敛性的问题,可以提高计算精度和计算效率,具有一定的工程适用性。     

结构可靠指标求解的一种新的迭代方法

由国际安全度联合委员会推荐的一次二阶矩可靠度方法是国际工程界公认的一种较好的可靠度分析方法（简称ＪＣ法）。这种方法具有迭代格式简单，收敛快的特点。但许多实际计算表明，这种方法只适用于功能函数非线性程度不高的情况，否则不收敛。本文提出一种新的迭代方法，这种方法迭代格式与ＪＣ法同样简单，一般情况下均能保证收敛．     

扁壳在任意法向载荷下的极限分析

本文利用对偶变分原理,从极限分析上限定理导得扁壳在任意法向载荷下的极限计算变分公式,通过把扁壳划分为有限元和叠代的方法获得极限载荷的近似值,通过算例显示了这种新的变分方法是解决扁壳极限分析的有效方法,可以用来对飞行器、船舶、防护罩等结构进行极限分析,具有一定现实意义。     

基于频率和模态保证准则的结构内部多缺陷反演

静响应(位移、应变等)在实际问题的反演分析中很难由安装在结构上的一组传感器记录得到,而结构的动力特性(频率、振型)和 动力响应(加速度、速度、动位移)在实际问题中较易通过传感器采集得到. 文中基于频率残差和模态保证准则构建了反演分析模型的目标函数,并结合频域内动力扩展有限元法和人工蜂群智能优化算法 的优点,扩展有限元法通过引入非连续位移模式在不重新划分网格的情况下通过改变水平集函数反映缺陷的数量、位置及大小, 避免了反演分析每次迭代过程中的网格重剖分,人工蜂群智能优化算法在每次迭代中都采用全局和局部搜索,找到最优解的概 率大幅增大并可很好地避免局部最优,同时,通过引入拓扑变量,将缺陷的数量纳入到反演分析过程中,迭代过程中可智能反演出缺陷的数目,建立了结构内部多缺陷(孔洞、裂纹)的反演分析模型. 通过若干算例的分析表明:建立的反演分析模型能够较为准确地探测出结构内部圆形、椭圆形以及裂纹状缺陷的数量、位置及大小,且算法具有较好的鲁棒性.      

An efficient iteration algorithm of mixed boundary/finite element method and application to free torsional vibration analysis of bodies of revolution

In this paper,the general formulation of a new proposed iteration algorithm of mixedBEM/FEM for eigenvalue problems of elastodynamics is described.Approximatefundamental solutions of elastodynamics are adopted in the normal mixed BEM/FEMequations.The accuracy of solutions is progressively improved by the iteration procedure.Not only could the awkwar dness of non-algebraic eigenvalue equations be avoided but alsothe accuracy of numerical solutions is almost independent of the interior meshing.All thesegive many advantages in numerical calculation.The algorithm is applied to free torsionalvibration analysis of bodies of revolution.A few cases are studied.All of the numericalresults are very good.