多体系统动力学子系统求解算法

为了降低系统求解规模,提高复杂多体系统动力学求解效率,研究了多体动力学子系统求解算法.通过分析系统拓扑构型,基于雅克比矩阵分块特性,提出子系统综合算法.算法采用半递归法建立含开环、闭环的子系统运动方程,将子系统质量和外力等效作用在内接基体上,解耦系统运动方程组装和加速度求解.此外,针对松散耦合子系统,提出了予系统分治算法.算法针对子系统间的耦合特点,分解复杂多体系统模型,实现子系统的独立求解.最后以机构动力学数值求解实例,验证上述算法的正确性和求解效率.     

基于局部DQ-RBF不可压缩N-S方程的数值求解

以RBF作为DQ方法的基函数,将迎风机制引入DQ-RBF中,建立了二维不可压缩黏性N-S方程数值求解模型,采用Levenberg-Marquardt算法求解非线性方程组.求解时分析了形状参数对求解精度的影响,改进了边界速度的处理方法.对平板Couette流及有限宽台阶绕流流动问题进行了数值求解.比较了本文方法和FLUE...     

数值求解不可压粘性流体定常运动的格林函数方法

本文提出了一种数值求解不可压粘性流体定常运动的格林函数方法.在本文中利用Stokes方程的基本解作为格林函数将求解不可压粘性流体定常运动的边值问题化为求解速度场和边界应力的非线性积分方程组,在解出速度场和边界应力后可直接计算流场中各点的压力;用有限元近似将积分方程离散化而进行数值求解。对于小雷诺数流动,只归结为求解边界积分方程,使求解区域减少一个维度。对于非线性问题,可用迭代方法求解,在每次迭代中只须解出边界点上的速度或应力。通过几个简单的算例,表明本文所提出的方法具有精度高、处理边界条件简单、通用性强的优点,并具有求解各种复杂流动的潜力。     

一种求解复波数域波动弥散方程的通用数值方法

波动弥散方程的求解是分析弹性波或声波问题时的一种比较常见的过程,求解所得的曲线即为弥散曲线。其方程形式一般是关于波数与频率的一个二元超越方程,当波数为复数时,方程实际为三元超越方程。这类方程没有解析解,由于方程的变量和参数涉及复数域,一般的数值求解方法也难以处理。本文提出了一种基于模值收敛性的数值方法,求解一般波动弥散方程复波数域内的解。首先论述了本方法的理论依据,然后介绍了其求解步骤,最后将本方法应用于求解复波数域波动弥散方程的具体算例。通过对比可以发现:本方法广泛有效地解决了复波数域内波动弥散曲线的求解问题,适用于不同形式波动弥散方程的求解。     

基于ODE求解器的高墩渡槽横向动力分析

本文将槽墩视为连续化的悬臂梁,将槽内水体用Housner模型简化,建立了高墩渡槽计算模型,并首次将基于ODE求解器的特征值算法应用于高墩渡槽的动力特性分析中,求解了边界条件含特征值的ODE特征值问题,还利用基于ODE求解器的振型叠加法,求解了高墩渡槽的横向弹性动力响应.     

有理宏单元法求解泊松方程

基于有理超限插值,提出了一种在求解域边界布点的全域求解数值方法——有理宏单元法。推导出了三角形及四边形单元的有理混合函数,划分单元各边的节点并选定各边上的插值函数,建立了三角形及四边形母单元的形函数。利用等参变换,将求解域影射到相应的母单元上,得到了求解泊松方程边值问题的有理宏单元方程组。通过将求解域划分为一个或多个宏单元,有理宏单元法可对任意形状的二维区域求解。作有理宏单元法解泊松方程边值问题的算例,验证了本文方法的有效性。     

对称模糊优化问题的直接求解法

通常,对称模糊优化模型的求解是一个迭代过程,要多次求解规划问题。本文根据对称模糊最优判决准则,给出一种直接求解法,只需求解一次规划问题,即可求出问题的最优解。     

Hermite梯度重构核近似配点法及其在功能梯度材料板的静力分析中的应用

由于直接配点法在求解边值问题时边界上的求解精度较低,本文提出了Hermite梯度重构核近似配点法（HGCM）来改进边界求解精度。重构核近似是无网格法中一种常用的近似函数,但是其在求解高阶导数时格式复杂且非常耗时。HGCM采用梯度重构核近似构建形函数的任意高阶导数,提高了计算效率;通过Hermite配点法构建离散方程,提高了边界求解精度。这种方法在求解对应变系数四阶偏微分方程的功能梯度材料板的静力问题时精度高,计算效率高,并可进一步推广应用于高阶偏微分方程描述的边值问题。     

无穷域势流问题的有限元／差分线法混合求解

提出了一种将有限元和差分线法相结合求解无穷域势流问题的算法。用两同心圆将求解域划分为存在重叠的有限和无限两个区域,在有限和无限域上分别用有限元和差分线法求解Laplace方程边值问题。用差分线法推导出的关系式修正有限元方程,求解该方程组从而得到原问题的解。本算法将求解无穷域问题转化为代数特征值问题和有限域内线性方程组的...     

带有襟翼偏转的三维机翼绕流的N—S方法解

给出了带襟翼偏转的三维机翼绕流的一种求解Ｎ－Ｓ方程的计算方法,采用区域求解算法和对接分区网络技术相结合的方法,有效地求解了绕此外形的复杂流动,区域求解算法中提出了一种满足通量守恒的内边界耦合条件,流场求解时采用中心差分的限体积方法对空间通量顶进行离散,采用显式推进方法进行时间方向的积分,数值算例表明本方法是求解带襟翼偏转的机翼绕流的有效方法。     

适用于IBM-PC微机的有限元程序

程序1:用波阵法求解有限元的线性方程的程序——BZH程序。BZH程序用BASIC语言编写而成,可用于求解有限元法的线性方程组。波阵法的特点是占用内存很小。求解过程中边组装边消去,充     

有限变形问题过程最优控制算法

在非增量算法的基础上,提出了用最优控制变分原理形成过程最优控制迭代求解的基本思路,并给出求解的基本控制方程。这一工作为有限变形力学问题的数值求解提供了一个新的处理方法。     

求解二维弹性问题的径向基差分法

常规的配点型无网格法在求解弹性力学问题中,存在求解精度差和纽曼边界条件处理等局限．为解决这一问题,通过利用流体力学中基于径向基构造的差分格式(RBF-FD),来求解弹性力学平面问题．同时,为了进一步提高求解精度,对纽曼边界条件采用Hermite插值进行处理．数值算例表明,该方法具备良好的收敛性,并有着较高的精度,可有效解决传统配点型无网格法精度差的问题．同时,也表明该方法可以应用于弹性力学问题的求解．     

平面粘弹性问题的辛求解方法

将弹性力学辛对偶求解方法与Laplace变换相结合,提出了一个求解粘弹性平面问题的新方法。首先利用Laplace变换,将粘弹性平面问题转化为一个准弹性问题,在辛弹性力学的框架下,利用分离变量和辛本征展开法对其进行求解,然后由逆变换得到原问题的解。为证明方法的有效性,求解分析了矩形域平面粘弹性圣维南问题,得到了令人满意的结果。     

一种几何大变形下的非线性气动弹性求解方法

非线性气动弹性的时域求解中, 涉及到非线性的流体动力学(CFD)和非线性的结构动力学(CSD)耦合问题. 基于Co-rotational理论, 推导了三维壳单元几何非线性下的切线刚阵和内力公式, 针对推进过程中的能量守恒, 引入预估-校正推进格式, 发展了一种近似能量守恒的非线性动态响应算法; 基于1/2时间步的交错耦合格式, 结合带有几何守恒律的双时间推进求解雷诺平均N-S方程的求解器, 发展了非线性气动弹性求解的高精度耦合格式. 通过结构几何大变形下的静力和动力分析验证了所发展的结构非线性求解器, 并通过AGARD445.6机翼的非线性气动弹性响应分析, 说明了所发展耦合求解方法的实用性.      

欧拉-伯努利梁的静态响应分析

欧拉-伯努利梁,即通常所称的工程梁或弯曲梁(不考虑剪切变形和转动惯量),其静态响应分析包括求解梁的挠度、转角、弯矩和剪力,亦即求解梁的变形和内力。求解梁的挠度的方法很多,诸如能量法、力矩面积法、差分法等,前不久有人撰文“用奇异函数法求解某些变截面梁的变形”。本文用传递矩阵法分析     

用材料力学解的切应力求解弹性力学问题

分析了切应力与正应力的关系,讨论了导出切应力公式的条件,提出按切应力求解弹性力学问题的方法.证明凡是σ_y与x无关的梁或者已知切应力为零的问题,按切应力求解都是可行的.用这种方法求解比传统方法方便,运算简单.该文给出了用切应力求解弹性力学平面问题的两个算例.     

Duffing方程强迫振动问题的插值摄动解法

本文用插值摄动法求解了一类弱非线性非自治系统有正阻尼的Duffing方程的初值问题。此法把一个求解非线性二阶常微分方程问题转化成为求解二次代数方程的问题,计算过程十分简单,并能反映振幅的变化,因而优于L－P解法。     

求解运动学问题的运动链法

详细阐述了求解理论力学的运动学问题的运动链法. 首先进行机构分析,然后画出运 动链图,根据运动链图确定求解方案,然后实施该方案. 列举了两个典型例题说明使用该方 法求解运动学问题的过程.     

结构模态分析的有限元并行模拟

以子结构模态综合分析为基础,提出一种求解大型结构特征值问题的并行解法．采用子结构模态综合算法,结构特征模态采用子空间迭代方式并行求解．这种子空间迭代法的子结构并行计算的实施是利用子结构的刚度阵和质量阵而不必完全组集系统刚度阵和质量阵求解综合系统的特征值问题．数值结果表明这种求解大型结构特征值问题的并行算法是可行有效的．