不规则边界上的近似超奇性边界积分方程

提出角点、棱边处不规则边界上的近似超奇性边界积分方程，其依据是边界上任一Cauchy或Hadamard主值积分可以由2个相应的近奇性边界积分的平均值近似表达。近似方程与光滑边界上的超奇性边界积分方程具有统一的表达形式，并且不需要正则化过程，有利于简化程序设计。给出的数值算例说明了近似方程的可行性与有效性。     

发动机燃烧室主动冷却管道的热-力耦合分析

为了准确高效地进行发动机冷却管道结构的热-力耦合数值分析,基于边界单元法提出一种冷却介质中结构对流换热过程的新型计算方法,计算过程中采用径向积分法将对流换热边界积分方程中的域积分转换为等效的边界积分,从而显著降低计算难度。采用改进后的边界单元法和有限元方法分别进行发动机燃烧室主动冷却管道处的热-力耦合分析,计算并获得了该处的温度场、位移场和应力场,发现了不同物理量随管道轴向的变化规律。通过比较基于两种不同数值方法的计算结果可以发现本文数值方法在显著降低计算模型复杂度的同时,取得了合理的计算结果。因此利用本文方法可以简便有效地进行发动机燃烧室主动冷却管道结构的热-力耦合分析。      

非紧致气动噪声半空间传播基本解的边界积分方法

为了获得低马赫数流动诱发的非紧致气动噪声在半空间内传播的基本解,结合复等效源方法和边界元方法建立了半空间精确格林函数的边界积分方程,当半空间边界为刚度型阻抗边界时可避免奇异性积分.同时基于等效源方法提出一种半空间二维非紧致圆柱声散射模型,推导了静止介质中声散射基本解的理论表达式.对静止介质中的二维圆柱声散射,数值解在研...     

考虑地面影响的半空间气动噪声数值计算方法

含非紧致边界的半空间气动噪声问题普遍存在于航空领域,在数值计算过程中不仅要考虑非紧致边界的散射效应,还需要考虑地面对声波的反射。基于Lighthill声模拟理论和广义格林函数方法,引入边界声压积分项和可压缩修正项可以考虑非紧致边界以及流体可压缩性对声波散射的影响。然而在半空间声场的计算过程中,包含沿无限大地面的积分项,为了克服沿地面积分所带来的巨大工作量,根据无穷大平板对声波的反射特性,推导得出频域下基于半空间格林函数的声学积分方程。同时开展了半空间中应力源噪声与圆柱绕流噪声的数值预测与分析,验证了该方法的准确性与高效性。     

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

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

实现边界元并行分布计算的一种新方法

 给出了实现边界元并行分布计算的一种新方法。它将边界元离散节点分为与计算网络所拥有的处理机台数相对应的若干组, 以此实现任务分配, 独立地计算形成每组节点对应的离散边界元积分方程, 最后利用迭代方法行求解分布的边界元方程组。以弹性问题为例, 在Inmos T800 Transputer 系统上实现了边界元计算。验证了本文方法的可行性与有效性。     

关于“跨音速流的边界积分法”一文的某些问题

 最近杨岞生教授和J.C.Wu教授提出了一种跨音速流的边界积分法。该方法以脉冲压强为基础建立非线性可压缩位流的边界积分表示式,利用边界元法给出了二维、三维机翼跨音速位流的计算方法。我们认为文献[1]在推导该方法的基本方程（即主管脉冲压强的方程）和求此方程的基本解这两个方面存在一些问题。     

无粘跨声速绕流的有限元法和界积分方程的某些进展

本文对无粘跨声速绕流的有限元法和边界积分方程法的目前进展,特别是对提高跨声速流算法效率起着主要作用的成果,例如普适元素和块结构有限元网格生成技术有限元解法的多重网格技术、SUPG有限元法、GRMES算法和可压流的完全边界积分方程法等方面作了回顾和评介。     

跨音速流中的一种边界元技术

边界元的名称源于C.A.Brebbian,边界元技术的特点是节点未知量均位于物面边界上,从而大大简化数值求解过程和降低对计算机容量的要求。本文将边界元技术推广应用于三维机翼的跨音速绕流的非线性问题,建立了一种适用于跨音速范围的带权余数表示式,并将其直接应用于全速位方程。为了保证冲波的不可逆特性,并使解在超临界流区域容易收敛,我们在全速位方程中增加有人工粘性项,而且将流场区域划分成许多小子域,并对每个子域应用格林定理。如果对每个子域我们所选用的插值函数对速位φ具有C~1连续性,则相邻元素交界面的面积分将相互抵消,使剩余的元素面积分(因而节点未知量)均位于机翼边界面上,从而导得边界元的积分方程。     

三维机翼可压缩层湍流边界层计算

 从三维时间相依可压缩边界层动量积分方程和平均流动能积分方程出发,将二维可压缩层、湍流边界层积分方程算法推广到求解有限翼展后掠机翼的边界层流动。利用四步Runge-Kutta时间步进方案数值求解积分方程组,并利用当地时间步长加速迭代收敛。文中分析了数值方法的稳定性与收敛性,并考查了横向耗散项对计算结果的影响。算例表明,能获得令人满意的三维机翼定常可压缩层、湍流边界层的计算结果。     

纤维增强及匀质层层叠板动力分析的解析与半解析方法

本文引入复合板条振型函数系列,提出二种层叠板解法。一种纯解析法适用于各种边界条件的矩形层叠板的动力计算;另一种半解析半数值方法适用于任意形状、任意边界的层叠板动力计算,比现有数值计算方法节省大量内存与时间。文中还给出了一系列常用层叠板计算公式。     

Hoff型夹层板振动分析的边界积分方程法

 本文利用静力问题的基本解建立了Hoff型夹层板振动分析的边界积分方程。在数值实现过程中,边界上和区域内的未知函数均釆用线性插值以提高计算精度。数值结果表明本文方法具有未知量少和计算精度高等优点。对于低阶频率,本文结果与精确解的误差不大于1％。     

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

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

含孔边裂纹板的弯曲断裂计算

 采用复变函数理论和边界配置方法 ,对含孔边裂纹板的弯曲断裂进行了分析计算。首先假设挠度的复变函数式 ,进而可以求板的内力。它们能满足一系列的基本方程和支配条件 ,仅板的边界条件需要考虑 ,并且可用边界配置法和最小二乘法近似满足。对孔双边裂纹问题进行了应力强度因子计算。数值算例表明 ,本文方法精度较高 ,计算量小 ,是一种有效的半解析、半数值计算方法。     

薄板弯曲问题的Trefftz方法

 本文从弹性力学广义变分原理出发,提出了一种求解薄板弯曲问题的Trefftz方法。为了让近似函数满足平衡微分方程,利用复变函数理论研究了试函数空间的构造问题。最后对矩形板给出了3个算例,计算结果表明该方法具有较好的收敛性和稳定性。     

叶片空间造型的计算几何方法

 本文把新近发展起来的计算几何中的几种典型方法（包括Bézier曲面、B样条曲面和Bézier-B样条混合曲面等）,同航空发动机叶片空间造型结合起来,形成一套计算机辅助叶片外形设计的新方法。扭曲叶片由三二次混合曲面片拼接而成,型面有完整的数学模型。因此有利于叶片三元流场性能和有限元强度计算分析以及与数控加工相结合。实例表明该方法简单、灵活、计算量小,对大弯度叶片设计亦行之有效。     

任意坐标系中的弹塑性矩阵

 本文推导了三维问题及平面应力、平面应变、轴对称问题在任意曲线坐标系中的弹塑性矩阵。由此最一般的表达式可以方便地得到任何特定坐标系中的弹塑性矩阵,供各类不同特点问题的有限元计算时使用。 文中推导得到的任意曲线坐标系中的本构方程,也同时为曲线差分法求解弹塑性问题提供了刚度系数显式。     

一种含分层复合材料层合板高阶理论及有限元模型

本文提出了一种新的复合材料县合板高阶弯曲位移模式。此模式考虑了横向剪切的影响,得到一种新的厚板有限元模型。对于含有分层的单元,通过改变其单元内横向剪切应力的分布,达到改变单元刚度的目的。因此,此模型适用于任何厚度的或含有分层的复合材料层合板。     

用矩量法计算任意形状目标雷达散射截面

 1.引言 在电磁散射计算问题中,过去研究较多的是高频情况,所用方法通常是几何绕射法、光学法等。用矩量法计基电磁散射问题,60年代提出了网格法。这种方法多用于尺寸较小的物体,如线天线、小平板等。l984年美国人用矩量法计算了飞机、导弹等的雷达截     

Vibrational behavior of isotropic plate structures in contact with a bounded fluid via unified formulation

This paper presents an analytical solution for free vibration analysis of thick rectangular isotropic plates coupled with a bounded fluid for various boundary conditions. In order to consider displacement theories of an arbitrary order, the Carrera Unified Formulation(CUF) is used. The eigenvalue problem is obtained by using the energy functional, considering plate and fluid kinetic energies as well as the potential energy of the plate. The Ritz method is used to evaluate the displacement variables, and the functions used in the Ritz series can be adjusted to consider any of the classical boundary conditions. The convergence of the solution is analyzed, and a validation of results considering open literature and 3D finite element software is performed. Parametric studies are carried out to obtain natural frequencies as a function of the side-to-thickness ratio, plate aspect ratio, fluid domain size, plate boundary conditions, and fluid-solid density ratio. Pressure and velocity in the fluid domain are evaluated in order to establish the consistency of the solution.Accurate results for thick plates are obtained with a much lower computational cost compared to that of 3D finite element solutions.