一种基于Hamel形式的无条件稳定动力学积分算法

时间积分算法是求解动力学系统的一个核心问题.动力学方程的时间积分经常会出现数值不稳定现象,有限元空间离散也通常会造成伪高频振荡,因而,发展解决上述问题的数值积分算法具有重要的理论价值.本文基于Hamel场变分积分子,通过新的数值积分算法的构造方法,提出了一种无条件稳定的Hamel广义α方法,具体内容包括:构造特殊的变分形式,利用变分积分子等工具,建立无条件稳定的数值积分算法;在相同框架下,提出更高精度的数值格式;结合活动标架法的特性,将算法的一般形式推广到李群空间,得到Hamel广义α文所提出的Hamel广义α方法是无条件稳定的,具有二阶精度并且能够快速过滤掉虚假的高频振荡.数值算例的结果显示,本文所提方法具备了传统方法的精度、耗散和稳定性优势,既适合一般的线性空间,也适用于李群空间,同时还可以发展高阶精度算法.本文发展了构造变分积分子的新模式.     

一类带高频耗散的动力学显式算法

带数值耗散的算法因其能有效过滤虚假高频响应的影响而备受关注。基于离散控制理论Z变换,提出一类带数值耗散的结构动力学显式算法。该算法采用CR法的速度和位移递推式,满足零振幅衰减率,且对线性系统和非线性刚度软化系统为无条件稳定,对刚度硬化系统则是条件稳定的;该算法由单个参数ρ控制数值耗散能力。通过对振幅衰减和周期延长的理论分析表明,系数a可调节算法的精度和非线性稳定区间,给出精度最优时系数的取值。对特定的系数取值该算法可转变为CR法。通过算例对线性系统和非线性系统的分析验证了新算法具有良好的精度、稳定性和数值耗散,表明新算法是正确有效的。     

结构动力响应数值算法耗散和超调特性设计

对结构动力响应数值计算问题提出引入多个自由参数来获得所希望的算法特性. 多参数的一个明显的好处就是在算法设计上有更大的自由空间. 利用这些自由参数获得了两个新的无条件稳定、有二阶精度的、有好的耗散和没有超调的单步时间直接积分算法. 在存在阻尼情况下基于有限差分分析理论证明了新算法的这些特性. 其中一个有高频渐进消去特性,且在有阻尼情况下与Houbolt方法相比对高频有更强的耗散. 另一个在低频极限无耗散,高频耗散可以用一自由参数控制. 超调分析结果显示两个新算法都不显示超调,而HHT方法不仅有速度超调,还有位移超调. 最后使用一些算例并通过与传统方法的比较数值地验证了理论分析结果.     

黏塑性本构计算的稳定性分析

提出本构方程计算方法的稳定性问题,针对黏塑性本构计算的显式精确算法的稳定性进行分析,发现该算法并非无条件稳定,使用小扰动方法给出了其计算稳定的必要条件,稳定性条件对数值计算中的时间步长提出限制要求。通过有限元算例验证了分析的正确性,计算结果也表明理论推导得到的稳定性公式能够准确预测满足计算稳定性条件要求的最大时间步长与各参数之间关系。     

瞬态热传导问题的精细积分数值流形方法研究

数值流形方法（NMM）因其特有的双覆盖系统（数学覆盖和物理覆盖）在域离散方面具有独特的优势,而精细时间积分法则具有精度高、无条件稳定、无振荡以及计算结果不依赖于时间步长等特点。发展了用于研究二维瞬态热传导问题的精细积分NMM。结合待求问题的控制方程和边界条件,并基于修正变分原理导出了NMM的总体方程,给出了求解此类时间相依方程的精细时间积分及空间积分策略,选取了两个典型算例对方法的有效性进行了验证,结果表明本文方法可以高效高精度地求解瞬态热传导问题。     

瞬态热传导问题的精细积分数值流形方法研究

数值流形方法(NMM)因其特有的双覆盖系统(数学覆盖和物理覆盖)在域离散方面具有独特的优势,而精细时间积分法则具有精度高、无条件稳定、无振荡以及计算结果不依赖于时间步长等特点。发展了用于研究二维瞬态热传导问题的精细积分NMM。结合待求问题的控制方程和边界条件,并基于修正变分原理导出了NMM的总体方程,给出了求解此类时间相依方程的精细时间积分及空间积分策略,选取了两个典型算例对方法的有效性进行了验证,结果表明本文方法可以高效高精度地求解瞬态热传导问题。     

与结构动特性协同的自适应Newmark方法

提出了一种与结构动特性协同的自适应Newmark方法,其参数可基于数值弥散和数值耗散最小化的条件用解析方法求得.对线性单自由度动力学系统,该方法的相位误差精确为零并且谱半径为1.对线性多自由度系统和非线性系统,该方法在所有二阶积分解法中最精确.数值结果验证了新提出格式的高精度和结论.      

无条件稳定的更新精细积分方法

提出将Pade逼近与精细积分方法中的指数矩阵运算技巧结合起来,建立了精细积分法的更新形式及计算过程,对该更新精细积分方法的稳定性进行了论证与探讨.结果表明,该更新精细积分方法是无条件稳定的,整个积分方法的精度取决于所取Pade逼近的阶数与高斯积分点的数量.数值例题也显示了该方法的高效率及其可行性.     

基于位移型Gurtin变分原理计算动力响应的逐步积分法

本文利用位移型Ｇｕｒｔｉｎ变分原理,在时间域上采用三次Ｈｅｒｍｉｔｅ插值函数进行离散,给出了一种计算结构动力响应的逐步积分方法。通过稳定性分析研究了该方法的稳定区情况表明,当１．６４≤θ≤２．０８时,该方法的数值计算精度很高,但是条件稳定积分格式。当θ≥４．１时,该方法是无条件稳定的积分格式,精度较高。     

基于拉格朗日插值的无网格直接配点法和稳定配点法

由于无网格法中大多数近似函数均为有理式,不具有Kronecker delta性质,因此难以精确地施加本质边界条件.边界误差较大容易导致整个求解域求解结果精度低,甚至引起数值不稳定现象.文章在无网格直接配点法和稳定配点法中引入拉格朗日插值函数作为形函数,构建了拉格朗日插值配点法(LICM)和拉格朗日插值稳定配点法(SLICM).由于拉格朗日插值具有Kronecker delta性质,可以像有限元法一样简单而精确地施加本质边界条件,提高这两种方法的数值求解精度.稳定配点法基于子域对强形式方程进行积分,可以满足高阶积分约束,即可以保证形函数在积分形式下也满足高阶一致性条件,实现精确积分.同时,进行子域积分还可以减少离散矩阵的条件数,从而提高算法的稳定性.进一步提高拉格朗日插值稳定配点法的精度和稳定性.通过数值算例验证这两种方法的精度、收敛性和稳定性,结果表明基于拉格朗日插值的配点法的精度优于基于重构核近似的配点法,拉格朗日插值稳定配点法的精度和稳定性均优于拉格朗日插值配点法.     

单点子域积分与多点子域积分

基于单点子域精细积分的思想,针对抛物线型热传导方程初边值问题,提出了多点子域积分的概念,推出了一种多点子域积分的FTCS格式。该格式为显格式,并证明其为无条件稳定。数值算例表明,多点子域积分的FTCS格式具有比单点子域积分的FTCS格式收敛速度快的特点。     

A priori stability estimates and unconditionally stable product formula algorithms for nonlinear coupled thermoplasticity

This article describes new a priori stability for the full nonlinear systems of coupled thermoplasticity at finite strains and presents a fractional step method leading to a new class of unconditionally stable staggered algorithms. These results are shown to hold for general models of multiplicative plasticity that include, as a particular case, the single-crystal model. The proposed product formula algorithm is designed via an entropy based operator split that yields one of the first known staggered algorithms that retains the property of nonlinear unconditional stability. The scheme employs an isentropic step, in which the total entropy is held constant, followed by a heat conduction step (with nonlinear source) at fixed configuration. The nonlinear stability analysis shows that the proposed staggered scheme inherits the a priori energy estimate for the continuum problem, regardless of the size of the time-step. In sharp contrast with these results, it is shown that widely used staggered methods employing an isothermal step followed by a heat conduction problem can be at most only conditionally stable. The excellent performance of the methodology is illustrated in representative numerical simulations.     

可倾瓦径向滑动轴承支承的转子系统瞬态响应计算

研究了可倾瓦径向滑动轴承支承的转子系统瞬态响应的计算方法,提出了可倾瓦径向滑动轴承非线性滑膜力计算的数据库方法,将可倾瓦与轴颈的结构和运动参数归结到三个参数中,建立了起可倾瓦的非线性油膜力数据库,通过精确插值到每一运动参数组合下轴承提代的油力和轴瓦上所受的油膜力矩。     

瞬态热传导方程精细积分方法中对称性的利用

采用精细积分法求解瞬态热传导方程时,对指数矩阵进行变换后使其具有对称性,利用这一特性可使存贮量和计算量降低一半。变换后指数矩阵的带宽特性不变,采用子域精细积分可进一步提高算法的计算与存储效率。     

Transient heat conduction analysis using the NURBS-enhanced scaled boundary finite element method and modified precise integration method

The Non-uniform rational B-spline(NURBS)enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this paper.The scaled boundary finite element method is a semi-analytical technique,which weakens the governing differential equations along the circumferential direction and solves those analytically in the radial direction.In this method,only the boundary is discretized in the finite element sense leading to a reduction of the spatial dimension by one with no fundamental solution required.Nevertheless,in case of the complex geometry,a huge number of elements are generally required to properly approximate the exact shape of the domain and distorted meshes are often unavoidable in the conventional finite element approach,which leads to huge computational efforts and loss of accuracy.NURBS are the most popular mathematical tool in CAD industry due to its flexibility to fit any free-form shape.In the proposed methodology,the arbitrary curved boundary of problem domain is exactly represented with NURBS basis functions,while the straight part of the boundary is discretized by the conventional Lagrange shape functions.Both the concepts of isogeometric analysis and scaled boundary finite element method are combined to form the governing equations of transient heat conduction analysis and the solution is obtained using the modified precise integration method.The stiffness matrix is obtained from a standard quadratic eigenvalue problem and the mass matrix is determined from the low-frequency expansion.Finally the governing equations become a system of first-order ordinary differential equations and the time domain response is solved numerically by the modified precise integration method.The accuracy and stability of the proposed method to deal with the transient heat conduction problems are demonstrated by numerical examples.     

Development of a family of explicit algorithms for structural dynamics with unconditional stability

A new family of explicit integration algorithms is developed based on discrete control theory for solving the dynamic equations of motion. The proposed algorithms are explicit for both displacement and velocity and require no factorisation of the damping matrix and the stiffness matrix. Therefore, for a system with nonlinear damping and stiffness, the proposed algorithms are more efficient than the common explicit algorithms that provide only explicit displacement. Accuracy and stability properties of the proposed algorithms are analysed theoretically and verified numerically. Certain subfamilies are found to be unconditionally stable for any system state (linear elastic, stiffness softening or stiffness hardening) that may occur in earthquake engineering of a practical structure. With dual explicit expression and excellent stability property, the proposed family of algorithms can potentially solve complicated nonlinear dynamic problems.