一维微裂纹的统计演化

 <正> 1.引言 由于组成统计微元的单个微裂纹,既不是古典的牛顿粒子,也不是量子的微观粒子,所以它们既不遵从古典力学的哈密顿方程,又不遵从量子力学的薛定谔方程,它们的运动规律一般由实验或先验模型来给出。本文从简化的、但抓住了物理实质的模型入手,给出了两种成核率函数,应用统计方法,通过计算分析,得出了一维微裂纹统计演化行为的物理图像。     

基于有限元/边界元耦合方法的管道进口声传播及声辐射模型

发展了一种基于有限元和边界元耦合方法的管道进口声传播及声辐射计算模型.该模型将整个声场分为内部有限域和外部无界域,分别用有限元和边界元方法求解控制方程,在两者之间的界面上使用具有物理意义的声阻抗参数进行匹配,并通过一种快速迭代方法实现全声场求解.这种迭代方法可以保证有限元刚度矩阵等带宽以及对称的特性不被破坏,有助于提高计算效率.该模型先得到了Levine-Schwinger标准解的检验,进而在无流动情况下对于简化的航发短舱进口管道模型进行了噪声辐射现象的数值模拟,最后基于计算结果分析了声衬对远场声辐射的影响.      

多转子系统振动分析的子结构耦合矩阵法

提出了一种用于多转子耦合系统振动分析的建模方法.将耦合的多转子系统在耦合处分开,把耦合力作为单转子系统上的外力,利用集总参数法和节点的自由振动微分方程通式建立单转子系统的振动方程.通过分析耦合部件作用在耦合节点上的力,得到了多转子系统耦合部件的耦合矩阵;利用耦合矩阵将单转子系统振动微分方程耦合在一起,建立了整个多转子系统的振动方程.通过分析多转子系统的振动方程,解出了多转子系统的固有频率、临界转速、幅频响应曲线、涡动轨迹等.通过两个算例,证实了多转子系统振动分析的子结构耦合矩阵法是可行的,计算结果是可靠的.      

一种简化的机载MIMO雷达杂波特征相消器

针对机载多输入多输出(MIMO)雷达杂波抑制问题,提出一种简化的杂波特征相消器(EC).根据杂波在空时二维平面的先验分布离线构造杂波子空间.以此替代由协方差矩阵特征值分解(EVD)得到的杂波子空间,从而将最优权简化为一个确知投影矩阵与目标信号空时二维导向矢量的乘积,避免了传统EC方法中复杂的协方差矩阵估计和EVD运算,...     

实数域中的投影法及其在单次快摄阵处理中的应用

阵列处理中计算简单的投影法避免了矩阵特征分解(ＥＶＤ)的庞大计算量,然而仍需进行复数计算。为了进一步减小计算量,提出一种实数域中的投影法,它将投影矩阵及空间谱的计算均简化为实数的计算,从而。使阵处理的计算量大大减少,而且使估计性能保持不变。计算机模拟与实际外场实验数据的计算结果表明了算法的有效性。     

一种基于有限元/边界元耦合方法的管道进口声传播及声辐射模型

本文发展了一种基于有限元和边界元耦合方法的管道进口声传播及声辐射计算模型。该模型将整个声场分为内部有限域和外部无界域,分别用有限元和边界元方法求解控制方程,在二者之间的界面上使用具有物理意义的声阻抗参数进行匹配,并通过一种快速迭代方法实现全声场求解。这种迭代方法可以保证有限元刚度矩阵等带宽以及对称的特性不被破坏,有助于提高计算效率。该模型先得到了Levine and Schwinger标准解的检验,进而在无流动情况下对于简化的航发短舱进口管道模型进行了噪声辐射现象的数值模拟,最后基于计算结果分析了声衬对远场声辐射的影响。     

双曲面相关定位的两种解析解法

本文研究了四站距离差(时间差)的空间定位问题,给出了空间定位的两种求解析解方法。一种由几何原理导出,另一种则通过定位方程的演变得到;两方法最终都将双曲面定位问题简化为一个二次方程求解。由此还可弄清距离差相关定位一些基本特性。     

大小叶盘结构连续参数模型和振动模态

给出了用于研究大小叶片整体叶盘结构固有振动特性分析的连续参数模型。模型用 Timoshinko梁模拟叶片 ,用平板模拟轮盘 ,并在叶片和轮盘间引入弹簧元件模拟叶片和轮盘间的耦合关系 ,简化了两者间的边界条件和连续性条件。采用 Gram-Schmidt方法生成正交多项式作为李兹容许函数进行模型的离散化 ,简化了所得到的频率方程求解过程 ,有利于求解过程的数值稳定性 ,并可方便地通过对方程系数性质的分析了解结构振动模态的性质。利用上述的分析原理和求解方法 ,进行了简化的大小叶片整体叶盘结构的计算 ,计算结果与有限元计算结果具有很好的一致性      

空间近程高速相对运动方程的状态转移矩阵

分析了空间近程高速相对运动的特点,建立了描述近程高速相对运动的线性微分方程组,并给出了状态转移矩阵与解析解。仿真结果表明:对于HEO-GEO近程高速相对运动,所建方程的相对距离计算误差仅为16 cm,不足Hill方程计算误差的15%。所给出的状态转移矩阵对于实现空间目标拦截等任务的实时高精度相对导航和制导算法具有一定的应用价值。     

阻尼振动系统复特征解的摄动分析

当振动系统的质量矩阵和刚度矩阵作小修改时,不重复求解大型广义特征问题,在初始的模态坐标变换基础上,用矩阵摄动法给出修改系统的振动频率和振型;然后利用修改系统的振动频率与比例阻尼参数求出修改系统的阻尼比。工程算例表明,这种方法计算简单,具有较高的精度。     

并行计算中一种非结构网格分割方法

将递归谱对剖分方法应用于流体力学并行计算中的非结构网格分割,以解决负载平衡和最小切割问题。为使用这种方法,计算了网格伴随图的离散Laplacian矩阵的第二特征矢,然后从该特征矢的分量引入网格的对剖分。特征矢计算中应用了Rayleigh商迭代,并进行了一些修正以使收敛强烈地偏向于第二特征矢及考虑逆迭代步中线性方程组的迭代求解。最后,通过非结构自适应网格上Euler方程分区计算的数值结果验证了所发展的网格分割方法。     

连续最小阶奇异滤波器

 本文讨论了当动态噪声统计特性未知时,奇异线性定常连续随机系统最小阶滤波器的设计问题。在系统部分观测量能精确测量的情况下,利用广义逆阵方法选择L矩阵,以消除动态噪声对降阶系统的影响,从而推导出缺动态噪声统计特性时的降阶奇异最优滤波器,其阶数为n-m+r。当量测方程奇异假设条件成立及引理有解时,本文证明了最小阶奇异滤波器必定存在。文中举例说明了这一降阶滤波方法的可行性。     

计算特征向量导数的若干方法的短评和一种改进模态法

 对特征向量导数计算的若干方法作了简短地评述,同时提出了一种改进模态法。这种方法具备Nelson法、改进Nelson法和直接扰动法的优点。该方法还可以具有多种形式,但本文仅介绍了一种同于直接扰动法的主要形式,这也是精度最好的一种形式。     

一种求解欧拉方程改进的通量分裂方法

本文以Steger-Warming提出求解Euler方程的通量分裂方法为基础,利用改进的特征值分裂方法,建立了一种连续的通量分裂方法。该方法克服了Steger-Warming通量分裂法存在的分裂后的通量在个别点导数不连续的问题,同时方法简单,便于推广。对于通量分裂后的Euler方程,采用迎风格式,建立了相应的二阶隐式有限体积方程。为了减少求解所得块矩阵方程组所需要的工作量,本文利用特征变换,在不改变原方程组离散化精度的前提下,将所得块矩阵方程组化为分离的代数方程组进行求解,从而大大减少了计算量。文中分别对一维管流,二维绕NACA-0012和NACA64A-10翼型在亚声速和超声速来流的跨声速流场进行了计算,并将所得结果与有关实验和计算结果进行了比较,均符合良好。 本文所建立的方法具有占机内存小,计算时间短的优点。在普通中小型计算机上亦可使用。     

两端支承输流管道内共振条件分析

讨论了两端支承输流管道发生内共振的条件.考虑管道变形的非线性因素,建立了两端简支的输流管道,当简支端简谐运动时的非线性运动方程.通过对其线性部分系数矩阵的分析,发现系统发生几种内共振时流速uo与质量比mi的关系.从理论上解释了流速诱发两端支承输流管道内共振的机理.     

两端固定输流管道在脉动内流作用下参数共振的最优控制

采用压电片作为控制执行元件,对两端固定输流管道在脉动流作用下的参数共振进行了最优控制.基于二次型性能指标最优控制原理,针对离散后的运动方程在零平衡点附近Jacobian矩阵的非时变部分设计了最优控制器.最后,通过数值仿真的方法对控制方案的有效性进行了验证.结果表明,该方案对输流管道在第一、第二振型次谐波共振区域及组合共振区域的参数共振的控制均有很好的效果.     

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

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

红外窗口不同冷却方式下的结构传热和热应力特性计算研究

对外冷喷流和内冷管流等不同冷却方式下的红外窗口传热和热应力进行了数值计算研究.窗口结构传热和热应力采用有限元方法求解,带外部冷喷流窗口外表面气动加热率则通过对带冷源项的N-S方程求解给出,当采用管流冷却时,管流液体温度和窗口结构温度通过耦合迭代方式统一求解,管流和壁面间的换热系数采用工程关联公式估算.研究表明,在达到来流总温以前,窗口各点温度和热应力随加热时间单调上升,各时刻最大温度发生在外表面;而最大热应力则发生在合金材料内部.两种冷却方式对比表明,外冷方式对于窗口整体温度和热应力的降低十分有效,但具体到局部重要部位,外冷方式效果不如内冷,内冷方式对于管道附近部位具有更好的降温和降低热应力效果.因红外窗口尺度限制,冷却管道流动雷诺数偏小,流动为层流态,这限制了冷却管换热效率的提高,因此建议增加管道数目和管壁粗糙度来强化冷却.     

航空发动机供油管裂纹失效分析

针对某航空发动机供油管钎焊焊缝边缘出现裂纹而导致的供油管渗油故障,对其断口及源区附近表面进行宏、微观检查,明确了供油管裂纹性质为高周疲劳裂纹,疲劳源区呈沿晶开裂特征。利用能谱分析手段发现裂纹源区存在外来Cu元素,通过对焊缝及开裂部位进行金相分析,发现钎焊料中的Cu元素沿基体晶界发生扩散,导致晶界粗化,在拉应力作用下,Cu与基体形成的脆性相开裂形成沿晶裂纹。晶界粗化和沿晶裂纹是供油管过早发生疲劳开裂的主要原因。建议将Cu基钎料更换为Ag焊料。     

Static/dynamic Analysis of Functionally Graded and Layered Magneto-electro-elastic Plate/pipe under Hamiltonian System

The 3-dimensional couple equations of magneto-electro-elastic structures are derived under Hamiltonian system based on the Hamilton principle. The problem of single sort of variables is converted into the problem of double sorts of variables, and the Hamilton canonical equations are established. The 3-dimensional problem of magneto-electro-elastic structure which is investigated in Euclidean space commonly is converted into symplectic system. At the same time the Lagrange system is converted into Hamiltonian system. As an example, the dynamic characteristics of the simply supported functionally graded magneto-electro-elastic material （FGMM） plate and pipe are investigated. Finally, the problem is solved by symplectic algorithm. The results show that the physical quantities of displacement, electric potential and magnetic potential etc. change continuously at the interfaces between layers under the transverse pressure while some other physical quantities such as the stress, electric and magnetic displacement are not continuous. The dynamic stiffness is increased by the piezoelectric effect while decreased by the piezomagnetic effect.