Lyapunov矩阵方程的一种数值解法及部分并行处理

本文设计了求解Lyapunov矩阵方程的一种新方法。所考虑的矩阵方程是 AX—XB=C(1)其中A,B,C分别是m×m,n×n和m×n的已知矩阵。 该方法首先是将系数矩阵A,B初等相似约化为三对角矩阵,即存在可逆矩阵U,V,使U~(-1)AU=A,V~(-1)BV=B,其中A,B为三对角矩阵。然后设计了矩阵方程AY—YB=C的公式解法,分三步: 1)求f(λ)=det(λI—A)的λ各次幂的系数a_0,…,a_m; 2)计算sum from i=1 to m (A_(m-i)-CB~(m-i)),f(B); 3)求解Y。解方程AY—YB=C的方法称为THR算法。 最后经逆变换获得原矩阵方程(1)的解X。 求解矩阵方程(1)的方法称为R—THR算法。该方法的计算量约为m~3+4/3n~3+7m~2n+5nm~2+m~2。 本文给出了R—THR的串行计算的数值例子,并给出了THR算法的并行计算格式。最后通过几种数值方法的比较,表明该方法是可行的,也是有效的。     

舱门机构运动的计算机模拟

1.构造及其原始数据 滑轨舱门收放机构由中、侧滑轨、三角支架、外伸臂、链条条、链轮、液压马达及减速装置等组成(图1)。一根中滑轨布置在43框前飞机对称平面内与飞机外形的下零纵线相贴合,两根侧滑轨分别布置在43框后机身的两侧,舱门的前端由一个接头铰支在三角支架上,三角支架的前支点与链条相连,舱门的后部和两个外伸臂固接,     

Semi-analytic Solution of Steady Temperature Fields During Thin Plate Welding

Usually, it is very difficult to find out an analytical solution to thermal conduction problems during high temperature welding. Therefore, as an important numerical approach, the method of lines （MOLs） is introduced to solve the temperature field characterized by high gradients. The basic idea of the method is to semi-discretize the governing equation of the problem into a system of ordinary differential equations （ODEs） defined on discrete lines by means of the finite difference method, by which the thermal boundary condition with high gradients are directly embodied in formulation. Thus the temperature field can be obtained by solving the ODEs. As a numerical example, the variation of an axisymmetrical temperature field along the plate thickness can be obtained.     

Hamilton正则方程半解析法的收敛和对称分析

近年来,Hamilton正则方程半解析法在工程问题上的应用越来越广泛,但至今未见有关这种方法收敛性和对称性问题研究的文献。基于Hamilton正则方程的半解析法理论,通过变分原理详细推导了Hamiltonian元素的固支和简支边界公式及对称边界公式。多个实例的数值研究表明：随着网格加密,Hamilton正则方程半解析法的收敛速度快于一般传统位移有限元法,对称解法的效率明显优于整体解法。     

C/SiC-W复合材料的制备与表征

采用化学气相渗透(CVI)结合溶液浸渍法制备C/SiC-W多元基复合材料,利用XRD和SEM对材料的物相组成及微观结构进行表征,并对制备过程的反应机理进行初步探讨.结果表明,W在复合材料中成颗粒团聚状,它不仅能渗入到复合材料的纤维束间,还渗入到纤维束内.制得C/SiC-W复合材料的密度为4.1g/cm3,开气孔率为12%.     

冲击响应谱绝对校准研究中冲击响应谱的计算方法

介绍了应用于冲击响应谱绝对校准中的冲击响应谱的计算方法,该方法的实现基于MATLAB编程语言,即通过MATLAB中的仿真直接求解单自由度二阶微分方程的数值解以获得冲击响应谱的绝对复现.     

CE/SE法模拟爆震波点火及传播过程

运用时空守恒元和解元(CE/SE)法模拟一维爆震波的形成及传播过程,分析产生爆震波的不同机理和DDT(爆燃-爆震转捩)过程的影响因素.采取激波点火和热点点火两种点火方式,分别对应SDT(激波-爆震转捩)和DDT过程,热点点火初始场采用线性温度分布和阶梯性温度分布,整个管道均为常压.用隐式梯形方法实现刚性源项积分.结果表明,SDT和DDT对应于不同的爆震波机理,并且DDT受点火区不同参数的影响,为研究爆震波点火和DDT过程提供了理论依据.     

层板孔隙结构中两个任意尺寸半球形珠状发汗的蒸发与燃烧速率分析

针对微重力条件下层板孔隙结构中两个相邻任意尺寸半球形珠状发汗的蒸发与燃烧过程,建立了相应的数学物理模型;在双球坐标系中求解分子扩散传质偏微分方程,得到发汗液珠的分子扩散组分分布的解析解和扩散速率表达式;分析得到在常温、高温条件下蒸发以及燃烧状态下发汗液珠的相互作用因子式,计算了相互作用的发汗液珠在不同条件下的发汗蒸发量以及变化趋势.     

助推-滑翔式弹道中段的近似解

研究了助推-滑翔式弹道中段的近似解。首先根据不同假设条件,给出以下四种情况的近似解:(1)自由飞行段的解;(2)只考虑气动力常迎角飞行段的近似解;(3)考虑引力和气动力的合成解;(4)平衡滑翔解。再由上述近似解合成助推-滑翔弹道中段的分段近似解。通过与数值解的对比,验证了近似解的合理性。研究结果可为新型助推-滑翔式飞行器在方案论证阶段时的弹道特性分析以及射程、峰值热流和过载等特征参数的估算提供依据。     

一个有理递归序列解的性质

文中分析了方程xn＋1=a＋bxn^2/1＋xn-2^2（其中a,6∈[0,∞）,a＋b〉0）的解的性质．并得山①当a=0,0〈b〈2时,方程的每一个解最终严格递减到0;②当b=0,0〈a时,方程的每一个非平凡解关于方程的平衡解振动。特别,当b=0、0〈a≤2时,方程的每一个正解收敛到平衡解x^-上.