大型稀疏复线性方程组双共轭梯度法

有限元复线性方程组的系数矩阵一般具有稀疏性和对称性的特点,全稀疏存贮方法就是利用这些特点,只存贮对称部分的非零元素,采用链表式管理,即节省存贮空间,又便于动态更改.在一般双共轭梯度法的基础上,本文利用广义变分原理对内积进行了重新定义,使双共轭梯度法求解复线性方程组更为有效.数值算例表明这种双共轭梯度法结合全稀疏存贮方案的求解算法在时间和存贮上都较为占优,可靠高效,能够应用于有限元线性方程组的求解.     

大型稀疏线性方程组改进的ICCG方法

有限元分析求解的大型线性方程组．其系数矩阵中只有少量的元素不为零，是一个稀疏矩阵，而且一般都是对称的。求解过程中只存贮对称部分的非零元素，这种全稀疏存贮策略采用链表式管理，使存贮结构可以动态更改。在带双门槛值ICCG方法的基础上加上对角元修正策略，产生了一种改进的ICCG求解法，它不仅节省了存贮空间，而且便于元素的查询、插入和删除，提高了求解效率。这里也给出了验证算例。     

隐格式并行直接求解方法研究

讨论了由CFD问题得到的大型稀疏线性方程组的并行直接法求解问题.介绍了三对角型方程组的SPP算法,将之推广来求解一般的带宽较窄的带状或者稀疏带状线性方程组.介绍一般稀疏的线性方程组的基于排序后再分解的并行求解方法的步骤及其中的一些关键问题.     

离散单元节点编号的优化

一、引言 用有限元素方法分析构件时,通常要产生一个高阶的总刚度矩阵[K],它是一个稀疏对称、带状分布的矩阵。在求解表示刚度矩阵[K]、节点力矢量{Q}以及节点位移矢量{U}之间的平衡关系的线性代数方程组:[K]{U}={Q}时,所需要的计算时间及计算机的存贮量,在矩阵大小一定的条件下,主要取决于总刚度矩阵[K]的带宽。计算时间和存贮量与带宽的平方成正比。     

高阶满阵线性代数方程组的快速解法及其应用研究

在非对称、非正定满阵线性方程组的求解过程中,本文以高斯迭代为基础,提出一种残差优化的迭代法.该法在迭代收敛过程中,残差单调下降,大大缩短了线性方程组的求解时间.用该法替代直接解法,增大了求解方程组的阶数,缩短了计算周期,对复杂飞行器气动、隐身一体化优化设计具有重要的工程价值.该方法还可以广泛地推广到其它应用领域.     

多介质数据库管理系统NAIMDB——模型与方法

采用面向对象的方法构造数据模型,把命名的线性表定义为简单对象。用表可以嵌套的性质构造复杂的多介质对象。对象分为3个层次,即原子对象,简单对象和复杂对象。简单对象在数据库中有实际的物理存贮;复杂对象在数据库中只存贮它的结构描述,即相应嵌套表的结构描述。     

环形过渡元素及其在轮盘振动分析中的应用

 本文发展了联结环形等参数元素和环形厚板元素的一类环形过渡元素,它对于轴对称体的结构动力分析具有重要意义。环形元素采用了半解析的方法,把三维问题降为二维问题来处理,大大节省了计算机存贮量及计算时间。文中对某实际涡轮盘的振动特性进行了计算,其结果是令人满意的。     

折叠翼变体飞机的传动机构设计和驱动方式优化分析

变体飞机能够在飞行过程中显著地改变外形和尺寸,从而使单一的飞行器能够实现在起飞、俯冲、巡航等多种飞行状态下都具备优良的飞行性能.折叠翼变体飞机的折叠传动机构对折叠翼的结构动力学和传力特性具有显著的影响.本文设计了一种轴/轴套式的折叠传动机构,利用ANSYS有限元软件建立了有限元网格模型,采用DYNA求解器,详细求解分析...     

二维非结构网格上的高阶间断有限元方法研究

对二维非结构网格上的求解欧拉方程的高阶间断有限元方法进行了研究.运用间断有限元基本理论,采用多项式基函数对解进行近似描述,欧拉方程的数值通量项运用高斯积分公式和Roe格式求取,时间推进采用传统四步Runge-Kutta方法.运用发展的间断有限元方法对NACA0012翼型在亚跨音速情况下的无粘流动和圆柱的低速无粘流动进行了数值计算,结果表明:间断有限元方法具有良好的收敛特性、很小的数值耗散和很好的激波捕捉能力.     

稀疏孔径ISAR机动目标成像与相位补偿方法

现代逆合成孔径雷达(ISAR)可实现多模式下对多目标的交替测量工作,但波束切换使得任一目标的方位孔径具有稀疏性,同时目标的机动性也使得目标的方位维信号具有高次调制特性,这均导致以连续采样信号模型为基础的传统ISAR运动补偿和成像方法难以直接应用。本文提出一种新的稀疏孔径(SA)ISAR机动目标相位补偿和成像方法。该方法首先建立了稀疏孔径ISAR机动目标的信号模型,然后通过求解2范数最大化的优化函数,实现对相位误差的精确估计,并且通过迭代处理提高了相位误差的估计精度。相位误差补偿后,根据压缩感知理论,构造了机动目标情况下的Chirp-Fourier字典,通过对稀疏优化问题的求解,能够较精确地从稀疏孔径信号中重构出机动目标的全孔径(FA)信号。仿真和实测数据处理结果证明了本文方法的有效性。     

有限元法求解QPNS方程模拟高速旋成体气动力和热

从ＦＮＳ方程中略去流向粘性导数项导出熵变量形式ＱＰＮＳ方程（拟抛物化Ｎ－Ｓ方程）其离散解能自动满足熵增不等式。利用ＧＬＳ（伽辽金／最小二乘法）有限元法构造出ＱＰＮＳ方程的弱解式,所形成的整体方程组用ＧＭＲＥＳ（广义最小余数法）求解。采用空间推进求解ＱＰＮＳ方程,能明显降低计算机内存和运算时间。算例给出典型旋成体的物面压强、剪应力、热流量分布以及流场中密度、温度分布。     

不对称推力完全自动补偿技术

研究了一种适用于多发飞机的不对称推力完全自动补偿策略,并针对发动机失效后的不对称推力飞行特性进行研究。首先,推导了横航向达到平衡状态的条件,提出了一种基于发动机转子转速信号的自动补偿控制方法,在保证补偿效果的同时,解决了自动补偿系统可靠性低的问题,该自动补偿控制方法将高压转子转速信号与低压转子转速信号进行逻辑运算,结合高压转子转速差控制副翼和方向舵偏转,再与主飞行控制系统共同作用使飞机达到平衡状态。然后,基于横航向主飞行控制系统具有滚转角保持功能的特点,优化自动补偿控制系统。最后,全包线内建立飞机达到平衡状态的线性化数学模型,设计了不对称推力完全自动补偿控制律。通过MATLAB/Simulink建模并进行仿真验证,仿真结果表明完全自动补偿控制方法对不对称推力飞行具有良好的补偿效果。     

用有限元法分析接触问题的一种新的力学模型及计算方法

阐明了用有限元法计算物体接触问题的一种新的力学模型:主动接触体网格中一个节点与被动接触体单元中任意点相接触的模型,放弃了过去常用的节点对模型,为网格划分带来了便利。在这种模型下,给出了粘式、摩擦滑移、开式等接触边界条件。按照修正势能原理导出了对称的接触控制方程,可求解接触体具有非线性的静态/动态的接触问题。本文给出了在一次求解迭代后判断各节点新的接触状态的方法,可进一步决定是否和如何修改接触控制方程。      

基于模态耦合方法的机翼颤振数值模拟

针对机翼气动弹性问题,发展了一种流固耦合模拟技术,将固体结构的广义位移表示为结构固有模态的线性组合,求解结构的运动,再与流体进行数据耦合、迭代计算.同时为了确定机翼颤振边界,发展了一种线性插值确定颤振边界方法.采用标准气弹机翼AGARD445.6算例,验证所采用的模态耦合技术和颤振边界确定方法的有效性及工程实用意义.     

榫头/榫槽接触问题边界元分析

本文用增量理论建立了二维弹性接触问题的边界元法方程,讨论了求解接触问题过程中的一些具体做法。文中分别对盘/片组件的燕尾型与枞树型榫头与榫槽的接触应力作了数值分析。与有限元及光弹试验结果对比表明,边界元法的结果是令人满意的。      

求解循环对称结构的伪单元法

 对于循环对称结构的求解,在引入循环对称条件时,需对有限元总刚进行行变换和列变换,这将改变总刚的带宽和非零元素的位置。本文提出了求解循环对称结构的伪单元法,用三种伪单元实现上述各变换项的改变量,无需形成完整的总刚矩阵,从而大大地提高了计算机的解题能力。文中用算例进行了说明。     

梁/液体耦合振动的边界积分对称列式法

 利用格林积分导出的液体内部扰动压力函数作为测试函数,通过变分得到一组具有对称系数的代数方程,将该方程同梁的有限元振动方程联立后,即可求解梁在液体中的耦合振动问题。本方法具有未知数目少,数据存储和计算效率高的特点。     

标量图像测速原理及数值检验

标量图像测速是一种新型的能够测量微小结构的湍流实验测量系统,其理论依据是Schmidt（Sc）数〉1的标量湍流里,标量场脉动的时间尺度和空间耗散尺度要小于对应的速度场脉动的时间尺度和空间耗散尺度,标量脉动比速度脉动具有更强的间歇性,标量场所含信息大于速度场所含信息,从高＆数标量场提取速度场是可能的。笔者详细介绍了标量图像测速原理,并在积分最小化法思想指导下根据标量输运守恒方程及限定条件,建立了积分最小化标量图像测速的数学表示。通过变分法,获得了积分空间每个点的三个速度分量的欧拉特征方程,对欧拉特征方程离散化即可在积分空间建立一线性稀疏方程组。采用变分迭代法对此线性稀疏方程组进行求解即可获得三维速度场,据此采用积分最小化法建立了标量图像测速图像处理系统。还介绍了四维标量场（三维空间加一维时间）的图像采集原理,并采用DNS数据对标量图像测速处理系统进行了数值检验,检验结果表明标量图像测速技术所求解的速度场和DNS精确解之间在结构上是十分相似的,标量图像测速可以正确地提取速度场。