三维瞬态导热温度场边界元计算及一阶数学奇异的处理

推导了三维瞬态导热温度场边界元解法的积分方程和离散代数方程，并编制了相应的计算机程序，采用坐标变换的解析消去法，处理边界元的一阶数学奇异点，从而提高了算法的数值精度。用本程序计算了正方体和圆柱环瞬态导热温度场，数值解与解析解的吻合程度令人满意。     

CFD空间精度分析方法及4种典型畸形网格中WCNS格式精度测试

计算流体力学中的精度阶数体现的是数值解随着网格加密向精确解收敛的速度,计算格式、网格质量和边界处理等都会影响计算精度。5阶非线性加权紧致格式(WCNS-E-5)是一种具有较高的精度和分辨率的有限差分格式,对强激波和复杂网格的适应性强。但是,格式的模板自适应选择器和网格的复杂性等都会对最终计算精度造成影响。为了分析精度,首先给出了两种通过数值解反推计算格式精度的方法:一种方法适合已知精确解的情况,另一种方法适合未知精确解的情况。然后在斜交网格、拉伸网格、拐折网格和扭曲网格等4种典型畸形网格中考察了WCNS-E-5格式的精度,考察算例涉及线性波动问题、等熵涡传播问题和超声速边界层问题等。计算结果表明WCNS-E-5在这四种网格中都体现出高阶精度。     

方腔流动气动噪声边界散射数值预测方法研究北大核心CSCD

方腔流动噪声问题因其流场变化复杂且剧烈而倍受关注,本文主要开展高亚声速方腔流动气动噪声数值预测方法研究。基于任意边界条件的格林函数解和Lighthill声模拟理论,提出可以考虑空间边界影响的气动噪声积分计算方法。数值模拟包含流动和噪声计算两部分,通过二阶精度的DDES模型进行流动数值模拟,边界积分方法计算散射声场分布。数值结果显示声场分布随时间呈现周期性变化,与流场的脉动及其脉动周期一致。观察点的声压级随频率逐渐下降且在谐波频率突然增大。本文计算结果与高精度计算气动声学方法计算结果相符,表明该方法合理、可靠,并且具有较高的计算效率。     

方腔流动气动噪声边界散射数值预测方法研究

方腔流动噪声问题因其流场变化复杂且剧烈而倍受关注,本文主要开展高亚声速方腔流动气动噪声数值预测方法研究。基于任意边界条件的格林函数解和Lighthill声模拟理论,提出可以考虑空间边界影响的气动噪声积分计算方法。数值模拟包含流动和噪声计算两部分,通过二阶精度的DDES模型进行流动数值模拟,边界积分方法计算散射声场分布。数值结果显示声场分布随时间呈现周期性变化,与流场的脉动及其脉动周期一致。观察点的声压级随频率逐渐下降且在谐波频率突然增大。本文计算结果与高精度计算气动声学方法计算结果相符,表明该方法合理、可靠,并且具有较高的计算效率。     

旋转变厚度圆盘的样条积分方程法

一、前言 对于旋转变厚度圆盘的应力分析,一般可用有限元法或差分法,但相对于边界元来说精度较差,若要提高精度就必须将单元分细,随之而来的是方程阶数的增加和计算时间的增长。为了克服上述缺点,本文提出了一种新的样条积分方程法,直接利用弹性力学二维问题的环基本解建立积分方程并用样条函数进行求解。这样只需划分少量单元,就可获得精度很高的数值解,可大大节约计算机内存和计算时间。     

固液相变问题的边界元法数值解

相变问题的数值计算方法有多种,本文用边界元法解应用工程中常见的相变导热问题。     

一种包含运动边界的高精度流场数值计算方法

为了准确快速地模拟运动边界的流场,提出一种反馈力源形式的包含运动边界的非定常流场数值计算方法。该方法采用完全正交的网格,以反馈力源作用点的运动来模拟边界的运动。采用物理量及其各阶导数在边界两侧的突跃修正中心差分格式,使之达到二阶精度,以此离散求解二维不可压Navier-Stokes方程。并且提出了与运动边界相适应的反馈力源构造方法及对边界上速度进行插值的方法。基于此数值计算方法,对低雷诺数的圆柱绕流、静止流体中的振荡圆柱以及昆虫振翅运动的二维非定常流场进行了数值计算,计算结果与以往的数值及实验结果非常吻合,表明本文方法与Peskin的浸入式边界方法在处理运动边界问题时具有同样的高效率,且精度高于浸入式边界方法。     

关于跨音速叶栅物面边界条件的数值试验

本文讨论了跨音速叶栅物面边界的数值处理方法,目的是寻求与我们所用内点格式相匹配的边界上流动参数的算法,以提高计算精度。本文首先推导了适合于物面边界计算的特征关系式,其次总结了五种边界处理方法,进行了数值试验。计算结果表明,配合我们的内点计算格式,以y向抛物外推的Denton修正法计算效果较好,能够较精确地反映跨音速叶栅带激波流场的物理特点。     

隐式通量分裂线松驰迭代法求解平面叶栅无粘绕流问题

本文工作是在 Mac Cormack通量分裂格式 [2 ] 基础上发展的一种有限面积通量分裂隐式格式 ,其特点是在求解隐式离散化方程时 ,采用往返扫描一次的 Gauss- Seidel线松驰迭代方法 ,避免了对时间步长的任何限制。为提高定常解精度 ,格式的显式右端项采用二阶精度的离散。在数值求解跨音速涡轮平面叶栅问题中 ,对壁面边界作了较仔细的隐式处理。数值计算表明本方法保持了 Mac Cormack格式具有的高收敛速率 ,(约 30个时间步即可达到定常解 )而每一时间步计算量约减少一半。数值结果与实验结果符合得很好。     

基于高阶间断伽辽金数值方法的浸没边界法

为了降低复杂边界条件下贴体网格生成难度，近年来基于笛卡尔网格的浸没边界法逐渐成为研究此类问题的主要数值方法之一，然而计算精度和计算效率仍然是此类方法目前面临的挑战。与传统的空间2阶精度的有限体积格式相比，空间精度为3阶或3阶以上的高精度方法具有空间精度高、数值分辨率高、数值耗散小的优点，而作为高精度数值方法之一的间断伽辽金方法在浸没边界方面的应用仍较少。本文将高阶间断伽辽金方法的高精度优势与浸没边界法无需贴体网格的优势结合起来，提出了适用于可压缩流动的高精度浸没边界法。其中边界条件采用体积惩罚方式实现，同时采用牛顿法迭代以及MPI并行提高计算效率，物面的数据重建采用插值点处的逆距离权重（IDW-IP）方法代替高阶格式下多项式插值方法。本文基于笛卡尔网格测试验证了二维定常和非定常情况下的数值模拟效果，并与传统贴体网格的计算结果进行了对比。     

An advanced five-unknown higher-order theory for free vibration of composite and sandwich plates

Accurate prediction of dynamic characteristics is quite critical to understand the strength of layered structures. Nevertheless, the existing five-unknown higher-order theories encounter difficulties to forecast accurately the dynamic response of sandwich structures. Therefore, a new five-unknown higher-order theory is developed for free vibration analysis of composite and sandwich plates, which possesses the same degree of freedom as those of other five-unknown higher-order theories. The developed model can meet beforehand interlaminar continuity conditions and the free-surface conditions of transverse shear stresses. To assess capability of the proposed model, analytical solution for such composite structures with simply-supported conditions has been presented by employing Hamilton’s principle, which is utilized for analysis of mechanical behaviors of composite and sandwich plates. Compared with the three-dimensional (3D) elasticity solutions, 3D finite element results and the results obtained from the chosen five-unknown higher-order models, the proposed model can yield accurately natural frequencies of composite and sandwich plates. Even for the thick plates, the higher-order frequencies calculated from the proposed model are in good agreement with the 3D finite element results. By studying effect of the thickness/length ratios on natural frequencies, it is found that the proposed model is adaptable to predicting natural frequencies of the sandwich plates with the thickness/length ratios between 1/4 and 1/100. In addition, some factors influencing accuracy of five-unknown higher-order models have been investigated in detail. Finally, by means of numerical analysis and discussion, some conclusions have been drawn as well, which can serve as a reference for other investigators.     

蜂窝夹层板动力分析的Hamilton体系与辛算法

将蜂窝夹层板动力分析从Lagrange体系改换为Hamilton体系。通过作者早已提出的一条简单而统一的新途径,建立了蜂窝夹层板动力学的相空间非传统Hamilton变分原理,并从该原理推导出相应的Hamilton正则方程、边界条件与初始条件。然后基于这种相空间非传统Hamilton变分原理,提出蜂窝夹层板动力响应分析的辛空间有限元-时间子域法,算例的数值结果表明,这种新方法的计算精度与效率都明显高于国际上常用的Wilson-θ法和Newmark-β法。     

基于三角剪切变形理论的对称层合板振动分析

根据三角剪切变形层合梁理论,推导出对称复合材料层合板的运动微分方程,利用逆复合二次径向基函数无网格配点法对运动微分方程进行离散,预测了对称复合材料层合板的自由振动特性。将不同材料参数、几何尺寸的层合板固有频率计算结果与相关文献中的结果进行对比,结果表明：用逆复合二次径向基函数离散的三角剪切变形理论在对称复合材料层合板自由振动分析方面具有一定的精度。     

Vibrations of circular sandwich plates connected with an elastic foundation

The results of analyzing axisymmetric forced vibrations of an elastic circular sandwich plate connected with an elastic foundation under different surface loads are presented. The hypotheses of a broken normal are accepted to describe the kinematics of a package that is nonsymmetrical over the thickness. The analytical solutions of a problem for plates with a light filler are obtained and numerical analysis is carried out.     

Calculation of the transonic dip of airfoils using viscous-inviscid aerodynamic interaction method

A numerical procedure for the calculation of the transonic dip of airfoils in the time domain is presented. A viscous-inviscid aerodynamic interaction method is taken to calculate the unsteady aerodynamic loads. In the present case the integral boundary layer equations are coupled with the Transonic Small Disturbance (TSD) Potential Equation. The coupling between structural motion and aerodynamic loads is carried out using State Space equation. It is solved by State Transition Matrix technique. Results are presented for NACA 64A010 and NLR 7301 airfoils with structural data from Isogai and DLR, respectively. Comparisons show good agreement with other numerical results. Certain deviations of experimental data taken from literature need more insight in the detailed test conditions.     

边界元法分析回转体的自由扭振问题

 边界元法用于分析非线性问题和动力问题时,在最终得到的积分方程中同时包含边界积分和区域积分。通常的做法是对区域积分项采用常值有限元离散,这需要用很密的网格才能达到较高的精度。     

航空板结构强度的边界元法

边界元即边界积分方程方法,可方便地用来解板弯曲问题。本文在局部极坐标下建立公式,采用具有较高协调性的函数插值方案,利用外点技术完全避免边界奇异积分。列举的一系列算例表明该方法输入数据少、机时省、精度高,是板结构强度分析的一个有效方法。     

面内弯曲载荷作用下两边简支两边固支复合材料层合板的屈曲

为扩充航空复合材料结构稳定性设计手段,基于Levy形式级数和有限差分法,提出了面内纯弯曲载荷作用下两边简支两边固支（SSCC）复合材料层合板临界屈曲载荷的数值解法。首先,应用有限差分法将挠度函数进行离散化处理;然后利用边界条件扩充方程组,最后将临界屈曲载荷的求解转化为广义特征值问题。通过数值算例分别在各向同性板和复合材料层合板上验证了本文方法,并与有限元方法（FEM）进行了比较。结果表明,本文方法具有很高的精度,为复合材料层合板的屈曲分析开辟了新的途径。     

论定常不可压粘流的积分方程解法

本文对定常不可压粘性绕流积分解法中的涡量迭代过程、涡量边界条件的提法以及压强系数的计算等方面都作了改进,使数值解更易于稳定收敛并提高精度。为了验证本文的方法,对低雷诺数下的圆柱绕流进行了计算,结果表明计算效率是高的,迭代过程是稳定的。将本文方法推广到三维流动没有原则上的困难。     

MLPG混合配点法分析复合材料层合板的弯曲问题

提出了分析复合材料层合板弯曲问题的无网格局部彼得洛夫-伽辽金混合配点法,层合板理论采用经典理论。通过移动最小二乘近似构造形函数,采用配点法建立系统平衡方程,推导了系统刚度矩阵和载荷力向量。采用配点法和罚函数法分别施加自然边界条件和本质边界条件。计算实例表明,MLPG混合配点法具有较高的计算精度,并且相比其它的无网格法具有更高的计算效率。