调和微分求积法权系数矩阵的一种显式计算式

简要介绍了调和微分求积法，导出了求一阶导数权系数矩阵的显式计算公式。利用该公式和其中反心对称的性能，可进一步提高计算效率。由于均匀网点有时不能给出可靠的解，本文导出了几种能出可靠结果的不等距网点公式，其中一种公式虽然用不同的方法导出，但结果与Ｇａｕｓｓ－Ｌｏｂａｔｔｏ方法等价，本文还证明了调和微分求积法权系数矩阵具有中心对称或中心反对称的性质（取决于导数的阶数），利用这些性质可以进一步减少计算工程     

BUCKLING ANALYSIS OF POLAR ORTHOTROPIC ANNULAR PLATES UNDER UNIFORM PRESSURES

ＢＵＣＫＬＩＮＧＡＮＡＬＹＳＩＳＯＦＰＯＬＡＲＯＲＴＨＯＴＲＯＰＩＣＡＮＮＵＬＡＲＰＬＡＴＥＳＵＮＤＥＲＵＮＩＦＯＲＭＰＲＥＳＳＵＲＥＳＹｏｕＲｅｎｃｈａｎ；ＷａｎｇＹｏｎｇｌｉａｎｇ；ＷａｎｇＸｉｎｗｅｉ（ＤｅｐａｒｔｍｅｎｔｏｆＡｉｒｃｒａｆｔＥ...     

CALCULATION OF THREE DIMENSIONAL BOUNDARY LAYER EQUATIONS USING DIFFERENTIAL QUADRATURE METHOD

ＣＡＬＣＵＬＡＴＩＯＮＯＦＴＨＲＥＥＤＩＭＥＮＳＩＯＮＡＬＢＯＵＮＤＡＲＹＬＡＹＥＲＥＱＵＡＴＩＯＮＳＵＳＩＮＧＤＩＦＦＥＲＥＮＴＩＡＬＱＵＡＤＲＡＴＵＲＥ　ＭＥＴＨＯＤＴａｎｇＤｅｎｇｂｉｎ；ＣｈｅｎｇＧｕｏｗｅｉ（ＤｅｐａｒｔｍｅｎｔｏｆＡｅｒｏ...     

高精度微分求积曲梁单元的建立与应用

首先由能量原理导出曲梁弯曲问题的控制微分方程，在此基础上应用微分求积法原理分别给出了曲梁内点和端点的微分求积方程，由此形成曲梁单元的刚度方程，从而建立了微分求积曲梁单元，并给出了曲梁结构刚度方程的边界条件。通过算例分析，得到了微分求积单元法结构离散时应使单元数量少的原则和求解精度与单元长度比基本无关的性质。与有限元方法的结果比较表明，本文导出的曲梁单元是一种具有很高求解精度的单元。》     

加筋板稳定性微分求积单元法分析

板和加筋板结构的稳定性特性,是设计人员十分关注的一个问题。本文首次采用新近提出的微分求积单元法分析了各向同性加筋板的稳定性问题,建立了微分求积梁单元和板单元,并给出了详细的分析过程。通过与现有结果的对比验证了所建立的分析过程和程序的正确性。计算结果表明:微分求积单元法具有简单、收敛速度快、计算量少和精度高等优点。     

三维非结构网格的生成及其应用

采用阵面推进法生成了高质量的三维四面体非结构网格,该方法的特点是在生成单元的同时引入结点,与其他预布结点的方法相比,其好处是能更好地控制网格单元尺度在整个计算域的合理２并提高网格质量。用Ｄｅｌａｕｎａｙ方法连接存储有控制参数的结点生成了背景网格。网格生成过程中采用了合理的数据结构,有效地提高了程序运行效果,最后给出了Ｍ６机翼的数值算例。     

微分求积单元法分析变曲率井内受径向约束管柱的非线性屈曲

基于Love的空间弯扭杆平衡方程,通过引入井壁约束条件导出了平面变曲率井内受径向约束管柱的平衡方程。采用微分求积单元和增量迭代法求解了曲率线性变化的井内管柱的非线性屈曲问题,通过与有限元计算结果的对比验证了所构建方法和编写程序的正确性,而且还表明微分求积单元法有方法简单、易于实施,计算量少、精度较高等优点。计算结果表明,等曲率井内管柱屈曲的临界载荷明显大于曲率线性变化井内管柱屈曲的临界载荷,另外,变曲率井眼的曲率变化对管柱弯矩、井壁约束力有显著的影响。     

基于一种高精度微分求积法的结构动力响应数值计算

采用了一种高阶精度的Pade逼进的时间域上格栅取点的时间步积分的微分求积方法，对双质点系及梁在强迫力作用下的振动特性进行了数值分析。计算结果表明．这种方法具有明显的高精度及低耗时并且对于二阶初值问题是无条件稳定的。在所考虑的时域内的动力响应过程中，系统的总能量是守恒的。     

离散奇异卷积法结构分析

离散奇异卷积法是近几年内提出的一种新型的、具有潜力的数值计算方法.目前该法已被成功地用于结构元件的力学分析.本文介绍了离散奇异卷积法的原理,尝试用该方法分析工程结构中常用的梁、矩形板结构的弯曲、自由振动性能,并将计算结果与解析法和微分求积法的所得结果进行比较.结果表明,离散奇异卷积法是一种简单、可靠的数值计算方法.     

编队飞行风洞实验研究

在1m非定常风洞中开展了两机编队飞行试验研究。前机采用尾支撑转接垂直叶型支杆与坐标架连接,可以实现相对位置（纵向、侧向和垂向间距）的精确改变;后机通过尾支撑连接到风洞的主支撑机构上,可以实现迎角的变化。采用内式六分量应变天平测量后机的气动力受前机尾涡流影响的变化情况,对后机的绕流场进行了PIV测量。试验中使用了2组模型,一组是简化的翼身组合体模型,另一组是翼身融合体飞翼布局模型。结果表明:当前机翼尖涡靠近后机翼面时,后机的升阻比变化较明显;当前机翼尖涡靠近后机翼尖时,后机可获得最大升阻比;前机迎角增大时,后机的升阻特性有较明显变化;当后机的迎角大于8°时,其升阻比基本不受前机影响。     

多外载联合作用下圆板的非线性弯曲

分析了多外载联合作用下圆板的轴对称非线性弯曲问题。分析从极坐标系下圆板弯曲的Ｖｏｎ－Ｋａｒｍａｎ方程出发，运用微分求积方法（ＤＱ法）导出了控制方程的ＤＱ形式；边缘径向位移和边缘力矩由两个统一的方程来表示，通过改变方程中的约束刚度和边缘载荷系数，实现了对任意边界条件的模拟；对最终得到的非线性方程组，用Ｎｅｗｔｏｎ－Ｒａｐｈｓｏｎ方法进行了迭代求解。文中给出了圆板受横向均布力、板心横向集中力、边缘均布径向力、边缘均布弯矩等四种载荷两两联合作用下的计算结果曲线，讨论了不同联合载荷对回板非线性弯曲的影响。与文献结果比较表明，该方法能满足各种边界条件，具有较高的求解精度。     

有初始横流阵列射流冲击的三维数值计算

对具有初始横流的阵列射流冲击的流场和温度场进行了三维数值计算 ,研究了射流孔与冲击靶面间距、射流冲击孔排列方式和射流与横流质量流量比等因素对换热特性的影响。结果表明 :( 1 )在不同的射流孔与冲击靶面间距下 ,对于给定的射流与横流质量流量比 ,顺排方式的冲击冷却效果要优于叉排方式 ;( 2 )射流孔与冲击靶面间距对阵列射流冲击的流动和换热特性影响是非常复杂的 ;( 3 )随着初始横流与射流质量流量之比增加 ,在相同的射流孔与冲击靶面间距下 ,顺排和叉排方式的冷却效果均呈现单调的下降     

压电圆浅球壳非线性微分求积单元法分析(英文)

用微分求积单元法分析了压电圆浅球壳在外加电压和外载作用下的非线性静力特性 ,首次给出了详细的公式和求解过程 ,并分析了几个典型算例 ,得到了非常精确的结果。基于本文的研究结果可以得出以下的几点结论 :微分求积单元法是一种有用的数值分析方法 ;压电圆浅球壳的几何形状理论上可以被控制 ;对于某些几何形状的压电圆浅球壳 ,当外加电压达到临界值时 ,即使没有外加载荷跳跃失稳也会发生。     

非结构网格生成技术

描述了非结构网格生成的过程,网格的生成方法基于Delaunay 三角形方法,并对物面附近网格进行拉伸以便获得在边界层内高拉伸比的三角形单元。网格的生成过程分两部分:(1) 按Delaunay 三角形方法生成适于无粘情况下的非结构网格。(2) 对近物面的部分三角形单元进行拉伸,使生成的网格能适用于粘流。使用上述方法本文构造了二维弯管内绕叶栅区域的非结构网格     

Experimental Investigation on Convective Heat Transfer Characteristics in Ribbed Channel

The heat transfer and pressure loss characteristics on a square channel with two opposite surfaces roughened by high blockage ratio ribs are measured by systematic experiments.Reynolds numbers studied in the channel range from 1 400 to 8 000.The ratios of rib height to hydraulic diameter (e/D) are 0.2and 0.33,respectively.The ratio of rib spacing to height (P/e) ranges from 5to 15.The rib orientations in the opposite surfaces are symmetrical and staggered arrangements.The results show that the heat transfer coefficients are increased with the increase of rib height and Reynolds number,though at the cost of higher pressure losses.When the rib spacing to height ratio is 10,it keeps the highest heat transfer coefficient in three kinds of rib spacing to height ratios 5,10 and 15.The heat transfer coefficient of symmetrical arrangement ribs is higher than that of the staggered arrangement ribs,but the pressure loss of the symmetrical arrangement ribs is larger than that of the staggered arrangement ribs.     

轴对称进气道外罩跨音速非守恒型全位势流AF3算法

本文针对轴对称进气道外罩跨音速绕流问题的特点,采用了贴体坐标网格;并对这类网格拓扑构造了三种Baker的隐式近似因式分解AF3有限差分迭代算法。通过理论分析与数值实验,找到了既能计算亚临界情形,又能计算超临界情形的AF3格式;并就此格式讨论了适合于进气道绕流特性的AF3迭代格式中间变量的提法,具体提出了格式实施的相应方法,研制了可供分析任意轴对称进气道外罩跨音速特性的计算程序。由于采用了任意曲线坐标系上的轴对称非守恒型全位势方程适合于用不同方法生成的计算网格,因此,本文的工作还有可能进一步推广到计算更复杂的轴对称进气道跨音速绕流问题。计算表明,本方法收敛快,结果好。     

高阻塞比肋化通道对流换热特性实验研究

为了降低涡轮叶片温度以保证其安全性,采用实验方法对高阻塞比肋化通道内流动换热特性进行了测量。针对顺排和叉排两种肋化通道考察了多个主要影响因素对换热特性的影响规律,研究表明:(1)对于高阻塞比肋化通道,当流体在层流区域流动时,肋片对流体在通道内部对流换热的强化作用要明显高于湍流状态;(2)无论是顺排还是叉排肋化通道,当肋间距比为10时,其对应着最大的对流换热系数和最大的流阻损失;(3)在实验几何参数范围内,顺排肋化通道的对流换热系数和流阻系数均高于叉排通道;(4)在高阻塞比肋化通道中,流阻特性的变化规律和流动状态的转变规律基本与水力光滑通道相同。     

一种基于块结构的非结构网格快速生成方法

以阵面推进法为基础,提出了一种用预制模块的方法实现在平面域内生成非结构网格的新方法,即模块方法。模块是由一组三角形构成,文中以三角形模块为例,介绍了网格模块的预制方法。网格的生成是根据空间步长的背景信息,从预存的网格模块中选取当地合适的模块,直接复制或变形到计算区域实现的。由于生成过程中每一步引入的是预制的网格模块而不是传统的单个三角形单元,因此可期望新方法要比传统方法快得多。算例显示了新方法的快速网格生成过程,该方法具有推广应用前景。     

基于移动代理的双重认证:一种高效的网格认证体系结构

网格技术作为下一代Internet技术,已逐渐成为计算机领域研究的重心之一。作为网格安全的重要组成部分,安全认证是一重要环节。因此,在客观分析了GSI（Globus security infrastructure）认证结构的基础上,针对其存在的问题与不足,采用了移动代理技术进行改进,提出了一种基于移动代理的双重（Two-step mobile agent based,TSMAB）网格认证体系结构。该结构通过移动代理技术实现了初级认证和高级认证两个阶段。通过对模型的设计与分析得出：TSMAB认证体系结构提高了原有GSI认证的效率。     

自适应计算网格的生成和优化方法

本文提出一种简单通用的自适应计算网格生成方法。配合文[1]的网格优化方法,能在生成自适应网格的同时,保证网格具有较好的光滑性和正交性。本文还阐述了单元加权函数W的取值和边界点的确定对自适应网格性质的影响。作为应用实例,本文将自适应网格用于二维无粘跨音速叶栅流场计算。结果表明,采用自适应网格能明显改善流场的计算精度。