三维移动脉动源格林函数的对称性及其应用

为进一步提高三维有航速船舶耐波性预报的计算效率,将对称性引入到移动脉动源格林函数法的速度势求解之中,提出了基于格林函数对称性的水动力快速求解措施。从源点对场点的对称特性出发,推导获得了调和函数的法向偏导数在船体对称线元和面元处的积分对称关系式,开展了对称性在简化边界元积分方程系数矩阵中的应用,进而构建了高效求解三维船舶耐波性能的理论计算方法。数值结果及算例分析表明:水动力计算中移动脉动源格林函数的对称性有利于降低计算量和减少系统内存,且适用于有航速船舶在波浪中的波浪干扰力、运动响应和自由面波形的数值模拟,可拓展至两船水动力干扰及复杂船型的频域水动力分析计算之中。     

基于RWG基函数的电场积分方程的奇异性处理

采用矩量法分析导体三维散射体时,基于RWG基函数的电场积分方程存在奇异性,如果直接使用数值积分,则准确性很低。为了得到准确的积分结果,将被积函数拆分为2部分,对于无奇异点的部分直接使用数值积分求解,而对于包含奇异点的部分通过积分变换简化被积函数,得到解析表达式,计算实例验证了这种方法的正确性。     

有限水深Kelvin源格林函数及其导数的快速计算

应用线性船舶水动力学理论,将有限水深Kelvin移动兴波源格林函数分解成3部分:简单Rankine源集合、局部扰动项和波函数项。对各部分在亚临界和超临界航速时的快速计算分别作了讨论,其中对计算费时的局部扰动项采用了复数中的最陡下降积分法,加快了计算速度,进而引出各部分偏导数的计算。最后,比较了格林函数的差分值与偏导数值,计算了薄船兴波阻力与表面兴波,讨论了不同面元网格划分的收敛性,证实了亚临界、超临界航速兴波分别以横波、散波为主。     

一种计算内外域声场的简单源方法

使用赫尔霍兹积分方程计算声压的一个缺陷是积分变量过多导致的计算量过大."简单源"法将内域和外域的赫尔霍兹积分方程相互耦合,消去了多余积分变量.然而经典简单源法只能用于静止介质.为此,利用格林函数、势理论和运动介质波动方程,得到了求均匀流动介质条件下声压辐射的"简单源"边界积分式.与利用变换域中的赫尔霍兹积分方程求运动介质条件下声压的方法相比,新方法边界积分方程更加简洁和易于数值求解.基于新方法推导了有限长导管声场的预报模型.通过与标准问题的比较验证了新方法的有效性.     

再生核的一种新的计算法及其递推性

利用函数带积分余项的Taylor公式很自然地给出了W2m[a,b]空间的内积;基于这个内积,用Green函数得到再生核简洁的表达式,并用矩阵讨论了再生核计算的递推关系。     

势流理论的一种数值方法在空泡形态中的应用

基于势流理论的假设,采取常规边界元法对具有自由面的定常平板和轴对称空泡流动数值方法进行了研究,运用Riabourichinsky镜像空化模型,建立了一种空泡自由面求解的优化迭代方法,以自由流线长度为参数求解模拟了定常无重工况下流经二维平板和轴对称圆盘的空化流动.数值计算结果与解析解的相对误差不超过3%,证实了轴对称体空泡外形的近似回转椭球体假说在空化器附近是不成立的.计算结果与同类文献数据的对比显示了边界元法进行空化流场计算具有着显著的优越性,所采用的迭代方法精度高、收敛速度快,所得结论对深入了解空化特性具有参考价值.     

水下航行体泵喷推进器设计与性能分析

为了快速、高效地设计出性能优良的泵喷推进器,采用全三维反问题计算和正问题计算相互迭代的方法对某水下航行体的泵喷推进器进行了设计。在不考虑流体压缩性和粘性条件下,将泵喷推进器内部有旋流动分解为周向平均流动和周向脉动流动,其中周向平均流动方程通过三维欧拉方程在周向平均后导出,周向脉动流动议程通过Clebsch变换得到控制方程;不考虑叶片厚度的情况下,用叶片中心面上的涡来代替叶片对水的作用,叶片形状根据边界条件设置迭代确定。然后,运用雷诺时均方法对所设计的泵喷推进器进行数值仿真,发现泵喷推进器在设计航速下水力效率为89%,泵喷转子叶片表面未见空化,推进器产生的推力略高于航行体阻力。该结果表明:所提方法能满足设计要求,由此验证了设计方法的有效性。     

舰船偏离航道中心运动时水压场的计算方法

本文给出了舰船以亚临界航速偏离航道中心运动时水压场的计算方法。在船体细长、水深较浅的情况下,通过渐近匹配展开法可以得到舰船运动在矩形航道中的扰动速度势控制方程。用纵向分布源涡代替舰船的存在,为满足边界条件需要使用映像原理。从文献知道,源强正比于舰船横截面面积的变化率,涡强需从一个积分方程求取。通过数值计算得到源涡强度分布后,即可求出航道中舰船的扰动压力分布。     

一种机载无源定位方法——干涉仪定位

通过分析无源定位相对于有源定位的优点,给出了干涉仪机载无源测向方法.在解释一维、二维、三维相位干涉仪无源测向原理的基础上,分析了干涉仪无源测向方法的误差,并分别采用多基线干涉仪技术和粗测向设备解决了干涉仪无源测向方法存在的主要问题-相位测量模糊,解决了测向精度和有限的不模糊视角之间的矛盾,即保证了高的测向精度,又能无模糊地测向.     

三维边界积分方程中的奇异积分的处理

针对三维边界积分方程中奇异积分的不同形式分别进行了讨论,用刚体位移原理消除强奇异因子,用线性位移函数消除弱奇异因子,在计算工作量大大减少的情况下,保证了原求解精度的要求,并通过数值算例表明了该方法的可行性。     

压电材料中的共线非等长多裂纹问题

建立了横观各向同性压电材料中的共线非等长多裂纹模型。将裂纹模拟为连续分布的广义位错,推导了单点广义位错的格林函数,并基于此得到了共线非等长多裂纹的奇异积分方程组,最后由数值计算给出了裂纹尖端的能量释放率。对奇异积分方程组的退化分析和对称条件下的数值结果的讨论,验证了本文结果的正确性。研究表明:当共线裂纹间距小于裂纹长度时,裂纹尖端场存在显著的干涉效应;裂纹越长,干涉效应越显著;电载荷系数越大,能量释放率也越大;相对于外侧尖端而言,共线裂纹内侧尖端对电载荷系数更为敏感。     

有限水深Kelvin源格林函数的传播特性及其应用

为探究有限水深环境下自由面兴波在远场的传播特性,应用线性水动力学的基本理论,从有限水深Kelvin源格林函数的积分表达式中提取相函数,采用同相分析法求得亚临界和超临界速度区下远场传播波系的等相线。结果表明,水深傅汝德数Fh1时远场传播波包括横波和散波两个波系,F_h1时远场波系仅由散波构成,点源兴波的影响范围可由Kelvin角的大小确定,当F_h趋近于1时Kelvin角迅速增大至90°。依据船型要素、两船相对位置和池壁几何参数,并计及船模的定常兴波经池壁反射作用后对自身和另一船模的干扰,提出了有限水深条件下并行航行两船模池壁效应的判别方法,并结合实例检查分析了某试验工况下的池壁效应。     

含腔体导弹雷达截面高频计算建模

在建立了简化含腔体导弹外形模型的基础上,运用曲面像素法、数值积分法、等效电磁流法等多种高频RCS分析方法,采用目标部件分解法计算导弹各散射源的RCS,根据目标各部分散射场间的相对相位关系,计算出导弹整体RCS,这种理论建模方法简单,物理概念清晰,计算量不大,其计算精度虽不如数值计算法、矩量法、时域差分法等精确算法,但结果可为工程设计和导弹探测技术研究提供部分理论依据.     

声场波数积分最大截止波数自动选取算法

为了提高声场波数积分传递函数矩阵法的计算稳定性与精确性，提出基于预估-校正思想的最大截止波数自动选取算法，具有简单可靠、附带计算量小的优点。自由空间球面波算例测试结果表明，传递函数自下而上（单向）求解顺序的计算稳定性相对较差，在声源上方易出现积分核函数非物理发散，需要在不同深度上采用不同的截止波数才能计算出正确结果；而上下对进（双向）求解顺序的计算稳定性与计算精度均较高。     

高空滑移条件对化学非平衡粘性流动的影响

本文用数值方法计算了不同高度滑移边界条件对钝体绕流流场特性的影响。采用化学非平衡粘性激波层控制方程,边界条件含壁面滑移、激波滑移和无滑移条件。计算结果表明,考虑滑移条件后,激波层较厚,激波处的电子数密度为一有限值。文中给出了计算结果,并与其它文献的结果作了比较。     

由路面谱重构路面不平度的AR模型法

将等级路面的不平度视为一零均值的噪声序列通过一个AR模型产生的随机序列,通过求解Yule-Walker方程得出AR模型参数,进而得到路面不平度。经数值算例构造出4个等级的路面不平度,验证了用AR模型构造路面不平度方法的可行性与可靠性。     

IDE集成GCC编译器的方法

为了缩短基于Windows平台的嵌入式软件IDE的研制周期,利用自由源码的编译器GNU GCC,实现IDE的编译功能。通过借助于Cygwin环境,以及后台编译命令执行进程到前台界面显示进程的单管道设计,实现GNUGCC从Linux到Windows的移植,以及编译器编译信息的实时显示,从而实现GCC编译器在IDE中的集成。该方法不仅降低IDE的研制难度,也提高其编译功能的可靠性,为IDE的编译器集成提供一条有效的捷径。     

Warranty costs: An age‐dependent failure/repair model

An age‐dependent repair model is proposed. The notion of the “calendar age” of the product and the degree of repair are used to define the virtual age of the product. The virtual failure rate function and the virtual hazard function related to the lifetime of the product are discussed. Under a nonhomogeneous Poisson process scenario the expected warranty costs for repairable products associated with linear pro‐rata, nonrenewing free replacement and renewing free replacement warranties are evaluated. Illustration of the results is given by numerical and graphical examples. © 2004 Wiley Periodicals, Inc. Naval Research Logistics, 2004.