波动方程差分格式的优化

用差分方法模拟非定常流动和气动声学问题时,重要的一点是使差分格式的色散关系尽量与原波动方程一致．文中总结了建立差分方程色散关系的各种方法,以积累误差为优化目标函数对差分格式进行了优化分析,有效地控制了差分格式长时间计算的精度与稳定性问题,证明时间积分采用四级Ｒｕｎｇｅ－Ｋｕｔａ法较之三级Ｒｕｎｇｅ－Ｋｕｔａ法更适合于将单波方程优化格式推广至方程组,以Ｏｓｈｅｒ－ＣｈａｋｒａｖａｒｔｈｙＴＶＤ格式和ＰＦＤＤ格式为例,给出有限面积格式的优化方法,并以实例证明其有效性．     

线性状态方程的精确计算

以解析法为基础,采用新的精确计算方法求解线性系统的状态方程．该方法属于一种２Ｎ类算法,它是把每一时间步长再进一步细分为２Ｎ个精细步长,然后在精细步长内利用泰勒展开式求解矩阵函数．与传统的Ｒｕｎｇｅ－Ｋｕｔａ方法比较,这种方法兼有解析法与数值法的优点,稳定性好,可以采用较大的时间步长以提高计算效率．     

求解结构振动非线性特征值问题的一个方法

本研究结构分析中动态有限元引起的非线性特征值问题。中提出将动态有限元法与Ｒｉｃａａｔｉ传递矩阵法相结合求解振动非线性特征值问题。该方法减少了普通动态有限元法中结构划分所需的节点数，大大扩大了在微型机上的求解问题的规模；对传递矩阵引入Ｒｉｃｃａｔｉ变换，从而比标准传递托 阵法减少了误差传递，提高了计算效率；采用Ｎｅｕｔｏｎ－Ｒｑｐｈｓｏｎ法代替行列式值试凑法计算频率值简单高效。最后，本给出了两     

变厚度圆柱壳在轴向压力作用下的稳定计算

应用传递矩阵法研究了变厚度圆柱壳的稳定计算问题.用解析法推导出了等厚度圆柱壳单元的精确传递矩阵,再应用传递矩阵原理建立圆柱壳稳定计算的特征方程.该法具有计算简单、节约内存的优点,可用于实际结构的设计.     

旋转薄壁圆柱壳的高节径振动特性以及篦齿结构的影响

本文利用传递矩阵法研究旋转薄壁圆柱壳的静态固有频率和行波特性,特别是其高节径振动特性以及壳体外壁周向篦齿结构（即密封齿）对固有特性的影响．首先基于Love壳体理论,引入科氏力和离心力,建立考虑旋转态的薄壁圆柱壳的动力学控制方程．然后,介绍了进行圆柱壳固有特性分析的传递矩阵法．对比分析分别由传递矩阵法和解析法得到的、三种边界条件（两端简支、两端固支和一端固支一端自由）下的旋转薄壁圆柱壳的高节径模态特性．最后,利用传递矩阵法对带有篦齿结构的旋转薄壁圆柱壳进行了计算,分析了篦齿结构对其高节径振动固有特性的影响．     

求解激波／边界层相互干扰问题的有限体积方法

提出了一种高效的求解二维激波／边界层相互作用问题的算法．应用该方法分析了超音速绕压缩拐角的流动情况,与有关实验数据比较,结果令人满意．本算法采用Ｒｕｎｇｅ－Ｋｕｔａ法对时间进行离散,湍流采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ模型．本算法可以推广到其他激波／边界层相互干扰的问题．     

丁二烯抽提（GPB）装置的模拟计算（三）──萃取精馏塔数学模型与算法

丁二烯抽提（ＧＰＢ）装置的第一、二萃取精馏塔底部均装有侧沸器，实际上是集液板。与一般塔板不同，模拟计算时不能简单地采用三对角矩阵法求解所列方程组，而必须用矩阵求逆法或其它方法求解，这就需要更多的内存和更多的计算时间。文中将数学模型进行变换，使之仍可采用三对角矩阵法求解。在此基础上，应用了二阶Ｒｕｎｇｅ－Ｋｕｔｔａ及三对角矩阵联合法计算，顺利收敛，达到满意的计算精度。     

复合材料锥柱结合壳的自由振动分析

用改进的传递矩阵法分析了变厚度的复合材料锥柱结合壳的自由振动问题，推广了传递矩阵法，导出了每一个壳元的封闭形式的传递矩阵．对于正交对称铺层的层合薄壳，其一阶微分方程可以采用传递矩阵法表示．     

材料力学中压杆临界力计算的矩阵法

本文从压杆的微分方程出发,推得了压杆的传递矩阵,用传递矩阵法对复杂的压杆进行了一系列的稳定计算。计算结果表明,本文公式正确,该方法具有精确度高、力学概念清晰及简便易行的优点,适于工程技术人员在微机上应用     

摄动力对静止卫星轨道根数的影响

给出了静止卫星在地球非球形引力、日月引力和光压作用下,以无奇点轨道根数为变量的摄动微分方程组．根据静止卫星轨道小倾角和小偏心率的特点进行了简化,用平均根数法求出摄动微分方程组的解．利用方程组的解计算各种摄动因素对静止卫星轨道根数的影响的量级估计,并用Ｒｕｎｇｅ－Ｋｕｔａ７（８）阶方法求出卫星运动微分方程组的数值解,通过解析解与数值解的比较,说明解析解的精度为１０－５量级     

用整体传递矩阵法计算旋转机械整机固有频率特性

工程中分析多转子系统固有特性常采用子结构传递矩阵法,由于具体结构的复杂性和多样性,用子结构法计算时处理繁琐,编制通用程序困难.为克服子结构传递矩阵法的缺点,给出了整体传递矩阵法计算过程.取各转子状态向量的集合作为系统的状态向量,将各转子通过耦合单元联接起来,该方法不会引入未知内力和位移,也不必建立分割处的平衡方程或变形协调条件.推导了各向同性耦合单元的传递矩阵,并进行了算例验证.由计算结果可见,整体传递矩阵法计算精度较高,它保持了传递矩阵法编程简单、计算工作量小和运算速度快的特点,为开发通用软件提供了有效的方法.     

传递矩阵法分析计算机床主轴刚度

应用传递矩阵法分析计算机床主轴的刚度 ,该法对梁或轴采用有限元处理 ,动态和静态结合 ,过程简单 ,计算精确     

精密电容器误差的网络并行计算

利用局域网和标准信息传递库构成网络并行计算环境。基于变动边界微扰法,实现了精密电容器误差的并行计算。改进了传统的Ｒｕｎｇｅ－Ｋｕｔｔａ方法。新方法在计算量、计算速度及稳定性上都优于原方法,并在实际计算中取得了良好的效果。最后,将用单机和用并行机进行计算的结果作了对比,讨论了与并行计算效率有关的因素。     

复合材料锥柱结合壳的自由振动分析

用改进的传递矩阵法分析了变厚度的复合材料锥柱结合壳的自由振动问题，推广了传递矩阵法，导出了每一个壳元的封闭形式的传递矩阵，对于正交对称铺层的层薄壳，其一阶微分方程可以传递矩阵法表示。     

受弹性约束的双层圆柱壳的振动分析

本文采用传递矩阵法对受弹性约束的双层圆柱壳进行自由振动和稳态响应研究,并运用复阻尼理论分析材料内阻尼对结构响应的影响。本文所用的分析方法对于承受各种荷载作用而只有任意齐次边界条件的圆柱壳结构都能适用,并可应用于求解多层圆柱壳的振动问题。     

强迫布鲁塞尔振子奇异吸引子的各种维数

用周期采样降低吸引子维数办法，计算和讨论了强迫布鲁塞尔振子的几种典型吸引子的各种维数。讨论了对于ｄｃ和ｄＬ关系的推测。指出Ｒｕｎｇｅ－Ｋｕｔｔａ差分方程使Ｄｑ或ｄｃ不收敛的原因及克服的办法。     

传递矩阵法解连续梁

本文引用传递矩阵法解连续梁。用矩阵形式表达梁的内力和变形,导出各种载荷下的跨间传递矩阵和结点矩阵,并用矩阵形式求解,使结果变得程序化。     

任意变厚度圆板的轴对称静力分析

采用传递矩阵法求解任意变厚度轴对称圆板的内力和位移问题，并导出了在任意荷载作用下环板和圆板单元传递矩阵的精确公式，最后，给出了一些计算结果．     

不同边界条件薄壁截锥壳的高阶振动特性研究

采用传递矩阵法研究了不同边界条件下薄壁截锥壳的高阶振动特性.基于Love壳体理论建立薄壁截锥壳振动微分方程,根据薄壁截锥壳子段间的状态向量,通过传递矩阵法得出整体传递矩阵,并用高精度的精细积分法计算固有频率,通过文献和有限元法进行验证,并分析了薄壁截锥壳在不同边界条件下的高阶振动特性.结果表明,不同边界条件下,采用传递矩阵法计算高阶固有频率与有限元法的计算结果基本一致.当轴向半波数增加时,频率明显增大;随着周向波数的增加,频率先减小后增大.固支-固支和简支-简支边界下在m=1和n=7处得到最小频率,固支-自由边界下在m=1和n=6处得到最小频率,三种边界下最小频率值分别为400.1、325.6和226.1 Hz;边界条件约束越多,最低阶固有频率越大.     

阶跃过渡应答法研究甲醇完全氧化反应动力学──反应机理及动力学参数的求取

在第一部分基础上，进一步研究了甲醇完全氧化反应过程中某些物种的行为，证实了甲醇在ＰｄＡｇ／Ａｌ２Ｏ３催化剂上的吸附是该反应的慢步骤，并提出了反应机理，用ＲｕｎｇｅＫｕｔａ法和修正的ＧａｕｓＮｅｗｔｏｎ法求取了各基元步骤的速率常数、活化能和指前因子等动力学参数。