一种无条件稳定的时域推进法与有限元混合分析瞬态热传导问题

本文将一种新的无条件稳定的时域推进法与有限元结合,用于分析瞬态热传导问题。根据微-积分型热传导控制方程,对时间变量在小区闻内插值积分,变初-边值问题为一系列离散时刻的边界值问题,再应用有限元法求解之。由于推进中的每一时刻的解,都严格满足原给定的初始条件,就消除了累积误差的影响。此外,这种时域推进法对时间变量的数值积分也优于数值微分。因而,与现行的基于微分型控制方程,采用时-空有限元模型同时离散,再逐步求解的直接积分算法相比,计算精度可大大提高。不难预料,当求解较长时间后的瞬态值时,本算法的优越性会更加明显。     

时域推进与边界元混合瞬态热传导问题

对时间变量在小区间内应用多项式插值,对空间变量采用边界元求解,提出一种新的分析瞬态热传导问题的计算格式。与现行的利用含时间变量的基本解所得到的边界元计算格式相比较,具有基本解简单易得、系数矩阵容易计算以及无重复积分等优点;与向前差分消去时间因子的边界元计算格式相比较,则在算法稳定性、计算精度以及时间步长的选取等方面具有明显的优势,并且很容易推广到求解非线性问题。算例表明,采用时域推进与边界元混合方法分析瞬态热传导问题有效可行。     

有限元的一种并行算法及实现

当前,由于有限元素法大都使用串行算法,处理时间较长。本文根据有限元素法的特点和采用的并行计算机结构,提出了一种并行有限元算法,能大大提高其处理速度,加快处理周期。且算法简单,实现方便。     

Navier-Stokes方程两层稳定有限元算法分析

分析了定常Navier-Stokes方程的两种两层稳定有限元算法。它们将局部Gauss积分稳定化技术和两层算法的思想充分结合,采用低次等阶有限元P1-P1或Q1-Q1对N-S方程进行数值求解。误差分析和数值算例结果表明,当粗、细网格尺度H=O（h1/2）时,它们与在细网格上的单层有限元算法具有相同的收敛速度,而两层算法却节省了大量的计算时间。相比之下,Simple算法具有更高的计算效率。而且进一步发现Oseen算法能够对小粘性系数N-S方程进行有效求解。     

一种时域稳定参数预报与最优估计兼容的参数辨识法

 本文导出一种基于时域稳定的参数预报与最优估计兼容的状态估计和参数辨识算法。这种算法将系统辨识分为三步:1.以辨识系统预报误差渐近稳定为准则的参数预报;2.以预报参数为条件的状态估计;3.对参数和状态的后验修正。即一般分割辨识算法（GPIA）与模型参考辨识算法（MRIA）相结合的兼容辨识算法。本文将这一算法用于飞行器气动系数和控制导数的辨识,并与GPIA算法所得结果相比较。     

一种中心型间断有限元MMALE方法

针对多介质可压缩流体动力学问题,提出了一种单元中心型二维MMALE(Multi-Material Arbitrary Lagrangian-Eulerian)方法.在拉氏步,流体力学方程组采用中心型间断有限元方法求解.对于混合网格,采用Tipton压力松弛模型更新物理量,用等参坐标法更新物质中心点坐标.界面重构采用一种...     

一种蜂窝夹层有限元

 在航空结构,特別是翼面结构中,随着复合材料的应用,蜂窝夹层结构应能得到更加广泛的应用。我们研制的金属及复合材料应力分析程序——YJS851中发展了一种蜂窝层元素,通过不同算例考核,其结果是满意的。现已用于某型飞机的打样计算。用美国的MSC/NASTRAN程序对同一模型作了计算,二者结果几乎没有差别。同时与MSC/NASTRAN相比,我们发展的程序使用方便,价格便宜。     

一种特殊梁的有限元分析

通过对含有裂纹的梁局部的特殊处理, 得出含裂纹单元的应力强度因子, 建立了该单元刚度矩阵, 进而建立了梁整体的有限元分析平衡方程。     

三维编织复合材料热物理性能的有限元分析

根据三维编织复合材料的细观结构,分别建立了三维四向和五向编织复合材料热物理性能的有限元模型.采用周期性的非绝热温度边界条件和位移边界条件,计算了三维四向和五向编织复合材料的整体等效热传导系数和热膨胀系数,计算结果同已有文献相比与实验值符合得更好.在此基础上,迸一步研究了编织角、纤维体积分数、编织结构等参数对材料热物理性...     

TC11钛合金压气机叶片锻后冷却过程瞬态温度场的数值模拟与试验研究

运用二维有限元法和试验方法,模拟、分析了TC11钛合金压气机叶片锻后冷却过程的传热现象,得到该叶片各个翼形截面上瞬态温度场的时变特性和分布规律.模拟结果为解决该类叶片锻后冷却过程翘曲变形提供了可靠依据.     

一种考虑结构瞬态互辐射的有限元法

研究结构元件之间瞬态辐射过程及特点 ,对于元件温度未知时的瞬态辐射提出了一种有限元计算模型。通过广义边界条件简化了元件互辐射的数学描述 ,最后得到了广义互辐射元。所建模型具有简单、方便、实用和精度高等特点 ,使复杂的互辐射计算化为统一模式 ,并通过试验与实例计算证明其有效性。     

一种新的有限元线性方程组的求解算法

本文采用链表数据结构存储有限元结构总体刚度矩阵,对一有限元结构分析系统进行了改进。算例表明,改进后的系统可使存储总体刚度矩阵所需内存及分解总刚阵时间分别下降到原系统的1/2～1/4,大大提高了计算机的解题能力和执行效率。该系统可在微机上进行大型复杂结构的受力分析。     

一种同位网格上的迎风有限元方法

通过Ｇｒｅｅｎ定理将对流项变量从微分算子中分离出来,从插值函数入手引入迎风格式,是对强对流问题有限元计算中对流项变量的一种新的处理方法。按这种方法采用局部斜迎风格式及速－压同位网格公式,构成了一种对高Ｒｅｙｎｏｌｄｓ数流体流动数值模拟比较有效的有限元方法。     

控制体有限元方法中的一种高阶插值格式

针对控制体有限元方法对流动方向不平行于坐标轴方向的回流问题数值模拟采用各种迎风格式都容易产生假扩散的问题,本文给出了一种由插值单元及其迎风相邻单元所确定的高阶插值格式。通过对Re=1000时三种不同形状的二维四边形空穴流模拟,表明该格式与局部斜迎风格式的适当组合可有效地减小控制体有限元方法对回流问题数值模拟的假扩散。      

一种可行的有限元网格优化法

在用有限元素法进行分析时的一个主要考虑是,对于某个特定问题,如何才能以尽可能少的计算工作量获得尽可能高的计算精确度,这就是在过去十年中迅速发展起来的有限元最优离散化问题,并已提出了不少优化准则和算法。本文的目的在于提出一项改善有限元离散化的有效方法,以避免现有方法中的一些困难和不足,即:既可以避免纯数学规划法中过于庞大的计算时间,又可以避免批处理方式的麻烦;它编制程序简便,也无需计算等高线的专用软件以及与计算机接口的图象显示设备。它从误差分析得到使离散化近乎最优的可行方向,再通过一维搜索得到近似于最优的离散化。数值例题表明,采用本文方法进行优化,得益是十分显著的。     

一种求取发动机状态变量模型的改进拟合法

采用的拟合辨识方法是直接对根据状态变量模型中各变量间的线性关系对系数矩阵进行拟合,并且实现对稳态过程和瞬态过程进行分步拟合,从而消除瞬态误差对稳态精度的影响.本方法有接近传统最小二乘法的精度,并且无需求解状态变量模型微分方程或使用迭代算法,只需要进行简单的矩阵运算即可得到各系数矩阵,因此计算效率远高于传统最小二乘法.同时,该方法也可与其他拟合方法结合使用从而进一步提高拟合精度.最后,应用该方法建立了某涡轴发动机的小偏差状态变量模型,通过与非线性模型仿真结果比较,验证了该方法的有效性,即使对象非线性较强,该方法也能得到精度较高的拟合效果.      

一种基于高阶剪切变形板模型的有限元算法

为了较好地模拟复合材料层合板的变形并得到较为精确的应力分布,本文基于Reddy提出的高阶剪切变形板模型建立了一种新型有限元。该模型保证横向剪应力沿板厚度方向呈抛物分布,在板的上下表面为零,从而很好地满足力边界条件。有限元算法构造的出发点是修正的虚功泛函,引入Lagrange乘子放松了挠度在单元边界之间的连续性,使得挠度的构造变得简单;由于通过常应力分片检查,所以有限元算法的收敛性得到了保证。此外,本算法的数值结果具有对单元网格形状不敏感、精度好和计算效率高的性质。