线性等参单元逆变换的几何二分法及其修正

在证明线性等参单元的一个几何特性的基础上,提出了线性等参单元逆变换的几何二分法。几何二分法简单易行,避免求解高维高次非线性方程组。在几何二分法基础上,提出一种修正措施以消除畸形单元对所求结果的影响,几何二分法及其修正措施克服了以往方法存在的缺点,两个实例验证了由几何二分法及其修正措施进行等参逆变换的正确性。     

一种基于面-面局部搜索的接触算法

接触-碰撞算法是影响动态接触问题数值分析精度的重要因素,发展健壮、精确的局部搜索算法对提高数值模拟精度具有重要意义。为解决点-面接触搜索算法存在的盲区问题,提出了一种基于面-面的局部搜索算法,并发展了相应的接触力计算方法。新的接触算法以面心坐标与特征长度表征接触片,进行预搜索快速排除不会发生接触的潜在接触对;利用主从接触片的投影与侵入深度关系确定接触片间的接触状态;利用等参逆变换解析给出接触点,避免了方程组的迭代求解。典型算例与工程实例计算表明,本文算法消除了接触搜索的盲区,具有很好的健壮性与计算精度。     

一种基于面-面局部搜索的接触算法

接触-碰撞算法是影响动态接触问题数值分析精度的重要因素,发展健壮、精确的局部搜索算法对提高数值模拟精度具有重要意义。为解决点-面接触搜索算法存在的盲区问题,提出了一种基于面-面的局部搜索算法,并发展了相应的接触力计算方法。新的接触算法以面心坐标与特征长度表征接触片,进行预搜索快速排除不会发生接触的潜在接触对;利用主从接触片的投影与侵入深度关系确定接触片间的接触状态;利用等参逆变换解析给出接触点,避免了方程组的迭代求解。典型算例与工程实例计算表明,本文算法消除了接触搜索的盲区,具有很好的健壮性与计算精度。     

等参元逆变换算法在渗流—位移耦合场分析中的应用

在渗流 -位移耦合场分析中应用等参元逆变换算法解决了两个场之间的数据传递问题 ,所提出的逆变换算法是基于解析性质的新算法 ,具有简便、高效和高精度的特点 ,便于在实际工程中应用。     

非协调四结点平面等参位移元新列式方法

本文导出一种等参协调元位移函数的新的表示方法，在此基础上建立起了构造等参非协调元的新方法，作为实例，构造出两个可以给出单刚显式的四结点平面非协调新单元。     

等参管单元及其在热传导问题边界元法中的应用Isoparametric tube elements and their application in heat conduction BEM analysis

采用边界元法（BEM ）求解实际工程问题时,很大一部分误差来自于离散误差。为此,本文基于Lagrange插值原理,提出了一种三维等参管单元边界元算法,该单元能很好地模拟管状结构的几何外形并对物理量进行高阶插值,大大地消除了离散误差。另外,当在边界元法中使用等参管单元时,提出了一种在等参平面内消除积分奇异性的方法。算例表明,本文算法具有划分网格少,求解精度高的优点。     

等参管单元及其在热传导问题边界元法中的应用

采用边界元法(BEM)求解实际工程问题时,很大一部分误差来自于离散误差。为此,本文基于Lagrange插值原理,提出了一种三维等参管单元边界元算法,该单元能很好地模拟管状结构的几何外形并对物理量进行高阶插值,大大地消除了离散误差。另外,当在边界元法中使用等参管单元时,提出了一种在等参平面内消除积分奇异性的方法。算例表明,本文算法具有划分网格少,求解精度高的优点。     

板问题混合变量等参Hamiltonian元的半解析解

本文给出板问题混合变量方程的一种半离散半解析方法－混合变量等参Ｈａｍｉｌｔｏｎｉａｎ元。该方法沿板厚方向未作任何有关位移和应力的人为假设，而是采用控制论中方法给出真解，所以可以求解任意厚度板问题。     

非协调四结点平面等参位移元新列式方法

本文导出一种等参协调元位移函数的新的表示方法,在此基础上建立起了构造等参非协调元的新方法。作为实例,构造出两个可以给出单刚显式的四结点平面非协调新单元。     

基于参变量变分原理的非均质多边形微结构材料弹塑性分析

为了更好地模拟复合材料及含夹杂非均质材料等的宏观弹塑性力学性能,简化有限元建模时间和减少有限元模拟计算量。本文基于参变量变分原理,提出了一种采用任意多边形弹塑性单元进行结构非线性分析的参数二次规划算法,给出了参变量最小势能原理以及最终的二次规划模型,并在有限元分析与优化设计软件系统JIFEX上进行了程序实现。数值算例证明了本文方法的正确与可行性。     

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

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

超椭球模型下结构非概率可靠性指标的迭代算法

迭代算法对于非概率可靠性指标的求解及其优化问题具有重要意义。本文基于不确定参数的超椭球描述,研究求解非概率可靠性指标的迭代算法。针对极限状态方程非线性情况较高时可能存在不收敛的问题,提出一个检测严重迂回振荡的判据,并在HL-RF迭代公式的基础上引入修正解,在一定程度上克服迭代不收敛的问题。数值算例验证了迭代算法的正确性和有效性。     

A spectral collocation method for compressible,non-similar boundary layers

An efficient and highly accurate algorithm based on a spectral collocation method is developed for numerical solution of the compressible, two-dimensional and axisymmetric boundary layer equations. The numerical method incorporates a fifth-order, fully implicit marching scheme in the streamwise (timelike) dimension and a spectral collocation method based on Chebyshev polynomial expansions in the wall-normal (spacelike) dimension. The discrete governing equations are cast in residual form and the residuals are minimized at each marching step by a preconditioned Richardson iteration scheme which fully couples energy, momentum and continuity equations. Preconditioning on the basis of the finite difference analogues of the governing equations results in a computationally efficient iteration with acceptable convergence properties. A practical application of the algorithm arises in the area of compressible linear stability theory, in the investigation of the effects of transverse curvature on the stability of flows over axisymmetric bodies. The spectral collocation algorithm is used to derive the non-similar mean velocity and temperature profiles in the boundary layer of a ‘fuselage’ (cylinder) in a high-speed (Mach 5) flow parallel to its axis. The stability of the flow is shown to be sensitive to the gradual streamwise evolution of the mean flow and it is concluded that the effects of transverse curvature on stability should not be ignored routinely.     

优化迭代步长的两种改进增量谐波平衡法

增量谐波平衡法(IHB法)是一个半解析半数值的方法, 其最大优点是适合于强非线性系统振动的高精度求解. 然而, IHB法与其他数值方法一样, 也存在如何选择初值的问题, 如初值选择不当, 会存在不收敛的情况. 针对这一问题, 本文提出了两种基于优化算法的IHB法: 一是结合回溯线搜索优化算法(BLS)的改进IHB法(GIHB1), 用来调节IHB法的迭代步长, 使得步长逐渐减小满足收敛条件; 二是引入狗腿算法的思想并结合BLS算法的改进IHB法(GIHB2), 在牛顿-拉弗森(Newton-Raphson)迭代中引入负梯度方向, 并在狗腿算法中引入2个参数来调节BSL搜索方式用于调节迭代的方式, 使迭代方向沿着较快的下降方向, 从而减少迭代的步数, 提升收敛的速度. 最后, 给出的两个算例表明两种改进IHB法在解决初值问题上的有效性.      

改进的速度迭代弹塑性大变形动力半显式算法

提出了基于速度迭代的弹塑性大变形动力半显式算法，并且在此基础上采用拟弯曲膜单元，结合有模卸载方式数值模拟了U形弯曲以及S形导轨回弹过程，着重比较了动力半显式算法与动力显式算法和实验的回弹结果．给出了动力半显式算法的迭代格式和收敛标准．     

Bruck formula for a perturbed Lipschitzian iteration of Lipschitz pseudocontractive maps

Introduction InRef.[1],BrowderdefinedpseudocontractivemappinginBanachspaceand mentionedtherelationshipofthisclassofmappingtoanimportantclassofmappingknown asaccretive(monotoneinHilbertspace)operator.ThemappingAisaccretiveiffI-Ais pseudocontractive.Themapp…     

The back and forth nudging algorithm applied to a shallow water model,comparison and hybridization with the 4D‐VAR

We study in this paper a new data assimilation algorithm, called the back and forth nudging (BFN). This scheme has been very recently introduced for simplicity reasons, as it does not require any linearization, or adjoint equation, or minimization process in comparison with variational schemes, but nevertheless it provides a new estimation of the initial condition at each iteration. We study its convergence properties as well as efficiency on a 2D shallow water model. All along the numerical experiments, comparisons with the standard variational algorithm (called 4D‐VAR) are performed. Finally, a hybrid method is introduced, by considering a few iterations of the BFN algorithm as a preprocessing tool for the 4D‐VAR algorithm. We show that the BFN algorithm is extremely powerful in the very first iterations and also that the hybrid method can both improve notably the quality of the identified initial condition by the 4D‐VAR scheme and reduce the number of iterations needed to achieve convergence. Copyright © 2008 John Wiley & Sons, Ltd.     

影响结构优化理性准则法优化效果的关键问题研究

深入仔细分析决定理性准则法优化效果与优化效率的关键问题——准则是否准确和迭代计算是否收敛。首次提出：在结构优化中,重量作为设计资源除有改善结构性能的一面之外,还会有作为载荷导致结构性能劣化的矛盾的另一面。虚功准则方程组忽略载荷导数,无法考虑重量矛盾作用的另一面,这种忽略不能视为合理近似。因为对于航空航天器、高精度天线和高速运行的机械、车辆等以自重和惯性载荷为主的一大类工程结构,这种忽略导致虚功法,得到的解离最优解相差甚远。导重准则法是严密推导的理性准则法,克服了虚功法准则不准的缺陷,优化效果大幅度提高。理性准则法最后归结为非线性准则方程组的迭代求解,由于工程结构优化的准则方程组难以满足其严格的收敛条件,而可采用步长因子法求解,可以证明使迭代收敛的步长因子一定存在,并可给出步长因子理论取值范围和实际取值方法。以十杆桁架考题和两个天线结构优化为例验证了以上论点。     

Convergence theorem of common fixed points for Lipschitzian pseudo-contraction semi-groups in Banach spaces

The purpose of this paper is to study the weak convergence problems of the implicity iteration process for Lipschitzian pseudocontractive semi-groups in the general Banach spaces. The results presented in this paper extend and improve the corresponding results of some people.     

Optimization of the nonlinear convergence rate of the SIMPLER Picard algorithm for steady incompressible internal flows through enforcement of continuity projections

We discuss the performance optimization of the Semi‐Implicit Method for Pressure‐Linked Equations—Revised (SIMPLER) Picard algorithm for steady incompressible internal flows. We discuss the nonlinear convergence of the Picard iteration as a function of the pressure and scalar potential continuity projections stemming from the SIMPLER algorithm, for three example problems. In particular, we discuss the choice of under‐relaxation method, and choice of under‐relaxation factors; the choice of the projection algorithm; and the required tolerance for the linear solve of the generalized Poisson equations for the pressure and scalar potential equations that arise from the projection operations. We conclude that the convergence of the nonlinear Picard iteration can be effectively controlled by the optimal enforcement of the continuity projections. Copyright © 2005 John Wiley & Sons, Ltd.