分数阶微分方程的Haar小波算法研究

在BPFs的Caputo分数阶微分算子矩阵的基础上,建立了Haar小波的分数阶微分算子矩阵,提出了一种有效的求解分数阶微分方程的Haar小波数值方法,并将该方法应用于线性和非线性分数阶常微分方程求解中.数值算例表明,该算法简单,数值精确度高,是一种高效的数值求解方法.     

奇异Hamilton系统矩阵的精细积分法

常规位移有限元的结构振动方程是n个二阶常微分方程组.采用一般交分原理推导,将结构振动问题引入Hamiltoil体系,将得到2n个一阶常微分方程组.精细积分法宜于处理一阶方程,应用于线性定常结构动力问题求解,可以得到在数值上逼近精确解的结果.对于非齐次动力方程,当结构具有刚体位移时,系统矩阵将出现奇异.本文借鉴全元选大元高斯-约当法求解线性方程组的经验,提出全元选大元法求奇异矩阵零本征解的方法,该方法可以简便快速地寻求奇异矩阵零本征值对应的子空间.利用Hamiltoil体系已有研究成果及Hamilton系统的共轭辛正交归一关系,迅速将零本征值对应的子空间分离出来,通过投影排除奇异部分,然后用精细积分法求得问题的解.数值算例表明,该方法对Hamilton系统奇异问题,处理方便,计算量小,易于实现,同时保持了精细算法的优点.     

求解结构动力学方程的一种辛格式及其优化

论文将四阶隐式高斯勒让德辛龙格库塔法应用于线性结构动力学方程,并对其进行了算法优化.针对n个自由度的动力学初值问题,先通过消元得到n阶线性代数方程组,利用其系数矩阵稀疏对称正定的性质,采用预处理共轭梯度法求解,其中预条件子由系数矩阵的不完全Cholesky分解得到.通过与中心差分法、Newmark-β法及Runge-Kutta法相比,论文方法在计算量未显著增加的前提下给出了更高的计算精度.     

数字体图像相关方法中基于迭代最小二乘法的物体内部变形测量

本文提出一种基于迭代最小二乘法的亚体素位移测量算法.该算法模型结合了符合实际情况的线性灰度变化模型和线性位移映射函数,采用仅需变形后体图像一阶灰度梯度的迭代最小二乘法计算亚体素位移和位移梯度.由于该算法采用的模型符合实际情况,因而具有更高的位移测量精度和更广泛的适用性.此外,由于仅需变形后体图像的一阶灰度梯度,无需计算...     

松弛方法在计算凝聚炸药爆轰问题中的应用

利用松弛近似,将非线性的凝聚炸药爆轰控制方程转化为线性的松弛方程组,并采用五阶WENO格式和五阶线性多步显隐格式对线性松弛方程组进行空间方向和时间方向的离散,由此建立具有高精度和高分辨率性质的计算凝聚炸药爆轰的松弛方法。建立的松弛方法可以避免求解Riemann问题及计算非线性通量的Jacobi矩阵,同时无需分裂处理反应源项。通过对凝聚炸药的平面一维定常爆轰波结构及球面一维聚心、散心爆轰起爆和传播过程的数值模拟,验证了所建立的松弛方法能够很好地计算凝聚炸药爆轰问题。     

平面亚声速定常势流的最小二乘有限元法

1.方法 平面亚声速定常势流的无量纲基本方程组是其中,x,y是以流场特征长度无量纲化的直角坐标;u,ν和α是以远方来流声速无量纲化的速度分量和局部声速;M_∞是远方来流马赫数;r是比热比。 边界条件视具体问题而定。我们用迭代法求拟线性一阶偏微分方程组(1)的数值解。在迭代的每一循环中,用最小二乘有限元法解线性化的方程组。我们试验了变刚度与等     

基于非线性状态空间辨识的气动弹性模型降阶

高维、非线性气动弹性系统的模型降阶是当前气动弹性力学与控制领域的研究热点之一.然而国内外现有的非线性模型降阶方法仍存在辨识算法复杂、精度有待提高等问题.本研究提出了一种基于非线性状态空间辨识的跨音速气动弹性模型降阶方法.首先,该方法基于非定常空气动力的单位脉冲响应数据,采用特征系统实现算法对非线性状态空间模型的线性动力学部分进行系统辨识.其次,引入状态和控制输入的非线性函数,采用优化算法对非线性函数的系数矩阵进行优化,进而得到考虑非线性效应的空气动力降阶模型.为了验证该降阶模型在预测跨音速气动弹性力学行为的精确性,本文以三维机翼为研究对象,分别从基于非线性降阶模型的气动力辨识、跨声速颤振边界计算和极限环振荡预测三方面进行了算例验证,并与现有的模型降阶方法进行了对比,进一步说明本文所提出方法的有效性.研究结果表明,该降阶模型对上述三类问题的计算精度与直接流-固耦合方法相吻合,可用于高效预测飞行器跨声速气动弹性力学行为.     

基于非线性状态空间辨识的气动弹性模型降阶

高维、非线性气动弹性系统的模型降阶是当前气动弹性力学与控制领域的研究热点之一.然而国内外现有的非线性模型降阶方法仍存在辨识算法复杂、精度有待提高等问题.本研究提出了一种基于非线性状态空间辨识的跨音速气动弹性模型降阶方法.首先,该方法基于非定常空气动力的单位脉冲响应数据,采用特征系统实现算法对非线性状态空间模型的线性动力学部分进行系统辨识.其次,引入状态和控制输入的非线性函数,采用优化算法对非线性函数的系数矩阵进行优化,进而得到考虑非线性效应的空气动力降阶模型.为了验证该降阶模型在预测跨音速气动弹性力学行为的精确性,本文以三维机翼为研究对象,分别从基于非线性降阶模型的气动力辨识、跨声速颤振边界计算和极限环振荡预测三方面进行了算例验证,并与现有的模型降阶方法进行了对比,进一步说明本文所提出方法的有效性.研究结果表明,该降阶模型对上述三类问题的计算精度与直接流-固耦合方法相吻合,可用于高效预测飞行器跨声速气动弹性力学行为.     

二维双曲型守恒律的一种五阶松弛格式

松弛格式是Jin和Xin提出的无振荡有限差分方法,其主要思想是将守恒律转化为松弛方程组进行求解.本文用逐维五阶WENO重构和显隐式Runge-Kutta方法对松弛方程组的空间和时间进行离散,得到了一种求解二维双曲型守恒律五阶松弛格式.所得格式保持了松弛格式简单的优点,不用求解Riemann问题和计算通量函数的雅可比矩阵.通过二维Burgers方程和二维浅水方程的数值算例验证了格式的有效性.     

多体动力学超定运动方程广义-α求解新算法

完整约束多体系统第一类Lagrange方程建模得到的运动方程是指标-3形式的微分-代数方程(differental-algebraic equations,DAEs).如果同时考虑速度约束,将得到超定运动方程,该方程是指标-2的超定微分-代数方程(over-determined differential-algebraic equations,ODAEs).基于结构动力学中常用的广义-α方法,将其拓展,求解包含速度约束的超定运动方程,相对于其他求解指标-2 ODAEs的算法,新的算法没有增加离散得到的非线性方程组方程的数目.通过数值实验验证算法,并说明其求解ODAEs不存在精度降阶的现象,仍然具有二阶精度,同时算法的数值耗散也是可以控制的.最后新方法与其他求解多体系统ODAEs形式运动方程算法的CPU时间进行了比较分析.