极限分析的无搜索数学规划算法

本文研究理想刚塑性介质极限载荷因子的计算方法。根据极限分权理论的上限定理,建立了计算极限载荷因子的一般数学规划有限元格式。针对这种格式的特点,提出了一个求解极限载荷因子的无搜索迭代算法。这个算法中采用逐步识别刚性、塑性分区,不断修正目标函数的方案,克服了目标函数非光滑所导致的困难。本文提出的算法建立于位移模式有限元基础上,有较广的适用范围,且具有计算效率高,稳定性好,格式简单易于程序实现等优点。     

含缺陷结构的塑性极限分析

结合极限分析中的数学规划理论和有限元技术,提出了三维含缺陷结构极限分析的数学规划方法,并采用罚函数法引入塑性不可压条件.对于考虑多组独立变化载荷联合作用的情况,提出了加载路径射线辐射求解方案,并基于这种射线辐射状的加载路径,推导了多组载荷联合作用下结构塑性极限上限分析的数崐学规划格式,编制了相应的有限元程序.文中的数值结果表明了该方法的正确性与有效性.     

吊耳的极限分析

将吊耳的外形及耳环的横截面作了简化,假定吊钩对耳环的横向压力均匀分布。选择了所有可能的机动模型,推导了相应的上限载荷公式,最后作了数值计算。     

极限分析的弹性补偿法及其应用

有限元法与数学规划法相结合，应用极限上、下限定理，将极限分析归结为求解最优化问题，是目前被普遍应用的极限分析方法，但是该方法受到计算能力的限制，难以应用到实际工程问题中。鉴于此，本文介绍一种基于线弹性分析基础上的简单的求解复杂结构极限栽荷下限、上限的方法——弹性补偿法，同时结合三维有限元分析，求解内压下三通结构的极限载荷。通过与弹塑性分析结果比较发现，简单的弹性补偿法能够很好的评估复杂三雏结构的塑性承载能力。     

极限分析和安定分析的并行算法

本文基于结构极限分析和安定分析的温度参数法，研究了线性规划的并行算法，温度参数法求解的关键是线性规划的求解效率，在分布式布储的多处理机系统上，采用了线性规划的分块单纯形法。对相当于大量工况的结构，子结构是一个很好的并行结构分析方法。     

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

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

运用有限元方法确定T型管节点的极限载荷

本文应用ADINA程序对支管轴向受拉和轴向受压的T型管节点模型实例进行了非线性有限元计算.在轴向受压的工况计算中,除考虑材料的物理非线性外,还考虑了几何非线性问题.根据对计算结果的分析,确定了两种工况的极限载荷,以及支管和主管在不同载荷步下的塑性变形区的扩展情况.     

孔洞受压的统一极限分析

采用俞茂宏统一强度理论分别就平面应变和平面应力总是对受压孔洞进行了极限分析，得到了各自的统一解形式。以往的基于Ｔｒｅｓｃａ、Ｍｉｓｅｓ、双剪屈服准则以及Ｍｏｈｒ－Ｃｏｕｌｏｍｂ强度理论的极限解均为本文统一解的特例。此解可以适应于广泛的不同性质的各灯材料，合理也得出不同材料的相应解。     

理想塑性结构极限与安定分析的数值方法

极限分析和安全分析的近代发展方向是寻找通用性强，计算效率高的数值方法。本文介绍将有限单元法和数学规划法相结合的、同时适用于极限分析和安全分析的统一数值方法，包括下限格式和上限格式。     

摩擦接触弹塑性分析的数学规划法

本文提出摩擦接触弹塑性分析的数学规划法,是弹性接触问题解法的推广,使之能处理比例加载时的塑性效应。利用所导出接触问题总刚度阵的逆,结合参变量线性互补公式,仍可消去由接触单元而引入的惩罚因子。于是,可以获得一个弹塑性问题、弹性接触和弹塑性接触问题的有限元参数二次规划法模型。所选的几个实例说明了方法的实用性。     

边坡三维极限平衡法的通用形式

基于以下假定条件:(1) 稳定系数定义为材料的强度折减系数;(2) 土体为刚体,底滑面服从Mohr-Columb强度破坏准则;(3) 微条柱底部法向力dNz的作用点处于条柱底部中点;(4)滑面剪力与底滑面和xoz平面交线的夹角为θ。本文建立了边坡三维极限平衡法的通用形式,通过给定不同的限制条件,可分别得到三维普通条分法 、三维简化毕肖普法 、三维简化简布法 、三维Spencer法 等三维极限平衡的具体算法。     

平面滑动型岩质边坡稳定性极限分析上限法

在极限分析理论框架体系的基础上,提出了平面滑动型岩质边坡极限分析上限法。该方法按照外力所做的功率等于内力所消耗的功率,即滑体处于极限状态时两功率相等的条件——虚功率方程,综合考虑作用在岩质边坡上的后缘裂缝静水压力、沿滑面扬压力、重力、水平地震作用力、锚固力等外力,按照强度折减法,推导得出了岩质边坡稳定性评价的极限分析上限解。锦屏一级水电站右岸泄洪洞引渠内侧工程边坡属于典型的多级台阶状岩质高边坡,该边坡岩体结构为层状结构,变形破坏模式为平面滑动。对引渠内侧岩质边坡稳定性分析结果表明,极限分析上限法对平面滑动型岩质边坡有较好的适用性。     

On the numerical convergence and performance of different spatial discretization techniques for transient elastodynamic wave propagation problems

We present a systematic investigation of several discretization approaches for transient elastodynamic wave propagation problems. This comparison includes a Finite Difference, a Finite Volume, a Finite Element, a Spectral Element and the Scaled Boundary Finite Element Method. Numerical examples are given for simple geometries with normalized parameters, for heterogeneous materials as well as for structures with arbitrarily shaped material interfaces. General conclusions regarding the accuracy of the methods are presented. Based on the essential numerical examples an expansion of the results to a wide range of problems and thus to numerous fields of application is possible.     

A transient finite element analysis of natural convection around a horizontal hot cylinder

This paper presents an advanced method for a 2-dimensional analysis of transient natural convection by finite element method. The present method, based on stream function—vorticity formulation, could get rid of numerical errors and constraint of perpendicular mesh subdivision, since we excluded a finite difference approximation of vorticity on no-slip boundaries. A considerable effect of upwind weighting function was examined. The method was successfully applied to a problem of natural convection around a horizontal hot cylinder.     

Variable penalty method for finite element analysis of incompressible flow

A new scheme is applied for increasing the accuracy of the penalty finite element method for incompressible flow by systematically varying from element to element the sign and magnitude of the penalty parameter λ, which enters through ?.v + p/λ = 0, an approximation to the incompressibility constraint. Not only is the error in this approximation reduced beyond that achievable with a constant λ, but also digital truncation error is lowered when it is aggravated by large variations in element size, a critical problem when the discretization must resolve thin boundary layers. The magnitude of the penalty parameter can be chosen smaller than when λ is constant, which also reduces digital truncation error; hence a shorter word-length computer is more likely to succeed. Error estimates of the method are reviewed. Boundary conditions which circumvent the hazards of aphysical pressure modes are catalogued for the finite element basis set chosen here. In order to compare performance, the variable penalty method is pitted against the conventional penalty method with constant λ in several Stokes flow case studies.     

A note on the stability and accuracy of Cn finite elements in steady diffusion-convection problems

The paper is concerned with stability and accuracy of an nth order Lagrangian family of finite element steady-state solutions of the diffusion-convection equation, and furthermore is concerned with the stability and the accuracy of on mth kind Hermitian family of finite element solutions. We discuss the stability of the numerical solution based on the fact that the characteristic finite element solution can be expressed approximately as a rational function of cell Peclet number Pe_c ( = uh/k). Moreover, it is shown that by eliminating derivatives and by using the interpolation method over elements a stable solution is obtained over the domain independent of Pe_c for P^1,3, and for P^2,5 the stable solution is obtained for Pe_c less than 44.4.     

A brief note on upwind collocation

Upwind collocation on Hermite cubics is compared to orthogonal collocation with respect to effective diffusion. The one-dimensional constant coefficient advection-diffusion equation is employed to this end. The effective diffusion coefficient is determined exactly and is found to be dependent on the nodal solution values. The effective diffusion coefficients of three other upwinding schemes are also presented. Upwind collocation is found to have an effective diffusion coefficient like other upwinding schemes plus an extra term which may enhance or reduce the non-advective flux, depending on the problem solution and point location within the domain.